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ABSTRACT 


Accurate  nomographs  and  computer  programs  for  designing  conventional  skin-stringer  and 
honeycomb  sandwich  construction  to  meet  sonic  fatigue  requirements  are  presented.  The 
nomographs  were  derived  from  the  results  of  sonic  fatigue  tests  of  30  skin-stringer  and  30 
honeycomb  sandwich  designs,  totaling  60  test  specimens  for  each  design.  Emphasis  is  also 
placed  on  assessing  the  effect  of  structural  curvature  on  sonic  fatigue  and  structural 
response.  Nomographs  for  determining  the  curvature  effects  are  included. 
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dimensionless  frequency  parameter  defined  by  Equation  (95) 
defined  by  Equation  (58) 
defined  by  Equation  (58) 

Poisson's  ratio 

in-plane  restraint  factor  -  Section  F.2.c 

generalized  coordinate 

3.141592654  (approx.) 

mass  per  unit  area  of  a  plate  or  shell 
mass  per  unit  volume  -  Section  F.2.c 

defined  by  Equation  (58) 

mass  per  unit  area  of  a  plate  inertia  element 

mass  per  unit  length  of  rib  inertia  element 

mass  per  unit  length  of  stringer  inertia  element 

denotes  bending  stress 

defined  by  Equation  (58) 

th  '  *  th 

stress  in  the  r  layer  of  a  honeycomb  plats  in  the  i  direction 

shearing  stress 

(GJ/L).  -  Section  D.2.a 

eigenvectors  of  the  mass  matrix  -  Section  D.4 

th 

generalized  coordinate  for  the  (m,n)  mode 

mean  square  pressure  density 

variable  of  integration  -  Section  C.l  .a 
subtended  angle,  a/r.  Section  F. 2. a 
coordinate  direction  -  Section  F.3.a 

defined  by  Equation  (58) 

generalized  coordinate  for  (m,n)*^  mode 


coordinate  direction 


dimensionless  frequency  parameter  defined  by  Equations  (121)  and  (132) 

undamped  natural  circular  frequency  of  the  fundamental  mode 

th 

undamped  natural  circular  frequency  of  the  i  n  mode 
undamped  natural  circular  frequency  of  the  (m,n)^’  mode 

Notation 

rectangular  matrix 
diagonal  matrix 
transpose  of  a  matrix 
inverse  of  a  matrix 

over-dot  denotes  differentiation  with  respect  to  time 

2 

d  (  )/3x9y,  etc. 

denotes  differentiation  with  respect  to  the  variable  « 

over-bar  denotes  a  mean  value  except  in  Section  D  where  it  denotes 
system  coordinates 

double  over-bar  denotes  stiffness  of  system 
evaluated  at  x 

Subscripts 

denotes  curved 
denotes  edge 
denotes  flat 
denotes  honeycomb 
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I  -  INTRODUCTION 


Previous  programs  directed  toward  eliminating  or  reducing  sonic  fatigue  effects  from  intense 
acoustical  noise  on  flight  vehicles  have  established  design  criteria  and  nomographs  (presently 
available  in  Reference  1).  The  nomographs  in  Reference  1  were  empirically  derived  from 
sonic  fatigue  tests  of  certain  of  the  more  popular  structural  designs.  The  tests  were  made 
generally  using  discrete  frequency  excitation,  and  the  resulting  fatigue  data  were  convert¬ 
ed  to  "random  fatigue  data"  through  the  use  of  Miles'  single-degree-of-freedom  theoiw 
(Reference  2)  and  the  Miner-Palmgren  cumulative  damage  rule(Reference  3  and  4).  The 
nomographs  have  been  helpful  in  designing  sonic-fatigue-resistant  structure  for  present-day 
and  near-future  aircraft. 


However,  lightweight,  sonic-fatigue-resistant  flight  vehicle  structure  is  essential  to  high- 
performance  aircraft.  To  meet  these  critical  demands  for  weight  reduction  and  increased 
resistance  to  sonic  fatigue,  existing  design  nomographs  must  be  refined.  These  design 
nomographs,  based  on  limited  empirical  data,  have  proven  to  be  conservative  and  somewhat 
restricted  in  application.  It  is  essential  that  the  reliability  of  the  available  design  nomo¬ 
graphs  be  refined  and  extended  in  application  by  conducting  additional  sonic-fatigue  tests 
using  random  noise  excitation  representative  of  power-plant  noise  environment.  Initially, 
the  refinement  program  would  naturally  involve  the  more  commonly  used  structural  config¬ 
urations:  skin-stringer  and  honeycomb  sandwich. 

The  first  step  in  the  sonic-fatigue  refinement  program  was  accomplished  by  the  research 
program  reported  here.  A  description  of  the  theoretical  and  analytical  studies  is  given  in 
Section  II;  the  experimental  investigations  are  described  in  Section  ill;  Section  IV  is  a 
description  of  the  correlation  of  the  test  results  and  theory;  development  of  the  design  nomo¬ 
graphs  is  given  in  Section  V;  and  a  description  of  the  test  facility,  data  collection  and  data 
reduction  system  and  test  procedure  are  described  in  detail  in  Appendices  I  and  II. 
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II  -  THEORETICAL  AND  ANALYTICAL 
A.  Introduction 

This  research  program  was  primarily  experimental  in  nature.  Consequently,  one  of  the 
main  objectives  of  the  theoretical  effort  was  to  suggest  the  significant  variables  and  the 
mathematical  expressions  governing  skin-stringer  and  honeycomb  sandwich  structural 
response  to  acoustical  excitation  representative  of  in-service  environment.  Also,  the 
utility  of  the  design  nomographs,  very  important  end  products  of  this  program,  depends 
upon  simplicity,  ease  of  use,  and  accurate  representations  of  the  test  results  from  which 
they  were  derived.  Therefore,  the  classical  theory  for  the  complete  description  of  the 
response  of  a  continuous  structure  such  as  an  aircraft  fuselage,  wing  flap,  or  control 
surface  as  described  in  Reference  5  is  too  complex  for  design  purposes.  Also,  it  will  be 
shown  in  Section  III  that  the  majority  of  the  structural  response  results  from  one 
predominant  mode,  usually  the  fundamental. 


B.  General  Theory 


From  Reference  6  ,  it  is  shown  that  the  mean-square  displacement  for  single-mode  response 
can  be  expressed  as  (for  the  i^  mode) 


irF.^(xj) 
2u  3  r). 


i  'll 


//  Fi<xA)  F;(*B)  Gp(xA'xB;  "PdAdA 

A  A 


(1) 


where 


Fj(x) 
M 

CJ 

V 


Gp(xA/XB;u>;) 


is  the  i*^  mode  shape 
is  the  generalized  mass 
is  the  undamped  circular  frequency 
is  the  structural  damping  factor 

is  the  cross-spectral  density  function  of  the 
excitation  pressure  at  x^  and  xg 

is  the  predominant  mode,  usually  the  fundamental 
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Equation  (1)  describes  the  structure  as  a  set  of  independent  plates  or  panels  responding  in 
their  first  mode  only. 


To  simplify  Equation  (1),  it  will  be  assumed  that  the  excitation  pressure  is  exactly  in  phase 
over  the  whole  panel.  The  mean-square  deflection  response  is  then  expressed  as 


— «  n  F.  (x,) 

V\T(t)  '  1 


2  <u.3  tj.  M.2 

i  'i  i 


f  F!(xA)dA 
A 


G  (c o.) 
P  ' 


(2) 


The  static  displacement,  W  ,  at  Xj  due  to  a  uniform  static  pressure  of  unit  magnitude  is 


W 

o 


F.(x,) 


(O 


,2m. 


/  F;<-A>dA 
A 


(3) 


2 


and  Equation  (2)  can  be  expressed  as 


W2(t)  =  <o.  G  (*>.)  W  2  (4) 

'  ~  i  p  t'  o  '  ' 

2  t;  .  r 

The  expression  commonly  used,  however,  for  the  mean-square  stress  in  terms  of  the  viscous 
damping  ratio,  5  ,  and  the  frequency,  f.,  in  Hz.  is 

t)  =  -L  U  (0  o  2  (5) 

45  1  ° 

where  oQ  is  the  stress  at  the  point  of  concern  due  to  a  uniform  static  pressure  of  unit  mag¬ 
nitude,  and  <|>  (f)  is  the  mean-square  pressure  density  function. 

C.  Skin-Stringer  Construction 

This  subsection  is  a  description  of  the  analysis  defining  the  response  of  a  single  bay  of  skin- 
stringer  construction  to  an  arbitrary  acoustical  excitation. 

Most  panels  have  edge  conditions  which,  in  practice,  are  between  simply  supported  edges 
and  clamped  edges .  The  basic  approach  will  be  to  determine  panel  response  for  both  clamped 
edges  and  simply  supported  edges.  These  two  edge  conditions  will  provide  a  bound  in  the 
actual  situation.  The  overall  motion  of  the  supporting  structure  will  be  presented  in  the 
next  subsection. 

1  .  Frequency  Analysis 

The  Rayleigh  energy  method  will  be  used  below  to  obtain  the  response  of  a  panel  vibrating 
in  its  fundamental  mode.  The  panel  material  is  considered  to  be  Hookean  and  isotropic, 
and  the  thickness  remains  constant  before,  during,  and  after  deformation. 

a.  Clamped  Edges  -  The  strain  energy  expression  for  a  thin  elastic  plate  of  length  a  and 
width  b  is 


/ 


Vw 

ax  ay 


o  o 
v  2 


^  +  A  1  -2(1-0 

3x2  ay' 


3 2w  d  2w 
3x2  3v2 


dydx 


/ 


and  the  kinetic  energy  for  the  vibrating  plate  is 

T  =  2  /  /  P  (7f)  dxdy 


(6) 


(7) 


where 

P  is  the  mass  per  unit  area 
D  -  Et3/1 2(1  -  v2) 

E  is  Young's  Modulus 

t  is  the  panel  thickness 

v  is  Poission's  ratio 

w  is  the  transverse  displacement 

For  clamped  edges  it  is  assumed  that  the  mode  shape  is  as  shown  in  Figure  1 . 


FIGURE  1 .  ASSUMED  MODE  SHAPE  FOR  A  SKIN-STRINGER  PANEL 
WITH  CLAMPED  EDGES 

The  first  mode  response  is  represented  mathematically  as 

w  (x,  y,  t )  =  [  1  -  cos  (2  rr  x/a)  ]  j  1  -  cos  (2n  y/b)  J  £  (t)  (8) 

When  Equation  (8)  is  substituted  into  Equations  (6)  and  (7),  the  expressions  for  the  kinetic 
energy  and  the  strain  energy  becomes 
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Basically,  the  kinetic  energy  and  strain  energy  are  expressed  as 


T  =  - 
2 

m  i  2  (t) 

(15) 

and 

u  =1 

2 

CM 

(16) 

where  M  is  the  modal  mass 

Kis  the  modal  stiffness  constant. 

The  modal  mass  constant  and  the  modal  stiffness  constant  from  Equations  (13)  through  (16) 
are  given  by 
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(17) 


M  —  — —  p  a  b 

4 


The  natural  circular  frequency  of  the  fundamental  mode  is 


co2  = 


and  from  Equations  (17)  and  (18): 

^  "4o  ,t 

9p  2,2 

a  b 

or  in  terms  of  the  frequency  In  Hz: 


(20a) 


f2  =  ±JL 


9aV  P 


(20b) 


b.  Simply  Supported  Edges  -  Using  the  same  approach  as  in  the  previous  subsection,  the 
mode  assumed  for  "simple  supports  is  of  the  form 


w  =  £  (t)  sin  (■  n—  )  sin  (-JLX  ) 
a  b 


The  modal  mass,  modal  stiffness  constant,  and  first  mode  frequency  are 


M  =  -  P  ab 
4 


1  +  _1 


co  -  n 


4  d  r  i 


(24a) 


JL  JL 

4  "p~ 


(24b) 
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2.  Stress  Response 

The  stress  analysis  of  a  single  bay  of  a  skin-stringer  panel  is  based  on  the  assumption  that 
the  elastic  characteristics  of  the  stiffeners  produce  negligible  effect  on  the  stress. 

a.  Clamped  Edges  -  The  expression  for  the  stress,  due  to  bending  at  the  surface  of  a  thin 
plate,  is  given  Ey  (see  Reference  7  ). 

M 

°x  -  -i-  (25) 

z 

where  z  is  the  section  modulus  and 


M  =  -  D 
x 


2  2 

3  w  .  .  aw 

— T-  +  — r- 

d  k  d  y2 


The  fundamental  mode  can  be  defined  as 


w  =  Wo  (l  -  cos  (JLLJL)}  (  1  -  cos  (  . /-  ))  (27) 

where  W  is  the  static  displacement  at  the  center  of  the  bay  due  to  a  uniform  static  pressure 
of  unit  n?cgnitude. 

After  substituting  for  w  in  Equation  (26),  the  bending  moment,  M  ,  is  given  as 

=  4  it  2DWq  — y  (cos  (-"-))(l  -  cos  (^~-^))  +  m ^(cos  (~£  ty)(l  -  cos  (—■-))  (28) 
A  a  b 

Substitution  of  Equation  (28)  into  Equation  (25)  yields  the  following  expression  for  the  stress: 

2  r 

w  D  ...  1  2rrx  ,,  2rry\ 

a  =  24  9 —  W  -y  cos  —  (1  -  cos  — r--2-) 


JUS) 
a  ' 


o  2 
a 


2?r  y  , ,  2 

cos  — j^-*-  (  1  -  cos  — 


The  point  of  interest  on  the  panel  is  the  center  of  the  long  side  along  the  fastener  row, 
where  the  moment  is  maximum,  i.e.,  at  x  =  0,  y  =  b/2,  and  x  =  a,  y  =  b/2.  At  this  point 
the  stress  is  given  by 

48  n2 D  W 


W  is  evaluated  for  this  case  as 

o 

W  -  F/K 
o 
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where  F  is  the  generalized  acoustic  pressure  force 
a  b 

F  =  f  f  PQ  (1  -  cos  (^!p))  (1  -  cos  dx  dy 

o  o 

and  Pq  is  the  acoustic  pressure  amplitude.  By  integrating  Equation  (32) 

F  =  P  ab 
o 


after  substituting  Equations  (18)  and  (33)  into  Equation  (31),  W  is  evaluated  as 
°  4„4R  \D  / 


By  combining  Equation  (30)  and  Equation  (34),  the  static  stress  due  to  a  uniform  static  pres¬ 
sure  of  unit  magnitude  is  obtained: 

—  ^  ^  /  b  1  /nrt 

(T)  r  (35) 

max  n 

An  estimate  of  the  mean-square  stress  response  of  a  lightly  damped  panel  exposed  to  random 
noise  excitation  can  be  obtained  by  substituting  Equation  (35)  into  Equation  (5) 


-2  =  36.  ,b*4 


b.  Simply  Supported  Edges  -  Assuming  that  the  elastic  characteristics  of  the  supporting  struc¬ 
ture  may  be  approximated  by  simply  supported  edges,  the  mean-square  stress  in  a  bay  of  skin- 
stringer  construction  is  developed  below,  using  the  same  general  procedure  as  described  in 
Section  2. a . 

The  fundamental  mode  for  a  simply  supported  panel  is  described  as 

w  =  Wq  sin  (-— )  sin  (-|^)  (37) 

Equations  (37)  and  (26)  are  used  to  evaluate  the  bending  moment.  Substitution  of  this  value 
into  Equation  (25)  results  in  the  expression  for  the  stress.  The  maximum  stress,  of  course, 
occurs  at  the  center  of  the  panel  and  is 


9  DW 

a  =6  n  — S- 

Xmax  t 


1  ,  „  1 
— 7T  +  V  — 7 
a  I  i 

a  b 


8 


The  static  deflection  at  the  center  of  the  panel  due  to  a  uniform  static  pressure  of  unit  mag¬ 
nitude  is 


W 


«  F/K 


where  F  = 


a  b 

J  f  s!n  IT sin  dy  dx 

o  o 


4ab 
v  2 


(39) 


K  is  given  by  Equation  (23). 


The  static  deflection  at  the  center  of  the  simply  supported  panel  due  to  a  unit  pressure  be¬ 
comes 


W 
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16 


rr6  D 


1 

“2 
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(40) 


The  maximum  stress  due  to  the  uniform  static  pressure  is 


a 

x 

max 


96 


42 

77  t 


The  mean-square  stress  response  of  a  lightly  damped  bay  of  skin-stringer  construction  consid¬ 
ered  to  be  simply  supported  by  the  stiffeners  and  excited  by  a  random  acoustical  noise  is  ob¬ 
tained  by  substituting  a  from  Equation  (41)  for  a  in  Equation  (5). 

X  o 

max 

Hence,  the  mean-square  stress  response  of  a  simply  supported  bay  of  skin-stringer  construc¬ 
tion  is  given  as 


_2  _  2304 


J_  +  JL 

2,2 

a  b  . 


i2 
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max 
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(42) 


D.  Vibration  Analysis  of  Supporting  Structure  using  Finite  Element  Techniques 
1  •  General  Theory 


The  analytical  model  of  the  supporting  structure  was  assumed  to  be  a  thin  cover  sheet  sup¬ 
ported  by  two  parallel  ribs  and,  perpendicular  to  the  ribs,  two  parallel  stringers,  with  the 
edges  of  the  system  completely  restrained  from  movement.  This  choice  was  prompted  by  the 
design  of  the  test  panels  used  in  the  experimental  phase  of  the  program.  The  coordinate 
system  and  general  arrangement  for  the  stiffened  panel  is  illustrated  in  Figure  2.  The  inter¬ 
section  points  of  the  ribs  and  stringers  were  selected  as  grid  points  for  measuring  the  defor¬ 
mation  of  the  system.  Each  grid  point  has  3  degrees  of  freedom:  two  rotations  and  a  trans¬ 
lation.  These  coordinates  were  ordered  as  shown  in  Figure  2. 
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The  equations  of  motion  represent  a  12 -degree -of- freedom,  spring-mass  system.  These  equa¬ 
tions  are  of  the  form 


I  h  f  I 
xx  j  xz  j  xy 

Cfc  i  w 

rTiir  ■ 

y«j  YY 


K  I  K 

XX  ,  xz 


K  I  K 

ZX  I  zz 


'  r K  ~^K 

yx  j  yz  j  yy 


for  free  vibrations.  The  above  equations  are  the  classical  eigenvalue  problem  of  linear 
algebra.  Details  of  the  formulation  of  the  mass  and  stiffness  characteristics  and  the  solution 
of  the  above  problem  are  given  in  the  following  sections. 


2.  Formulation  of  the  Stiffness  Matrix 


a.  Supporting  Structure  -  The  stiffness  matrix  for  the  system  will  be  treated  in  two  parts: 
supporting  structure  and  cover  plates.  The  supporting  structure  stiffness  matrix  will  be  syn¬ 
thesized  as  a  network  of  beam  elements.  The  stiffness  characteristics  of  the  system  are  de¬ 
veloped  using  the  equilibrium  method.  This  method  is  well  described  in  Reference  8  . 

Beginning  with  the  beam  element  shown  in  Figure  3,  the  load-deflection  relation  for  the 
ith  element  is 

P1  1  K,  I  KI2  1  [d.  1 

-  -  - 1 - -  (44) 

P2  .  K21  I  K22  .  d2  . 


where  at  station  j  on  the  element 
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For  simplicity,  it  was  assumed  here  that  there  was  no  coupling  between  the  bending-torsion 
coordinates  of  the  beam  element. 

The  elements  were  ordered  as  shown  in  Figure  4  ,  and  it  is  necessary  to  rotate  the  coordinates 
of  the  elements  to  match  the  system  coordinates.  Such  a  transformation  is  easily  accomplished, 
and  is  represented  as 

T 


K. 


K. 


Ir 


I  L. 


(46) 


for  the  i^  member,  where  the  overbar  denotes  the  element  stiffness  in  system  coordinates 
and  the  transformation  matrix  [  T  ]  ,  is 


cos  a . 

i 

0 

sin  a . 


0 

1 
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-sin  a. 

i 


cos  a. 

i 


(47) 


for  the  coordinate  transformation,  as  illustrated  in  Figure  5. 

As  shown  in  Figure  4,  the  elements  are  oriented  by  the  direction  of  an  arrow  which  points 
from  station  1  to  station  2  for  each  element.  Equations  (48a)  through  (48d)  are  the  compati¬ 
bility  equations  for  grid  points  A  through  D,  respectively  (system  coordinates  are  denoted 
by  an  overbar): 


1 

JA 


d 

fd 


d2. 

2 

v 

3 

d2_ 

d2 

3 

d2 

5 

_dl  . 

d2 

1 

d2 

9 

d2 

d2 

6 

d2_ 

8 

_dl_ 

2 

6r 


12 


i_dij 


12 


10 


(48a) 

(48b) 

(48c) 

(48d) 


The  boundary  conditions  for  clamped  edges  are  given  by 


J  r 


(49) 


for  r  =  1,  2,  4,  5,  7,  8,  10,  11 
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For  equilibrium  of  forces  at  each  grid  point  the  following  equations  are  obtained: 


(50a) 

(50b) 

(50c) 

(50d) 


where  the  load-deflection  relations  for  the  i^  element  in  the  system  coordinates  are 

P, 
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-J  i 
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(51) 


After  combining  the  compatibility  equations  (48a-d),  the  boundary  conditions  (49),  and  the 
load-deflection  equation  (51),  the  equilibrium  equations  (50a-d)  become 


where 
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After  grouoinc)  the  coordinates  as  implied  by  Equation  (43)  where 
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and  ufrer  dropping  fk*  overbad,  the  stiffness  matrix  for  the  supporting  s‘tychw  becomes 
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<6[(/3L)6- 

<6(0L)I2) 

(0) 

(0) 

(6(/3L)6) 

(0)  (-608  l)12)  (0) 


iei)5})  (o)  (-6(0D6) 

(6[OL)n  -(/3L)|2j)  10) 

(0)  (6|OL)]0-(/3L)6j) 


(55c) 
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K 


l 


K 

yy 


<t4  +  tI2  +  4[(*l2)2+  0^)3j)  (?OL2)3)  (-t  12)  (0) 

(2{/St.2)3)  (T5  +  t6+4[(/9L2)3»(j9L2)7])  (0)  C"t6) 

(-*,2)  (0)  (tn  +  tu  +  4[(/3L*?)j  +OL2)9J)  (2(^L2)9) 

|^(0)  (-t6)  (2(/3L2)9)  (^6+  T10  +  4f(^L2>8+(^L2)9]) 

(6[OL)3-()9L)2])  (-6(6L/3)  (0)  (0) 

(6(/3L)3)  (6[(/3L)7-(/3L)3])  (0)  (0) 

(0)  (0)  (6[<0L)9-(/3L),j)  (-6(/3L)9) 

!  (0)  (0)  (6(/SL)9  (6[OL)8-OL)9]) 

(I2[/32  +/33  +/34  +  /3]2j)  (-I2J83)  (-12/5,2)  (.0) 

(-12/33)  (J2[^3+/95  +  06  +  /®7])  (0)  (-12  fi6) 

(-120, 2)  (0)  (12[0, +09  +  0,,  +  0,2J)  (-I2/59) 

(0)  (-I2/36)  < -12/39)  (12[/36  +  /^8  +  ^9  +  ^]0]) 


(55e  ) 


(550 


and  the  symmetry  condition 


K.. 

= 

K.. 

>0 

(55S ) 


b.  Cover  Plates  -  The  approached  used  for  deriving  the  stiffness  matrix  of  a  cover  plate 
element  is  given  in  Reference  9,  and  a  brief  derivation  is  given  in  the  appendix  for  com¬ 
pleteness.  For  the  coordinates  in  Equation  (43)  and  the  nomenclature  ir>  Figure  5,  the 
stiffness  matrix  for  the  complete  cover  sheet  is 
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(Vjaj  * +«4  +a5j)  (0)  (a4+o5)  (0) 

(0)  (2|«2 +a3  1  «5 +a6j  )  (0)  (a5+a6) 

(0)  (2[°4  +  fl5  +a7  'aQ  j)  (0) 

(G)  <«5+a6)  (0)  (2[«5 +  a6 +«*8 +«9|) 


(  56a ) 


(P,  +P2-P4-P5)  (30[(S/b)5  -(  «/b)2])(cr4  +  <r5)  (  -6(  5/b  )4  ) 

(30[(5/b)5  -  (5/b2  )  ]  (  P2  +  P3  -  )  (-6(5  A  )5  )  (ct5  +  ct6) 

-<ff4  +  <75)  (6(5A)5)  (p4+  <Vp7-p8)  (30[(.5A)8-(sA)5J) 
(6(5A)5)  -  (<r5  +  a6)  (30[(SA)g  -  <*/b)5])  (  P5  +  P 6  -  P g  -  P ?) 


(56b) 


(203  +  +  04+/|])  (/32  +  ^5)  (0)  (0) 

(/32  +  ^5)  (2[/32  +  ^3^5  +  i36])  (0)  (0) 

(0)  (0)  <2[VVV08]) 

(0)  (O'  +  , S8)  (2[/85+/36  +  ^8+/39]) 


(56c  ) 
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(  *2  +  *  5  ~  *  1  ~  ^  4  ^  ^-T2  ~T  5^  (  30 1  ( 5/ a  -  (  S/a)5j  )  (6(5/a)^) 

(  r  2  +  T  5  )  (X3+X6-X  2~*  5>  (-6(5/°)5)  (30[(5/a)5-(5/a)6J) 

(30[(S/a)4-(S/a)5  ])(6(S/a)5)  (  X 5  +  X Q  -  X  4  -  X  ?)  (  -  t,_  -  t  g  ) 

(-6(S/o)5)  (30[(8/a)5-(6/a)6l)  (x5  +  cg)  <  *6  +  *9 - *5  ’  V 


(56d ) 


(  «,  +k2+ k4 +  k5)  (-X2-X5)  (-M4-/i5)  (y5  +  72(5 /ab)5 

(-X2-X5)  («2  +  *3  +  'W  (>/5  +  72(S/ab)5) 

(  y5  +  72(  5/ab  )5 )  <  *4  +  *5  +  *7  +  *g  )  ("X5-A8) 

(V5  +  72(S/ab)5)  (-a5-a8)  (K5+'t6+K8+K9) 


(56e ) 


where  for  the  i 


.th 


panel  element 


D.  = 

i 


(57) 


P.  =  D.(a/b).  5.  =  — (  wD). 

i  ii  1  16 

P.  =  3  (a/a).  +  10(  5/b ). 

a .  =  3  (a/a).  +  2(5  /b  ). 

*;  =  3  (/3/b ).  +  10(  5/a ). 
t.  =  3  (/3/b).  +  2(  5/a  ). 

«.  =  y.  +  7? (5/ab  ).  +6(«/a2).  +6  (/3/b2). 

A.  =  y.  +  72  (  5/ab  ).  +  6  (  £/b2). 

^  .  =  y.  +  72  (  5/ab  ).  +  6  (a /a2  ). 


(58) 


3.  Inertia  Matrix 

The  Inertia  characteristics  of  the  system  is  derived  below  by  considering  the  rigid  body 
motion  of  a  patch  of  structure  as  shown  in  Figure  6.  The  displacement  of  a  point  in  the 
surface  of  the  patch  at  a  point  (*,z  )  is 


(59) 
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The  kineric  energy  of  the  patch  of  structure  is 


a2  „  b2 


U  l 


PD  d  x  d  z  +  — 
K  2 


al  "bl 


+  1 


/ 


Ps^s 


d  x 


“bi 


/ 

-a, 


Vs 


dz 


(60) 


where 

y^(t)  =  8(t)-z9^(t)  -  displacement  of  rib 

y_(t)  =  8,(t')  +  x9(t)  -  displacement  of  stringer 

o  y  z 

p 

P  -  mass  per  unit  area  of  the  plate 

Pr  -  mass  per  unit  length  of  the  rib 
p 

S  ~  mass  per  unit  length  of  the  stringer. 

By  integrating  the  expression  for  the  kinetic  energy  one  obtains 

T  =  ;[(mp  +  mR  +  mS  )Sy2  *  ^  P1  <mp  +  mS  )0Z2  +  5  P2(mp  +  mR,6x2 


+  Po(m  +mc)5  9  —  P  j  ( m  +  mD)5  9  -  — m  PQ  P .  9  9 
o'  p  5  wy  z  4'  p  R  y  x  7  p  3  4  x  z 


where 


mp  =  Pp(^i  +  a2 ) (bj  +  b2 ) 


m 


R  Val  +a2) 

S  =  PS^bI  +b2^ 


mass  of  plate  element 
mass  of  rib  element 
mass  of  stringer  element 


(61  ) 


Pj  bj  ~  bj  h>2  +  b2 
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2  2 

P2  =  °1  “  °la2  +  °2 

P3  ~  b2  -  bl 


P4  =  a2  -  a, 


From  Equation  (61  ),  it  follows  that  the  inertia  forces  at  the  coordinates  are 

M  ~  —  P«  (m  +  mD  )*9  -  —  m  P~  P /O  -  —  P,  ( m  +  mD )  5 

x  3  2'  p  R  x  ^  p  3  4  z  2^  P  R  y 

Mz='^mpP3P4ex  +^Pl(mp  +  mS)0z  +"P3(mp  +  "S^’y 

Py  =-'-p4(mp+mR)0x  +^P3(mP  +  mS)ez  +(mp+mR  +  ms)5y 


Then,  with  reference  to  Equation  (43)  and  Figure  5,  the  inertia  matrix  is 


t 

r, 

1 

1 

XX 

-<P2(mp  +  mR,)i 

o 

L 

4 

c 

(m  P-iP^). 

4  P  3  4'. 

1 

xz 

= 

= 

0 

zx 
L 


xy 


-l(p4(mp  +  mR)). 


yx 
L  :J 


zz 
L  ^1 


-(Pi  (m  +  m,  )  ). 

3  i  p  oi 


ZY 


-( P^(m  +  m  )  ) 

2  J  P  °  1 


Yz 

il  L  Ol 


(62a ) 
(62b) 
(62c) 

(63a) 

(63b) 

(63c) 

(63d) 

(63e  ) 
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( m  + 
P 


mR  +  mS)i 


(63f) 


where  i  ”•  1,  2,  3,  4  denotes  grid  points  A,  B,  C,  and  D,  respectively. 

4.  Method  of  Solution 

The  general  method  of  solution  of  the  equations  of  motion  (43)  is  presented  below.  Begin¬ 
ning  with  the  equations  of  free  undamped  motion  of  a  spring- mass  system 


MX  +  k]  X  =  0 


one  obtains  the  classic  eigenvalue  problem  of  linear  algebra,  upon- assuming  simple  harmonic 
motion 


’x]  =  eiwf  [x  • 


The  relation  (65)  is  substituted  into  the  equations  of  motion  (64),  and  the  eigenvalue  problem 
is  obtained  in  the  form 


K  X  =  co  2  M  : 


From  a  computational  standpoint,  it  is  advantageous  to  solve  an  eigenvalue  problem  of  the 
form 


Q  Z  =  X  Z 


where  Q  j  is  a  real  symmetric  matrix  -  that  is,  it  is  advantageous  to  transform  Equation 
(66)  to  the  form  of  Equation  (67)  before  solving  the  eigenvalue  problem  (66)  directly. 


The  following  method  is  used  to  perform  the  above  mentioned  transformation.  The  eigen¬ 
values,  [M|  ,  and  the  eigenvectors, [  ]  ,  of  the  mass  matrix  are  computed.  The  mass 

matrix,  [Mj  is  then  expressed  in  terms  of  its  eigenvalues  and  eigenvectors  with  form 

r  i  r-  i  i  r  iT 

M  -  0  M  $  (68) 

LJLniJLJLmJ  ' 

M  from  Equation  (68)  is  substituted  Into  Equation  (66),  and  the  eigenvalue  problem  becomes 

k  i  [  x  =  o) 2  r <t>  ]  r m  i  r  iT  r  x  i  (69) 

LJLj  LmJL^ILmJLJ 

After  premultiplying  Equation  (69)  by  ,  letting 

r  i  r  -iT  r  ir  i 

H  -  <J>  K  $  (70) 

J  L  mj  L  JLm. 

and  making  the  transformation 

Y  ]  -  [  $mlT  [  X  ]  (71) 


25 


the  eigenvalue  problem  (69)  becomes 
H  Y  =  w2  [  Y 

Now,  by  expressing  the  eigenvalues,  £  m],  of  the  mass  matrix,  M  , 

[»]'  MM  r- 

where,  obviously,  G..  =  \  /  M..  ,  Equation  (72)  becomes 

["][*>  -2[«Wv] 

M'1  [h][g]''[g][y]-  "2[g][4 

Finally,  after  making  the  transformation 

[2HG][9] 

the  final  form  of  the  eigenvalue  problem  (64)  is 

[«r  [«][«]"  “2w 


The  eigenvalue  problem  (77)  is  solved  by  the  Jacobi  method,  and  the  resulting  eigenvectors 
are  transformed  back  to  the  original  problem  (Equation  64)  by  the  relation 

[*]-  [*-][«]■'  [4  •  (78) 

Of  course,  the  eigenvalues  of  Equation  (77)  are  invariant  under  the  above  transformations. 
It  should  be  noted  that 

m"  ["][«]"  ■  [°r  [•-]t[k]  <79) 

is  indeed,  a  symmetric  matrix.  A  more  detailed  discussion  of  the  above  theory  is  given  in 
Reference  10.  The  above  problem  is  programmed  for  the  IBM  Remote  Access  Computing 
System.  The  program  is  written  in  IBM  System/360  Basic  Fortran  IV  Language.  The  coordi¬ 
nates  are  numbered  as  shown  in  Figure  2  and  the  elements  are  numbered  as  shown  in  Figure 
4.  Use  is  made  of  the  IBM  furnished  subroutines  "LOC",  "EIGEN",  and  "NROOT".  The 
subroutines  are  described  in  detail  in  Reference  1 1  .  The  program  accepts  basic  input  data 
and  fills  the  inertia  and  stiffness  matrix  directly.  The  output  format  consists  of  natural  fre¬ 
quencies  and  the  corresponding  mode  shapes.  As  shown  in  Figure  7,  the  mode  shapes  are 


ordered  as 

0xA 

°xB 

°xC 

°xD 

°zA 

0zB 

9zC 

0zD 

SyA 

SyB 

5yC 

5yD 

rotation  about  the  x  axis  at  the  i^1  grid  point 
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VIBRATION  ANALYSIS  CF  A  ST  I F  F  ENEC  PANEL 
FRE  CUE  NC Y  «  0.11396C  C4  FZ. 

■PC CF.  SFAPE 

C.11095C  C0-C.11O95C  CC  C. 110951  0C-C.11C95D  CO 
-0.172470-01-0. 1 7247C-C 1  C  .  1 7  2~4  7  C  -  C  t  C. 172A70-C1 
C.48723C  00-0.48723C  00-0.4e723C  CO  0.4B723D  CC 
F  P  ECLENC  Y“  C.11372C  OA  FZ. 

FCCE  SHAPE 

C.1117CC  00  0.II17CC  CC  0.I117CC  CC  C.1117CD  CC 
-C.18050C-01  C.18C5CC-01  0 . 1 EC5 C C-C l-C . 1 BC5CD- C I 
C.48T03C  CC  C.A87C3D  00-C.A87C3C  C0-C.Ae7C3C  OO 
F  R  ECLE  NC  Y  *  0.69184C  C3  HZ. 

FCCE  SHAPE 

-C. 220190  00  C.22CJ9C  00  0.22C1SC  CC-0.22019C  CC 
-C.21795C-C1-C.21795D-01-C.21795D-C1-C.21795C-01 
0.AA837C  C0-0.AAe3eC  CC  C.44B37E  CC-C.AA838D  OC 

wEauTfrrv*  o.isvrnr  ciftt 

FCCE  SHAPE 

-C. 223280  CC-C. 223280  CC  C.2232eC  CO  C.22328C  00 
-0.237A6C-C1  0.23746C-C1-C.23746C-C1  C.237460-C1 
C.AA6750  CO  0.44674C  CC  0.44675E  CO  0.446740  CC 
FRECUENCY*  C.62562C  03  FZ . 

FCDE  SHAPE 

-C.16A72C  CO  0.164  72C  CC-C.16472C  CC  C .  It  A  720  CC 
-C.22C50C-01-C.22C5CC-C1  0.22C5CS-C1  C.22C5CD-C1 
C.A7157C  CC-C.A7157C  CC-C.A7J57C  CO  0.A7157C  CO 
FRECLE  NC  Y  =  C . t2 360C  03  HZ. 

FCCE  SHAPE 

-C.163A1D  C<P0".1634  1C  CC-C.16341E  CO-0.1t3AVC  CC 
-C.2A6AA0-C1  C  «  2 At  4  AC-0 1  C . 2 46 AA 0-0 1 -0 . 2A6 A AC-C 1 
C.A7190C  CC  O.A7190C  CC-C.A719CC  CC-C.A719CD  OC 
FRECUENCY*  0 . 3367  1C  C3  FZ. 

FCCE  SHAPE 

-c.8t2e6C-jn  0.8628^0^0  i_c.  86 2860-ci-r.E6 2 e6C-oi 
0i2  ie2A'C-0 1  0. 21 E24C-C 1  C  .  2  1 8  2  AO-C 1  C. 218240-01 
-C.A92C1C  00  0.492C1C  0C-0.A92C1C  CC  C.A9201D  CC 
F  RE  CUENC  Y  «  0.33CS1C  C3  FZ. 

FCCE  SHAPE 

0.e5A12C-01  0.85412C-CI-C.B5412C-C1-C.85412D-01 
-C.33122C-01  0.33122C-C1-0.33122C-C1  0.331220-C1 
C.  49 1540"' CO  C.A9T5ACCC  C.A915AC  00  0.49154D  CO 
FRECLENCY=  0.169510  C3  HZ. 

FCDE  SHAPE 

C.3C7A6C-02-C.3C7A7C-02  0 . 3C 75 E C- C 2-0 . 30 7 5 ED- C 2 
-0.A99720  CC-C.A9972C  OC  C.49S68C  CO  C.49968D  CC 
-0.17C19C-C1  C.17C19C-C1  C.  170 100-0 t-C .  1 7C 1C0-C1 
FREQUENCY*  0.16836C  C3  FZ « 

FCCE  SHAPE 

C. 12AA60-C1-C • 12AAtC-01-C . 12445D-01  0.12AA5D-C1 
0  A9C66C  CC  0.49067C  CC  C.A9071C  CC  C.A9C71D  CC 
C . 952A  1C-01— 0. 9  52A1C-C1  0 . 9 5 2  A 3C- C 1- C . 9 52 A2D- C 1 
FRECUENCY«  0.13569C  C3  FZ. 

FCDE  "SHAPE 

-0.8222AC-C2-0.8222AC-C2-C.8222AC-C2-C.e222AC-CZ 
C.A9676C  CO— 0.496  76C  CC-0.A9t7tC  CC  C.A9t7tO  CC 
C.  5625 30-01  C.56252C-01-0  .56  25  2C-01-0.562  5  30-01 
FREC'ENCY=  C.12830C  C3  FZ. 

FCOE  SHAPE^ _ 

-C.39885C-01-0.39FE5C-C1  0.39EE5C-C1  0.398650-0 
-C.A05570  CC  C.4C557C  0C-C.AC557C  CO  C.AC557C  CC 
-0.28970C  CC-C.2897CC  OC-C.2897CO  C0-C.2eS7CD  CO 


FIGURE  7.  TYPICAL  COMPUTER  OUTPUT 
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0  .  -  rotation  about  the  z  axis  at  the-  if^'  grid  point 

Sy.  -  deflection  in  the  y  direction  of  the  5^’  grid  point 

i  -  1,  2,  3,  4  denotes  A,  B,  C,  D,  respectively. 

The  mode  shapes  are  normalized  to  the  length  of  the  eigenvector  so  that 
4 

Z  (<£  +  ®2i  +  Sli  >  =  1*0  (80) 

i  =  1  1 


In  Figures  8a  through  81  the  mode  shapes  for  a  typical  stiffened  panel  are  sketched  with  the 
deformations  highly  exaggerated.  It  is  seen  that  three  modes  will  respond  to  normal  forcing 
ie.,  Figures  8d,  8h,  and  81.  These  modes  represent  rib-bending  stringer-torsion,  rib-tor¬ 
sion  stringer  bending,  and  coupled  rib-stringer-bending,  respectively.  It  is  difficult  tostate 
which  mode  will  be  the  most  responsive  to  excitation,  since  the  panel  modes  described  in 
Section  C  must  also  be  taken  into  account.  However,  the  analysis  does  explain  some  of  the 
coupling  effects  which  arise  due  to  the  substructure  motion.  If  ne  stiffeners  are  uniform ar.d 
equally  spaced,  the  rib-torsion  stringer-bending  mode  wili  probably  be  the  most  responsive 
of  the  substructure  modes.  In  any  event,  the  computer  program  included  in  Appendix  IV 
should  yield  the  substructure  modes  for  any  particular  design.  It  should  be  noted  that  any 
element  may  be  deleted  from  the  model  by  setting  its  modulus  of  elasticity  to  zero  in  the  in¬ 
put  data . 


E.  Honeycomb  Sandwich  Construction 

For  design  purposes,  a  somewhat  simplified  method  is  proposed  to  relate  the  stress  response 
to  panel  parameters  and  acoustical  excitation.  In  other  words,  simple  plate  theory  can  be 
used  to  suggest  the  significant  panel  parameters  and  theii  mathematical  arrangement  in  the 
empirical  development  of  the  design  nomographs.  The  applicability  of  simple  panel  theory 
is  demonstrated  by  the  correlation  between  measured  and  computed  frequencies  and  stress 
response.  Also,  tests  conducted  by  Bal lentine  (Reference  12),  Sweers  (Reference  13)  and 
Mead  (Reference  14  )  indicate  that  core  shear  does  not  affect  the  first  mode  response  of  flat 
honeycomb  panels,  provided  that  the  core  has  a  specific  weight  of  2  Ib/ft^  or  greater. 

The  analysis  to  be  presented  in  this  section  is  based  on  a  theoretical  model  of  the  honeycomb 
sandwich  panel  edge,  as  shown  in  Figure  9-  This  theoretical  model  is  compared  with  me 
actual  edge  geometry.  Although  the  edge  taper  is  reduced  to  a  step-discontinuity  in  thick¬ 
ness,  it  is  believed  that  this  effect  will  be  negligible. 

Physically,  the  honeycomb  sandwich  panel  model  is  replaced  by  a  simple  plate  with  equiva¬ 
lent  stiffness.  This  approach  presents  analysis  of  the  fiat  honeycomb  panel  on  the  basis  of 
simple  plate  theory.  To  idealize  the  model,  several  simplifying  assumptions  are  introduced 
os  follows: 

o  The  panels  are  flat  and  rectangular. 

o  The  edge  condition  corresponds  to  clamped  and/or  simply  supported  edges. 

o  The  panels  behave  as  thin  plctes.  (Restoring  force  is  due  to  bending  rather  than 

membrane  action.) 
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c  .  f  -  691  .8  Hz 


FIGURE  8.  MODE  SHAPES  FROM  COMPUTER  OUTPUT  (CONTINUED) 


29 


i.  f  -  169.5  Hz 

FIGURE  8  (CONTIN 


Theoretical  Model 


'ff  <<{*?/*  ft' rmTz&z- 


//S7777 

rtSSSSSSSSS 


I 


JZZZ2ZZZZZZL 


sJMCESSsSSsiSn: 


l 


Neutral  Axis 
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Actual 


fe  IS  THE  EDGE  THICKNESS  ONLY  WHEN  METAL  PAN  IS  USED. 


FIGURE  9  THEORETICAL  AND  ACTUAL  MODELS  OF  TAPERED  EDGE  GEOMETRY 


o  The  panel  stiffness  is  derived  from  the  facing  sheet  only, 

o  The  panel  responds  to  random  excitation  only  in  the  first  mode. 

I  .  Frequency  Analysis 

Estimation  of  the  first  mode  vibrational  frequency  for  clamped  edge  panels  will  be  consid¬ 
ered  first,  with  simply  supported  panels  presented  subsequently.  * 


a-  Clamped  Edges  -  By  changing  the  flexural  rigidity  only,  Equations  (20a)  and  (20b)  are 
modified  to  obtain  an  estimate  of  the  first  mode  frequency  of  flat  clarnped-edge  honeycomb 
sandwich  panels  with  thin  facing  sheets,  i  .e. , 
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(81) 


where 


h 

t 

A 


is  the  core  thickness 

is  the  minimum  facing  sheet  thickness 


is  (1  -  v2) 


and 


h<f2>  =  ?  <±> 


2  Ethi 


2A  P 


(82) 


b.  Simply  Supported  Edges  -  Equations  (24a)  and  (24b)  are  modified  by  replacing  the  flex¬ 
ural  rigidity  of  a  plate  with  that  of  a  honeycomb  sandwich  panel  with  thin  facing  sheets  to 
provide  an  estimate  of  the  first  mode  frequency  for  simply  supported  edges.  From  Equation 
(24b)  the  fundamental  frequency  is 


(83) 


2.  Stress  Response 

Estimations  of  the  dynamic  stresses  at  the  edges  and  in  the  facing  sheets  of  flat  honeycomb 
sandwich  panels  when  exposed  to  random  acoustical  noise  is  presented  below  for  clamped 
and  simply  supported  edges. 


a.  Clamped  Edges  -  Equations  (25)  through  (35),  modified  to  reflect  the  honeycomb  sand¬ 
wich  parameters,  are  used  to  estimate  the  maximum  stresses  at  the  edges  of  the  honeycomb 
sandwich  panels.  Equation  (25)  is  rewritten  as 

6  M 

h(ffx)  =  (84) 
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where  t  is  the  metal  edge  thickness  of  the  flat  honeycomb  panel  (see  Figure  9).  Equation 
(35)  then  becomes 


/  V  _  12  /  b  x2 

h(ax)  -  (-) 


max 


(85) 


for  the  maximum  edge  stress  resulting  from  a  uniformly  distributed  static  pressure  of  unit 
magnitude. 
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The  mean-square  edge  stress  response,  considering  light  damping,  for  a  honeycomb  panel 
exposed  to  random  acoustical  noise  excitation  is,  by  combining  Equations  (5)  and  (85), 


36.  ,  b,4  i 

SR2' 


(86) 


An  estimate  of  the  maximum  facing  sheet  stress  is  obtained  from  Equations  (25)  and  (28) 
evaluated  at  the  center  of  the  panel  (x  =  a/2,  y  =  b/2) .  The  resulting  expression  for  the 
facing  sheet  stress  is 


8  »2P 

thb2 


<-)2  + 
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(87) 


Equation  (34)  is  an  applicable  expression  for  W  .  Therefore,  the  maximum  facing  sheet 
static  stress  resulting  from  a  uniform  static  load  of  unit  magnitude  is 


(88) 


and  the  mean-square  facing  sheet  stress  response  for  a  lightly  damped  panel  exposed  to  ran¬ 
dom  acoustic  excitation  is 


36  f.  a4  ffif)  [(b/q)2  +  y 
n3  8  h2t2  R2 


(89) 


b.  Simply  Supported  Edges  -  it  can  be  shown  that  an  estimate  of  the  facing  sheet  stress 
response  can  be  obtained  from  Equations  (25),  (26),  and  (37),  evaluated  at  x  =  a/2, 
y  =  b/2,  as  follows: 


h<°>f  = 


n2  D  W 
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th  b2 


(b/a)z  +  v 


(90) 


where  W  is  evaluated  using  Equation  (40). 

The  maximum  static  facing  sheet  stress  resulting  from  a  uniform  static  pressure  of  unit  mag¬ 
nitude  is  obtained  by  combining  Equations  (40)  and  (90)  as 
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(91) 


An  estimate  of  the  mean-square  stress  response  in  the  facing  sheet  for  a  lightly  damped,  flat 
honeycomb  sandwich  panel  exposed  to  random  acoustic  pressure  excitation  can  be  obtained 
from  Equations  (5)  and  (91).  The  result  is 
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1  .  Introduction 

Based  on  both  theoretical  and  experimental  evidence,  it  has  been  shown  that  curvature 
greatly  influences  the  dynamic  response  characteristics  and  fatigue  resistance  of  aircraft 
structure.  Hence,  cuivature  effects  on  dynamic  response  and  fatigue  applications  are  very 
important  in  the  analysis  of  flight  vehicle  structure. 

Considerable  analytical  work  has  been  done  for  this  research  program,  and  where  possible, 
the  object  has  been  to  provide  design  information.  The  effects  of  curvature  on  the  natural 
frequencies  and  stress  response  of  skin-stringer  and  honeycomb  sandwich  cylindrical  panels 
have  been  studied.  The  relations  developed  and  reported  herein,  for  estimating  natural  fre¬ 
quencies  and  stress  response,  must  be  considered  as  the  first  step  in  a  long-range  program  to 
accurately  determine  curvature  effects.  A  combined  theoretical  program  to  extend  the  ideas 
presented  here  and  an  experimental  program  involving  a  number  of  test  specimens  of  various 
configurations  and  radii  of  curvature  is  required  to  obtain  an  accurate  definition  of  the  dy¬ 
namic  response  characteristics  for  both  skin-stringer  and  honeycomb  sandwich  construction. 


a.  Frequency  Analysis  -  Much  work  has  been  done  in  the  past  few  years  in  the  area  of  shell 
dynamics  (Reference  15)  •  Little,  if  any,  of  this  work  is  capable  of  providing  preliminary 
design  information  (i.e.,  without  extensive  mathematical  analysis  and  computing  faculties) 
in  the  response  of  cylindrical  panels  to  acoustic  excitation.  Indeed,  the  calculation  of 
frequencies  is  no  simple  task  for  any  set  of  specified  boundary  conditions.  The  basic 
approach  has  been  to  use  a  Rayleigh-Ritz  technique  to  estimate  the  fundamental  frequency 
of  a  thin  cylindrical  panel  with  clamped  boundaries.  Assumptions  are  introduced  in  order 
to  provide  an  estimate  suitable  for  design  use.  Also,  several  analytical  methods  are  com¬ 
pared  in  order  to  illustrate  the  effect  of  boundary  conditions  on  the  fundamental  frequency. 


The  approach  used  to  estimate  the  fundamental  frequency  of  a  clamped  cylindrical  panel  is 
the  same  as  described  in  Reference  12  •  The  characteristic  equation  is  determined  from  the 
equations 


For  flat  panels  or  shallow  cylindrical  panels,  the  frequency  of  flexural  vibration  is  approxi¬ 
mated  by 


=  0. 


(94) 
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After  substituting  the  expressions  for  Ggg  and  Jgg,  the  frequency  parameter,  X  ,  is  expressed 


as 


=  ^b3  (1  -v2)  w2 

t2  E  2  2 
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(95) 


where 


a 
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P 

E 

v 

r 

A  =  £ 

b 


s  the  panel  arch  length 
s  the  panel  length 
s  the  panel  thickness 
s  the  panel  mass  density 
s  Young's  modulus 
s  Poission's  ratio 
s  the  panel  radius  of  curvature 


L  =  £,  <f> 

t 


The  curved  panel  coordinate  system  is  shown  in  Figure  10.  The  M's  and  N's  are  defined  in 
Table  I  for  both  clamped  and  simply  supported  edges.  Values  for  the  M‘s  and  N's  are  given 
in  Table  II  for  both  boundary  conditions. 


Substitution  of  the  values  from  Table  II  for  the  fundamental  clamped  modes  yields 
A2  =  41. 7  A  +  (25.2+  L2  2)/A  +  41.7/A3 


(96) 


b.  Stress  Response  -  The  stresses  in  the  panel  are  (Reference  15). 
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TABLE  ! 


DEFINITIONS  OF  M.  AND  N. 

i  i 


Clamped 

Simply  Supported 

(m) 

1 

1/2 

M2(rr.) 

a  R  b  (a  8  b  -  2) 

(£b)2/2 

m 

M3(m) 

(  8  o)4 
x  H  m  1 

(/3mb>4/2 

N]  (n) 

1 

1/2 

<N 

2 

0n  Yn  °  ( 0n  0  “  2) 

n  n  n  n 

(  Xna)2/2 

N3(n) 

( y  a)  4 
'  n  ' 

(yna)4/2 

TABLE  II 

NUMERICAL  VALUES  OF  M.  &  N. 


m,  n 

M2(m),  N2(n) 

M3(m), 

N3(m) 

Clamped 

S.S. 

Clamped 

S.S. 

1 

12.303 

4.93 

500 

48.6 

2 

46.050 

19.72 

3,  800 

778 

3 

98.905 

44.4 

14,618 

3930 

4 

171.586 

78.9 

39,  943 

12,430 

5 

263.998 

123.1 

89,  135 

30, 350 
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The  displacement  functions  u,  v,  and  w  have  been  approximated  in  this  analysis  by 


U  Umn  XL(x)  YnW 

V  =  y  T  v  X  (x)  Y'  (y) 

w  =  yy  W  X  (x)  Y  (y) 
mn  m  '  nv/ 

For  clamped  edges  the  assumed  modes  are 

X  (x)  =  cosh  (8  x)  -  cos  (8  x)  -  a 
Yn(y)  =  cosh  (  Vny)  -  cos  (  y  ny)  -  0n 


(98) 


sinh  (/3  x)  -  sin  (  0  x) 
sinh  (  y  ny)  -  Sin  (  y  ny) 


(99) 


For  skin-stringer  construction,  stresses  are  highest  at  the  center  of  the  longer  edge.  This 
stress  is  normal  to  the  edge  direct  stress.  In  the  case  where  the  straight  edge  is  the  long 
side,  the  maximum  stress,  a  ,  is 


d  v  +  _ t_ 

dy  2 


,  2 
d  w 


(100) 


After  substituting  for  v  and  w  from  Equations  (98)  and  (99)  and  after  evaluating  at  the  center 
of  the  long  side,  the  following  expression  for  the  maximum  stress  is  obtained: 


O  =  y  2  X  (b/2)  W 

y  (i-^2)  n  m  1 


mn 
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where  Vmn  is  determined  from  the  frequency  Equation  (93)  in  terms  of  W '  ,  as  given  below: 


G21G13  ‘  G23  (G11  "  X  mn  Jll) 
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The  term  2K/(  y^t)  can  be  defined  in  terms  of  the  dimensionless  ratios  A,  L,  and  fp  as 
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The  root  mean  square  response  tc  an  arbitrary  acoustic  input  at  resonance  is  taken  from 
Reference  1  2  and  is 
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where 

-2 

Qmn  Is  the  mean  square  generalized  acoustic  force 
M  is  the  generalized  mass 
5mn  is  the  viscous  damping  ratio 
<i>  is  the  natural  frequency. 


(104) 


The  ratio  of  curved  panel  to  flat  panel  root  mean-square  stress  is 
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The  parameters  which  are  affected  by  curvature  in  Equation  (105/  are  and  K  (K  --0  for 
a  flat  panel).  Hence,  the  itic.tr  rcilio  reduces  to 


a  W  f  2K  AL 
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From  Equation  (1C4),  fot  roci  mean-square  response,  the  response  ratio  reduces  to 


(106) 


(10?) 


when  it  is  assumed  that  damping  and  input  force  are  identical  for  both  configurations,  and 
that  the  generalized  mass  is  unaffected  by  small  curvature.  The  stress  ratio  now  becomes 


(108) 
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The  frequency  ratio,  using  the  results  of  Equation  (96),  for  the  first  damped  mode  is 


i 

0,024  AL02 
A4  +  0.61  A2  +  1 


009) 


Data  from  tests  described  in  Reference  12  show  that  the  analytically  determined  constant 
0.024  in  Equation  (109)  should  be  0.006  to  match  the  test  frequencies.  This  difference  is 
attributed  to  the  fact  that  the  test  panel  boundaries  were  not  completely  clamped. 


An  approximate  solution  to  Equation  (102)  for  the  first  mode  is 


2KAL 


Yn 


=  C, 


Az  +  0.0336 
A4 +  9.62  A2  +  1 
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(110) 


and  the  final  result  for  the  stress  ratio  is  obtained  by  substituting  Equations  (109)  and  ( 1 1 0) 
into  Equations  (108): 
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c.  Effect  of  Boundary  Conditions 

Several  methods  of  analysis  were  compared  in  order  to  calculate  the  natural  frequencies  of 
a  cylindrical  panel  wi'h  various  boundary  conditions.  Figure  11  presents  a  comparison  of 
three  different  methods  of  solution  for  the  completely  clamped  panel,  illustrates  the  analogy 
of  the  arch  to  the  panel,  and  presents  the  exact  solution  for  the  simply  supported  panel 
(from  Referenced).  The  methods  used  for  the  clamped  edges  were  (1)  a  double-finite 
Fourier  transform  solution,  (2)  a  Galerkin-Fourier  transform  solution,  and  (3)  a  Rayleigh- 
Ritz  solution.  The  panel  parameters  for  this  example  were  selected  to  match  available 
experimental  data  (Referenced).  Figure  11  clearly  illustrates  the  relative  importance  of 
clamping  along  the  length  as  opposed  to  clamping  the  arc.  The  experimental  points  in 
Figure  11  are  for  an  ideally  clamped  panel.  Clearly,  the  analytical  methods  leave  much 
to  be  desired  when  compared  to  the  experiment.  After  much  study,  it  is  felt  that  this  dis¬ 
crepancy  is  due  to  imperfect  restraint  of  the  inplane  motion  of  the  clamped  edges  in  the 
circumferential  direction.  The  importance  of  this  relaxed  restraint  will  now  be  investigated. 

The  approach  used  is  a  one  term  Rayleigh-Ritz  analysis  with  a  free  parameter,  £ ,  included 
with  the  assumed  mode  for  the  inplane  motion,  v,  which  is  selected  so  as  to  make  the 
fundamental  frequency  a  minimum. 

The  expression  for  the  kinetic  energy  and  the  strain  energy  are  given  below  for  the  classi¬ 
cal  cylindrical  panel;  however,  an  additional  strain  energy  term  is  required  to  include  the 
elastic  effects  of  the  edge.  The  kinetic  energy  is 
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A  solution,  in  the  form  presented  below,  is  assumed: 


u(x,y)  =Tyi  U  X  '  (x)  Y  (y) 

'  p  mn  m  n'7/ 
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w(x,y)  X  (x)  Y  (y) 

'  mn  m  '  n7 

where 

Xm(x)  &  Y  (y)  are  given  by  Equation  (99)  and 

y‘(y)  ‘  0  -f)«in  <i*)  +  f  Y  (y)  018) 

a 

^  is  a  factor  proportional  to  the  degree  of  inplane  restraint,  o  1  £  -  1 


After  assuming  a  one  term  solution  for  u,  v,  and  w  in  Equations  (117)  substituting  this  into 
Equation  (115),  evaluating  the  integrals,  and  differentiating  as  indicated  below 
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=  0 


yields  a  system  of  linear  algebraic  equations.  The  integrals  in  Equations  (119)  are  evaluated 
for  the  most  part  in  Reference  16.  The  system  of  algebraic  equations  has  the  form 
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where 
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and  the  integrals  are  evaluated  as 
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The  eigenvalue  problem  (3  20)  is  solved  by  selecting  a>  for  specified  f  so  that  the  determ¬ 
inant  of  the  coefficient  matrix  vanishes.  The  edge  restraint  parameter,  £  ,  is  selected  so 
2 

that  a>  is  a  minimum.  The  dimensionless  inpiane  spring  constant,  (Ka/C),  is  difficult  to 
determine  in  practice,  but  its  affect  on  the  fundamental  frequencies  of  the  previous  example 
(Figure  11)  is  illustrated  in  Figure  12  where 


q2  =  _Pgbj  ( l.,r  y}\  w2 
Et2 


021) 


The  referenced  experimental  data  in  comparison  with  the  various  analytical  studies 
indicates  that  aircraft-type  panel  construction  does  not  realize  the  degree  of  in-plane 
restraint  represented  in  the  classical  clamped  panel  analysis.  An  analytical  solution  that 
more  nearly  determines  the  fundamental  frequency  of  an  aircraft-type  curved  panel  is 
obtained  by  relaxing  the  in-plane  restraint  along  the  length  through  appropriate  changes 
in  the  energy  representation  and  the  assumed  mode  shapes. 

Having  compared  available  experimental  data,  a  fundamental  frequency  estimate  for  a  curved 
panel  clamped  on  one  side  to  a  frame  with  screws  or  rivets  is  obtained  by  assuming  a  ya lue 
for  Kb/C  of  zero.  If  the  panel  is  clamped  between  two  frames  with  screws  or  rivets, 

Kb/C  =  1.0  gives  the  best  frequency  estimate.  A  more  rigidly  clamped  panel,  of  course, 
requires  a  higher  assumed  spring  constant. 


3.  Cylindrical  Honeycomb  Sandwich  Panels 


The  curved  honeycomb  sandwich  panel  design  techniques  used  in  this  study  relied  almost 
entirely  upon  analytical  methods.  The  stress  data  for  the  curved  sandwich  panel  design  pre¬ 
sented  in  Reference  12  were  used  where  possible. 

Plumblee  (Reference  12  )  analysed  the  curved  sandwich  panel  with  clamped  boundaries  using 
a  Rayleigh-Ritz  method.  The  flexural  modes  were  of  primary  concern;  however,  the  in-plane 
and  shearing  motion  was  also  included.  Wallace  (Reference  17  )  used  Donnell's  method 
(Reference  18  )  for  including  shear  deformation  of  elastic  panels  to  derive  the  natural  fre¬ 
quencies  of  simply  supported  sandwich  panels .  Face  sheet  stress  was  calculated  and  conclu¬ 
sions  were  drawn  that  the  maximum  facing  sheet  stress  occurs  at  the  center  of  the  panel .  The 
maximum  shear  stress,  for  simply  supported  edges,  that  occurs  at  the  middle  of  the  largest 
boundary  was  calculated  . 

Experience  shows  that  fatigue  failure  of  curved,  sandwich  panels  must  be  considered  to  occur 
in  the  thin  edge  member  or  in  the  facing  sheets.  Hence,  the  effect  of  edge  member  restraint, 
and  thus  stress,  can  only  be  taken  into  account  if  clamped  or  elastic  boundary  conditions  are 
used.  Also,  for  design  purposes,  only  a  ratio  of  the  change  in  edge  stress  due  to  chcnges  in 
curvature  will  be  investigated.  This  approach  will  yield  trends  which  indicate  the  effect  of 
panel  parameters  on  fatigue  life  due  to  acoustic  excitation. 

The  above  reasoning  allows  use  of  the  vibration  analysis  of  a  cyl indrical ly  curved  honey¬ 
comb  sandwich  panel  with  clamped  edges,  as  reported  in  Reference  12  .  This  analysis, 
though  too  complex  for  design  purposes,  provides  a  basis  for  a  simplified  solution  and  also 
an  accurate  check  for  this  simplified  solufion-  In  essence,  the  analysis  is  a  free  vibration 
analysis  of  a  curved  sandwich  panel .  The  basic  assumptions  are  that  the  materials,  includ¬ 
ing  the  core,  are  linearly  elastic,  homogeneous,  and  orthotropic.  The  radial  (or  normal) 
displacement  does  not  vary  through  the  shell  thickness.  Normals  to  the  middle  surface  be¬ 
fore  deformation  remain  normal  to  the  middle  surface  after  deformation  (i.e.,  the  core  shear 
displacement  is  linear).  The  facing  sheets  are  thin  and  resist  only  membrane  stresses.  The 
core  resists  only  transverse  shear  stresses.  The  original  analysis  (Reference  12  )  included  in¬ 
plane  and  rotatory  inertias;  however,  only  transverse  inertia  will  be  included  in  this 
analysis. 


A  typical  curved  sandwich  panel  is  shown  in  Figure  13.  The  actual  panel  is  modeled  as  a 
panel  of  uniform  thickness  in  order  to  avoid  complications  in  the  analysis  due  to  the  tapered 
edges  used  in  practice.  The  sandwich  pane!  configuration,  coordinate  system,  and  shell 
element  reactions  are  shown  in  Figures  14  and  15.  While  the  actual  frequencies  calculated 
for  sandwich  panels  using  this  method  may  be  high,  the  ratio  of  curved  panel  to  flat  panel 
frequencies  will  be  sufficiently  accurate  to  provide  reliable  design  information. 


a ■  Frequency  Analysis  -  The  Rayleigh-Ritz  method  is  used  for  the  frequency  analysis  of  the 
clampea  cylindrical  honeycomb  panel  described  above.  The  strain  energy  in  each  layer  of 
the  panel  is  expressed  in  terms  of  the  midsurface  displacements  and  rotations.  Mode  shapes 
which  satisfy  the  clamped  boundary  conditions  are  introduced,  and  finally  the  energy  density 
is  integrated  over  the  volume  of  the  body. 


ft*) 

The  strain  energy  density  in  the  r  layer  of  the  sandwich  panel  is 


u  =  1  L  o.  J 

r  o  2  r  1 
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The  convention  adopted  for  the  stress  nomenclature  is  from  Reference  12  . 
ment  relations  used  are 
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In  the  facing  sheets,  z  in  Equations  (123)  is  replaced  by  h|  (it  is  assumed  that  there  is  no 
rotation  in  the  face  sheets). 


In  matrix  form  the  strain-displacement  relations  are  represented  as 
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where  Ui^  represents  the  displacements  (i.e.;  u^  =  uf  U2  ~  v,  Ug  =  w,  u^  -  h  ,  u^  - 
hj  4>  )  and  J  jA.^  is  an  operator  matrix  for  the  r*h  layer . 

In  Reference  12  the  mode  shape  was  assumed  to  be  a  series  of  clamped  -  clamped  beam 
functions,  and  it  is  shown  that  only  one  term  is  necessary  to  accurately  predict  the  funda¬ 
mental  frequency.  The  displacement  functions  for  the  (m,  n)th  mode  are 
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where  Xm(x)  and  Yn  (y)  are  given  by  Eauation  (99).  The  prime  denotes  a  derivative  with 
respect  to  the  variable.  The  assumed  displacements  are  written  in  matrix  form  as 


- 

UK 

n  k 

^  m  n  k 

(126) 


The  strain  energy  density  in  terms  of  the  generalized  coordinates  is  obtained  by  substituting 
Equation  (126)  into  Equation  (124)  and  then  using  Equation  (122)  . 


The  strain  energy  density  for  each  layer  is  integrated  over  the  layer  thickness  and  the  contri¬ 
butions  summed.  Finally,  the  integration  over  the  panel  area  is  performed. 

For  the  core  the  elastic  constants  are 


,C. .  —  G 
!  44  yz 

and  all  others  are  zero. 
The  elastic  constants  are 


1C55  =  Gxz  <’27> 

The  facing  sheets  are  assumed  to  be  isotropic  and  of  the  same  material. 


2C1 1  =  3C11 

2C22  ~  3C2 2  '  H/(!  - 

2C 1 2  =  3C12  = 
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0  28) 


where  it  has  been  assumed  that  tj  =  tg. 

The  stiffness  matrix  resulting  from  using  the  strain  energy  expression  and  Castigliano's 
theorem  is  presented  in  terms  of  the  following  dimensionless  parameters 

g  =  t/h 

s  =  0  -  ^2)Gyz/E  (129) 

c  =  G  /G 

x/  yz 

The  number  of  variables  that  must  be  considered,  when  the  parameters  (129)  are  used,  is 
reduced  by  two. 

The  kinetic  energy,  neglecting  the  core  mass  is 

T  =  pt2ab  w2  (130) 
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A  =  a/b 


The  eigenvalue  problem  is  then  cast  in  the  form 
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For  completeness,  the  elements  in  the  stiffness  matrix  are  listed  below 
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The  nondimensional  frequency,  fl  ,  is  related  to  the  natural  frequency  by  the  relation 

fl2  =  -P <u2  (132) 

E 

The  nondimensional  terms  A/^m)  and  ^(njare  given  in  Tables  I  and  II,  and  for  the  funda¬ 
mental  mode 

(0,b)  =  (y,a)  =  4.73 


b.  Stress  Response  -  It  is  possible  to  solve  the  eigenvalue  problem  and  obtain  stress  in 
terms  of  the  generalized  coordinates.  Then  by  using  the  generalized  damping  force,  the 
generalized  acoustic  forcing  function,  and  the  generalized  mass  in  LaGrange's  equations  of 
motion,  the  actual  stresses  can  be  obtained.  Experience,  however,  implies  that  this 
method  is  not  an  extremely  accurate  indicator  of  stress  response.  But  in  cases  where  no 
better  method  is  available  these  equations  must  be  used.  Since  there  is  only  a  limited 
amount  of  experimental  data  available  for  acoustically  induced  core  shear,  then  the  core 
shear  stress  must  be  calculated  using  this  method.  A  solution  for  the  fundamental  frequency 
is  first  necessary  to  obtain  the  stresses. 

It  is  recognized  that,  for  a  flat  panel  (0  =  0),  the  eigenvalue  matrix  uncouples  into  a  2  x  2 
matrix  and  a  3  x  3  matrix  from  examination  of  Equation(131) .  The  3x3  matrix  contains  the 
coupling  terms  between  w,  0  and  0  while  the  2x2  matrix  contains  the  coupling  between 
0  and  v,  the  in-plane  motions.  For  small  angles,  the  u-w,  v-w,  and  v-0  coupling  is 
small  and  does  not  significantly  affect  the  frequency.  This  has  been  verified  numerically. 
Therefore,  for  small  angles  the  frequency  can  be  determined  from  the  3x3  (w,0,0) 
matrix.  After  eliminating  higher  order  terms,  the  first  mode  frequency  of  a  shallow 
cylindrical  curved  sandwich  panel  with  clamped  edges  is 

=  ?!  +i£w!a!_  h 

A  A ^  L2  S2 

S,  =  1  +A4  +  40.7-&-,  (1  +  cA2) 

1  scL2 

S  =  |  +40.7  —2-—  ( 1  +-£,  )  +  1655. 

2  scL2  A2 

Restrictions  on  the  use  of  this  equation  are 

0  <  .35  radians 
g  >  10 
s  <  .01 
L  >  100 

If  these  limits  are  exceeded,  the  results  will  not  be  reliable. 

In  order  to  determine  stress  or  stress  ratio,  it  is  necessary  to  determine  the  generalized 

coordinates  U  ,  V  ,0  ,  and  'I'  in  terms  of  W  .  Some  simplication  is  necessary  so 
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that  these  coordinates  can  be  expressed  in  fairly  simple  terms.  The  nature  of  the  generalized 
coordinates  was  studied  as  the  pertinent  parameters  were  varied.  Table  III  shows  the  varia¬ 
tion  of  the  generalized  coordinates  and  natural  frequency  with  variation  in  the  subtended 
angle  0  for  a  given  panel  configuration.  For  angles  less  than  0.35  radian,  t//  and  0are 
relatively  unchanged,  indicating  an  independence  of  the  shell  rotation  angles  with  panel 
curvature.  Also,  the  variation  of  u  and  v  is  nearly  linear  with  0,  for  0  <  0.3  radian. 

These  observations  aid  considerably  in  the  evaluation  of  the  generalized  coordinates. 

TABLE  III 

EFFECT  OF  SUBTENDED  ANGLE  ON  THE  GENERALIZED  COORDINATES 


G 

W 

VV11 

hi*n 

hi<tn 

vn 

un 

ft 

0.0 

.994 

-.0244 

-.0119 

0.0 

0.0 

.04430 

0.05 

.994 

-.0244 

-.0119 

.00329 

.00051 

.04464 

0.1 

.9945 

-.02445 

-.01197 

.006605 

.001017 

.04566 

0.2744 

.9965 

-.02418 

-.01225 

.01888 

.002693 

.05339 

0.5 
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1 .0 
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A  =  a/b  =  1.3939 
L  =  b/h  =  1031.25 
g  =  t/h  =  11.625 
s  =  (1  -  s' 2 )G  =  .0005007 

-  yz 

E 

c  =  Gxz/Gyz  =  1.98895 


When  evaluating  '{'and  <I>  in  terms  of  W,  all  coupling  with  U  and  V  can  be  neglected  as 
demonstrated  in  the  above  example.  The  expressions  for  'I'  and  <I>  are  found  to  be 
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In  this  evaluation,  v  was  assumed  to  be  0.33  and  is  representative  of  most  metals  and  metal 
alloys  used  in  the  fabrication  of  flight  vehicle  honeycomb  sandwich  structure. 


With  the  evaluation  of  'I'  and  <1*  complete,  the  first  two  expressions  of  Equation  (131)  are 
used  to  determine  U  and  V.  The  following  relationships  were  determined: 
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(135) 
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U=- 0J074QW  K  w 

A(1  +  0. 1 A  ;  J 

The  expressions  for  stress  in  terms  of  displacement  are  determined  by  substituting  Equation 
(124)  into  Equation  (122) .  In  expanded  form  the  face  sheet  tensile  stress  and  core  trans¬ 
verse  shear  stresses  are 
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Substituting  the  mode  shapes  of  Equation  (125)  into  Equation  (136)  gives  the  stresses  in  terms 
of  the  generalized  coordinate  W,  the  panel  physical  parameters,  and  the  assumed  mode 
shape.  These  generalized  stress  response  equations  are 
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where  the  primed  variables  are  derivatives  with  respect  to  the  argument.  The  conclusion  to 
be  drawn  from  this  analysis  is  that  there  are  two  regions  which  are  tensile  stress  critical: 
the  center  of  the  panel  and  the  panel  edges.  For  a  flat  panel  with  some  edge  fixity,  the 
maximum  edge  stress  occurs  at  the  center  of  the  long  boundary.  Generalizations  of  this 
type  are  not  possible  with  curved  panels.  The  stresses  in  Equation  (137)  are  only  for  sand¬ 
wich  panels  with  uniform  thickness.  An  approximation  must  be  made  in  the  determination 
of  stress  at  the  edges  for  designs  of  the  type  shown  in  Figure  9. 
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The  change  in  tensile  stress  at  the  center  of  the  panel,  caused  by  curvature  i 
tion  is  evaluated  as 
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Later  in  this  section,  it  is  shown  that,  for  sinusoidal  excitation, 
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From  Equation  (133)  the  ratio  of  frequencies  is 

=  1+0.002  ~  j2 

Substitution  of  Equations  (134),  (135),  and  (140)  into  Equation  (138)  results  in 
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Similarly  the  ratio  of  axc/o^  is 
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The  ratio  of  a  to  o  for  the  curved  panel  is  given  as 
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Another  item  of  iterest  in  structural  design  is  the  ratio  of  stress  in  the  outer  facing  sheet  to 
the  stress  in  the  inner  facing  sheet.  This  ratio  is 
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In  order  to  obtain  a  first  order  evaluation  of  the  effect  of  curvature  on  edge  stress,  the 
effect  of  curvature  for  a  simple  panel  is  examined.  In  Reference  12,  the  ratio  of  stress 
in  a  curved  panel  to  stress  in  a  flat  panel  was  found  to  be,  assuming  sinusoidal  excitation. 
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If  the  frequency  ratio  of  the  sandwich  panel  is  used,  then  the  ratio  of  stresses  becomes 
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In  examining  the  transverse  shear  stress,  it  is  determined  that  the  maximum  value  of  t 
occurs  at  y  =  a/2,  x  «0.22b  and  x  «0.78b.  These  coordinates  were  determined  by 
calculating  the  x  and  y  which  produce  maxima  in  X'and  Y .  The  maximum  in  tyZ  occurs  at 
x  -  b/2,  yw 0.22b  and  y  «0.78b.  For  clamped  edges,  the  transverse  shear  is  zero  at  the 
center  of  the  panel  and  aiong  its  boundaries.  Actually  the  greatest  transverse  shear  occurs 
somewhere  between  the  edge  and  22  percent  of  the  span  as  evidenced  by  the  core  fatigue 
failure  experienced  by  HC-15A.  (See  Table  X,  Section  III.)  It  should  be  pointed  out  that 
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because  of  the  type  loads  resisted  by  the  honeycomb  core,  panel  curvature  does  not  affect 
transverse  shear  stress.  Only  the  natural  frequency  is  changed.  Therefore,  for  sinusoidal 
loading 


r*zc  =  Vc  =  (_£LL> 2  (147) 
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It  is  necessary  to  determine  the  generalized  response  to  an  acoustic  forcing  function  before 
estimates  of  stress  magnitude  can  be  determined.  Lagrange's  equation  which  describes  the 
motion  of  the  generalized  coordinate,  W  ,  is 
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For  sinusoidal  excitation,  the  response  is 
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If  the  forcing  function  is  random  acoustical  noise  (or  if  the  bandwidth  of  the  forcing  function 
is  at  least  three  times  the  bandwidth  of  the  system  at  resonance)  then  the  response  is 
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For  the  case  of  sinusoidal  response,  the  ratio  of  W  to  W  r  becomes 
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The  generalized  force  and  modal  mass  do  not  change  with  frequency  of  the  mode  and  it  is  a 
valid  assumption  that  the  modal  damping  ratio  remains  essentially  constant  for  small  changes 
in  frequency.  In  the  same  manner,  the  generalized  displacement  ratio  for  random  excitation 
is 


(152) 


The  generalized  force  for  the  case  of  sinusoidal  plane  wave  excitation  at  normal  incidence 
is 
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b  a 

Qmn  =  f  f  P(x,  y,  t)  Xm(x)  Yn(y)  dxdy 
o  o 


while  the  generalized  force  for  random  acoustical  noise  excitation  is 


r  r  -  2 

•  i  n  *- 


omn 


J  J  Po  (x,  y,  t)  Xm(x)  Yn(y)  dxdy 


o  o 


—2  2 
where  is  in  units  of  (psi)  /Hz. 

For  uniform  pressure  distribution 
Qn  =  0.691  a  b  P 

Then 

Q  =  0.691  a  b  P 
oil  o 


(154) 


(155) 


The  modal  mass  for  the  sandwich  panel  with  clamped  edges,  assuming  contribution  from  the 
face  sheets  only,  is 


M  r  2  P  t,  a  b 
mn  / 


(156) 


Substitution  of  Equations  (156)  and  (155)  into  Equation  (150)  gives  an  expression  for  rms 
modal  response  to  a  random  normally  incident  pressure  wave.  The  result  is 
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The  stress  response  can  be  estimated  by  combining  Equations  (1 57),  (133),  (134),  (135),  and 
(137). 

The  equations  are  useful  only  for  the  full  depth  honeycomb  sandwich  portion  of  the  panel 
(i.e.,  the  equations  are  not  to  be  used  for  the  thin  tapered  edge). 
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Ill  -  EXPERIMENTAL 


A.  Introduction 

The  objectives  of  the  experimental  investigation  were: 

o  To  increase  the  accuracy  of  existing  skin-stringer  and  honeycomb  sandwich  design 
nomographs. 

o  To  extend  the  range  of  application  of  existing  skin-stringer  and  honeycomb  sandwich 
design  nomographs. 

To  accomplish  the  intent  of  the  investigation,  high-intensity,  random-noise,  sonic  fatigue 
tests  of  120  specimens  were  conducted.  Tests  were  made  on  30  skin-stringer  designs  and  30 
honeycomb  sandwich  designs,  2  each  to  check  repeatability.  The  test  specimen  design 
covered  a  wide  range  of  parameters .  The  skin-stringer  design  parameter  limits  were: 

skin  thickness  0.020  to  0.  100  inch 

rib  thickness  0.027  to  0.  125  inch 

rib  spacing  3.00  to  10.00  inches 

aspect  ratio  1.50  to  3.00 

Honeycomb  sandwich  parameter  limits  were: 

facing  sheet  thickness 
doubler  thickness 
core  thickness 
edge  thickness 
overall  size 
aspect  ratio 

A  description  of  the  test  specimens,  test  procedure,  skin-stringer  tests,  and  honeycomb  sand¬ 
wich  tests  is  given  in  the  following  subsections. 

B.  Test  Specimen  Design 

Detailed  drawings  of  the  test  specimens  are  presented  in  Appendix  Vi;  however,  a  brief  de¬ 
scription  of  them  is  given  here.  A  few  photographs  and  line  drawings  are  also  presented. 

The  test  specimens  were  designed  using  the  applicable  design  nomographs  in  Reference  1  . 
Fastener  edge  distances,  etc.,  were  selected  in  accordance  with  accepted  aircraft  standards. 
Standard  manufacturing  techniques  and  processes  were  used  in  test  specimen  fabrication.  All 
ribs  and  stringers  were  aluminum  alloy  extrusions.  In  a  few  instances,  chemical  milling  was 
required  to  obtain  the  desired  flange  and  web  thickness.  Bonding  and  finishing  processes 
were  checked  against  process  standards.  Qualify  control  procedures  were  used  at  all  times. 
All  test  specimens  were  made  of  7075-T6  clad  aluminum  alloy. 

1  ■  Skin-Stringer 

The  skin-stringer  test  specimens  were  designed  as  shown  in  Figure  16,  with  the  details  shown 
in  Table  IV.  The  number  of  ribs  and  stringers  on  each  test  specimen  was  determined  by  the 
stringer  spacing  and  aspect  ratio  of  each  bay.  For  example.  Figure  1  7  is  a  photograph  of  a 
typical  nine-bay  test  specimen  showing  the  skin  side  and  the  stringer  and  rib  side. 


0-008  to  0.040  inch 

0.015  to  0.090  inch 

0.25  to  0.82  inch 

0-025  to  0.  1 15  inch 

21  x  21  inches  to  37  x  61  inches 

1.0  to  1.7 
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FIGURE  16.  SKIN-STRINGER  TEST  SPECIMEN  DRAWING 


FIGURE  17.  TYPICAL  SKIN-STRINGER  TEST  SPECIMEN 
SHOWING  SKIN  AND  BACK-UP  STRUCTURE  SIDES 
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Each  test  specimen  was  designed  in  accordance  with  Figure  31,  Reference  1 .  MS20426AD 
rivets  were  sized  and  spaced  as  specified  by  the  design  nomograph  for  fasteners.  Figure  35, 
Reference  1 . 


2.  Honeycomb  Sandwich 

The  honeycomb  sandwich  test  specimens  were  designed  using  Figure  41,  Reference  1 .  Fig¬ 
ure  18  is  a  drawing  showing  test  specimen  details.  Figure  19a  is  a  photograph  showing  the 
general  appearance  of  a  typical  honeycomb  sandwich  test  specimen.  An  enlarged  photograph 
of  crushed  honeycomb  core  as  it  appeared  at  the  edge  closure  of  test  specimen  HC-29  isshown 
in  Figure  19b.  Table  V  is  a  listing  of  the  design  details  for  all  thirty  designs.  There  were 
four  test  specimen  sizes,  21"  x  21",  21"  x  25",  25"  x  37",  and  37"  x  61",  and  four  edge 
configurations,  fiberglass  closure  pan,  aluminum  alloy  closure  pan,  crushed  core  edges,  and 
tapered  doubler. 


Doubler  widths  were  determined  at  first  using  a  conservative  rule  of  thumb  based  on  static 
test  results.  After  testing  four  honeycomb  sandwich  beams  described  in  Section  Ill.D,  a 
rule  of  thumb  for  designing  to  meet  sonic  fatigue  requirements  was  established  at  about  the 
half-way  point  in  the  test  program.  Specimen  configuration  HC-15  and  all  subsequent 
doubler  widths  were  established  from  the  new  rule  of  thumb. 

Honeycomb  core  densities  were  determined  using  the  core  shear  design  nomograph  in 
Section  V.C . 


C.  General  Test  Procedure 

The  various  tests  were  conducted  using  the  following  procedure: 

1 .  Modal  Frequency  Studies 

Mode  shapes  for  each  test  specimen  design  were  determined  at  the  beginning  of  the  program. 
The  test  specimen  to  be  observed  was  fastened  to  a  picture  frame  mounting  fixture.  This 
mounting  fixture  was  then  placed  over  two  electromechanical  speakers  arranged  as  shown  in 
Figure  20.  Cork  particles  were  sprinkled  on  the  surface  of  the  vibrating  test  specimen  to 
produce  a  Chladni  pattern  for  each  mode  of  concern.  See  Figure  21.  When  possible,  mode 
shapes  for  each  test  specimen  design  were  determined  through  the (3, 3)  mode. 

2.  Test  Specimen  Instrumentation 

After  a  thorough  study  of  test  specimen  response  characteristics,  strain  gage  locations  were 
selected.  Figure  22  shows  typical  uniaxial  gages  mounted  on  skin-stringer  rest  specimens. 
Typical  uniaxial  gages  mounted  on  honeycomb  sandwich  test  specimens  are  presented  in 
Figure  23. 

A  strain  gage  rosette  was  mounted  at  the  center  of  one  of  the  honeycomb  sandwich  test  spe- 
mens  to  estimate  the  direction  and  value  for  the  maximum  principal  stress.  See  Figure  24 
for  strain  rosette  location. 
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DETAIL  A 


FIGURE  18  .  ALUMINUM  HONEYCOMB  SANDWICH  TEST  SPECIMEN  DETAILS 
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DOUBLER  (TYP  2  PL) 


(a)  Typical  Test  Specimen 


(b)  Crushed  Honeycomb  Core  -  At  Edge  Closure  of  HC-29 


FIGURE  19.  HONEYCOMB  SANDWICH  TEST  SPECIMEN 
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SPEAKERS  INSIDF  BOX 


FIGURE  20.  TEST  SPECIMEN  MOUNTED  FOR  MODAL 
FREQUENCY  STUDY 


FIGURE  21.  CHLADNI  PATTERN  DEVELOPED  FOR 
(1,3)  MODE,  STR-6 
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RIVET  CENTER  LINE 

- 0T - 

CENTER  BAY 
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PANEL  A 


RIVET  CENTER  LINE 

Or - 

CENTER  BAY 

□  » 

I 

I 


PANEL  B 


FIGURE  22.  STRAIN  GAGE  LOCATIONS  FOR  SKIN-STRINGER  PANELS 
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FIGURE  24.  HONEYCOMB  SANDWICH  PANELS 

SHOWING  STRAIN  GAGE  ROSETTES 
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TABLE  V 


3.  Damping  Studies 

An  avarage  damping  ratio,  <5  ,  was  determined  for  first  mode  response,  Jn  every  case  the 
test  specimen  used  was  equipped  with  four  strain  gages.  The  first  mode  was  excited  using 
the  electromechanical  speaker  arrangement  described  previously.  Decaying  strain  oscilla¬ 
tions  were  used  to  determine  S  by  the  log-decrement  method. 

4.  Sinusoidal  Frequency  Sweeps 

Test  specimen  A  and  test  specimen  B  of  a  given  design,  strain  gaged  as  described  above, 
were  mounted  to  a  picture  frame  mounting  fixture.  The  complete  assembly  was  then  mounted 
on  the  progressive  wave  test  section  of  the  High  Intensity  Structural  Test  Facility.  Figures 
25  and  26  pre:.3nt  a  typical  test  arrangement. 

Grazing  incidence  sinusoidal  freauency  sweeps  in  the  range  of  50  -  2000  Hz  were  made  at 
a  sound  pressure  level  of  140db  to  determine  the  test  specimen  strain  response.  Strains 
from  all  six  gages  were  monitored.  They  were  then  studied  and  used  to  shape  the  broad¬ 
band  acoustical  excitation  for  the  fatigue  tests. 

5 .  Broad-Band  Acoustical  Noise  Test  Spectra 

The  test  specimen  picture-frame  assembly  was  removed  from  the  progressive  wave  test  section 
and  replaced  by  aone-inch  thick  plywood  panel.  The  test  sound  pressure  spectrum  was  then 
shaped  out  of  the  presence  of  the  vibrating  test  specimens  and  in  accordance  with  the  strain 
responses  obtained  from  the  sinusoidal  sweeps.  Spectrum  shaping  was  necessary  to  concen¬ 
trate  the  acoustical  energy  in  a  desired  range  of  frequencies. 

The  plywood  panel  was  removed  when  the  desired  spectrum  level  and  shape  were  obtained. 

6.  Fatigue  Tests 

The  test  specimen  and  its  picture  frame  mounting  assembly  were  reinstalled  on  the 
progressive  wave  test  section  ready  for  fatigue  testing  to  take  place.  Inspections  were 
made  following  test  segments  ranging  from  five  minutes  at  the  beginning  to  one  hour  after 
ten  fiours  of  totol  exposure  to  acoustical  excitation. 

A  detailed  description  of  the  test  procedures,  test  facilities,  date  collection  and  reduction 
systems,  and  modal  frequency  study  test  setup  is  given  in  Appendix  I. 

D.  Test  Results 

Results  of  the  skin-stringer  and  honeycomb  sandwich  panel  tests  are  described  in  the  follow¬ 
ing  subsections, 

I .  Skin-Stringer  Test  Results 

a.  Modal  Frequency  Studies  -  Mode  shapes  for  each  skin-stringer  design  were  determined, 
when  possible,  through  the  3,3  mode.  A  summary  of  these  studies  is  presented  in  Table  VI. 

b.  Damping  Ratios  -  In  order  to  obtain  an  average  damping,  seveial  decaying  strain  oscil¬ 
lations  of  the  type  shown  in  Figures  27a  and  27b  were  recorded  for  each  design.  An  arith¬ 
metic  average  of  the  damping  ratio  was  calculated  for  each  design.  Table  VII  is  a  listing 
of  these  average  damping  latios  determined  by  the  log-decrement  method. 
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[FIGURE  25.  HIGH  INTENSITY  TEST  FACILITY  WITH  TEST  SPECIMENS  MOUNTED  FOR  TESTING 


TABLE  VI 

SUMMARY  OF  SKIN-STRINGER  MODAL  FREQUENCY  INVESTIGATION 


m 

I 


n  — 


DESIGNATION 

m,n 

U 

1,2 

1,3 

Modal  Frequency 
2,1  2,2  2,3 

-  Hz 
3,1 

3,2 

3,3 

STR-1 ,STR-2,STR-3, STR-21 

142 

249 

358 

417 

526 

637 

777 

941 

1071 

STR-4 

375 

521 

743 

1006 

1196 

1445 

1638 

2056 

2349 

STR-5,STR-22 

328 

481 

716 

892 

1121 

- 

- 

- 

- 

STR-6 

133 

194 

300 

346 

455 

633 

528 

876 

1032 

STR-7 

156 

249 

313 

357 

481 

564 

709 

830 

918 

STR-8 

468 

776 

1092 

1331 

1683 

2050 

- 

- 

- 

STR-9,STR-24 

162 

192 

384 

268 

326 

- 

500 

535 

876 

STR-1 0,STR-23 

198 

320 

470 

428 

513 

679 

711 

905 

1068 

STR-1 1 

113 

169 

238 

264 

320 

381 

527 

563 

640 

STR-1 2 

220 

331 

417 

531 

641 

774 

846 

1184 

1267 

STR-13 

245 

348 

- 

642 

798 

- 

1054 

1282 

- 

STR-1 4 

139 

192 

277 

310 

344 

501 

628 

690 

812 

STR-1 5 

540 

750 

1160 

930 

1700 

2211 

- 

- 

- 

STR-1 6 

230 

361 

507 

613 

715 

910 

1056 

1320 

1670 

STR-1 7 

183 

287 

396 

451 

611 

739 

811 

- 

1235 

STR-18 

290 

435 

615 

705 

863 

1122 

1311 

1571 

1910 

STR-1 9 

310 

461 

633 

715 

931 

1242 

1337 

1682 

2100 

STR-20 

154 

236 

322 

439 

509 

633 

- 

- 

~ 

STR-25 

335 

571 

706 

815 

992 

1390 

1903 

1811 

2250 

STR-26 

360 

500 

745 

813 

1071 

1405 

1480 

- 

2443 

STR-27 

260 

385 

546 

637 

781 

1025 

1150 

- 

1778 

STR-28 

285 

437 

604 

966 

855 

- 

1296 

1647 

1950 

STR-29 

216 

320 

445 

507 

643 

852 

960 

1131 

1505 

143 

216 

302 

351 

428 

656 

763 

- 

- 

(a)  Decaying  Strain  Oscil lotions,  STR-l 


(b)  Decaying  Strain  Oscillations,  STR-14 


FIGURE  27.  LOG-DECREMENT  DAMPING-SKIN-STRINGER 
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TABLE  VII 


SUMMARY  OF  SKIN-STRINGER  DAMPING  RATIOS 


Test  Specimen  Configuration  Average  Damping  Ratio, % 


STR-1 

1.5 

STR-2 

1.8 

STR-3 

1.4 

STR-4 

1.0 

STR-5 

1.3 

STR-6 

1.0 

STR-7 

1.6 

STR-8 

1.7 

STR-9 

1.7 

STR-10 

1.8 

STR-1 1 

1.1 

STR-1 2 

1.3 

STR-13 

1.7 

STR-1 4 

1.8 

STR-1 5 

1 .5 

STR-1 6 

1.8 

STR-1 7 

1.7 

STR-1 8 

1.3 

STR-1 9 

1 .4 

STR-20 

1.3 

STR-21 

1.5 

STR-22 

1.3 

STR-23 

1.2 

STR-24 

1.5 

STR-25 

1 .3 

STR-26 

1.7 

STR-27 

1 .5 

STR-28 

1.2 

STR-29 

1.4 

STR-30 

1.3 
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c.  Frequency  Sweeps  -  Figures  28  through  30  are  plots  of  strain  response  to  sinusoidal  ex¬ 
citation  at  1 4o  db.  These  frequency  sweep  plots  are  representative  of  the  strain  responses 
for  all  the  rib-stringer  designs. 

d.  Broad-Band  Acoustical  Test  Spectra  -  The  acoustical  noise  test  spectra  were  shaped 
su'jTTniat  the  bandwidth  of  the  excitation  was  at  least  three  times  the  bandwidth  of  the  test 
specimen  response.  Figures  31  through  33  are  representative  test  spectra  for  single  mode 
response.  Figure  32  is  typical  fora  multimode  stress  response. 

e.  Fatigue  Tests  -  Sixty  skin-stringer  test  specimens,  2  each  of  30  designs,  were  fatigue 
tested  until  one  or  mote  cracks  formed  in  the  center  bay  of  the  test  specimen.  Table  VIII  is 
a  summary  of  the  test  results.  Figure  34  shows  the  relationship  between  overall  nominal 

rms  stress  along  the  panel  fastener  row  and  the  number  of  positive  crossings  (cycles  roughly) 
to  failure.  Ninety-five  percent  confidence  limits  were  computed  as  described  in  Appendix 
II  and  as  shown  in  Figure  34. 

Magnetic  tape  recordings  were  made  of  the  signals  from  the  strain  gages  positioned  along 
the  fastener  row.  These  strain  data  were  analyzed  using  a  narrow  band  filter  (nominal  2Hz.) 
to  determine  the  characteristics  of  the  strain  response  to  the  acoustical  excitation.  Figures 
35  through  38  are  the  narrow  band  analyses  of  typical  single  modes  and  multi-modal  responses. 

In  addition  to  the  narrow  band  analysis,  a  10-second  sample  of  the  fastener  row  strain  was 
analyzed  using  a  probability  analyzer  to  determine  the  probability  density  of  instantaneous 
strain  and  probability  distribution  of  strain  peaks.  Figures  39  through  42  are  typical  results 
of  rhe  statistical  analyses. 

2.  Honeycomb  Sandwich  Fatigue  Tests 

a.  Mode  Shapes  -  Results  of  the  modal  frequency  studies  are  listed  in  Table  IX.  It  can  be 
seen  that  very  few  of  the  modes  above  the  (2, 1)  mode  could  be  excited.  Honeycomb  panel 
rigidity  coupled  with  the  insufficient  driving  force  from  the  electromechanical  speaker  sys¬ 
tem  is  believed  to  be  reasons  why  the  higher  modes  were  not  detected. 

b.  Damping  Ratios  -  Damping  ratios  for  the  first  mode  were  determined  as  described  for  the 
skin-stringer  test  specimens.  Typical  decaying  strain  oscillations  are  shown  on  Figure  43. 
Arithmetic  average  damping  ratios  are  listed  in  the  last  column  of  Table  IX. 

c.  Frequency  Sweeps  -  Each  test  specimen  was  subjected  to  frequency  sweeps  made  at  a 
discrete  frequency  sound  pressure  level  of  145  db.  These  sweeps  were  accomplished  to  de¬ 
termine  test  specimen  strain  response  characteristics.  Typical  strain  response  curves  for  the 
honeycomb  sandwich  test  specimens  are  presented.  Figures  44  through  46,  to  show  that  a 
majority  of  the  strain  was  concentrated  in  the  first  mode. 


d,  Broad-Band  Acoustical  Noise  Test  Spectra  -  Broad-band  acoustical  noise  test  spectra 
were  shaped  to  fit  the  discrete  frequency  sweep  strain  response.  Typical  spectra  are  shown 
in  Figures  47  through  49. 

e.  Fatigue  Tests  -  Sixty  flat  honeycomb  sandwich  test  sepcimens  (designs  are  described  in 
Table  V)  were  exposed  to  high-intensity,  broad-band  noise  until  a  fatigue  crack  was  formed. 
Table  X  is  a  summary  of  the  results  of  the  test.  Figure  50  shows  the  overall  rms  nonfinal 
stress  along  the  fastener  row  at  the  edge  of  the  panel  and  the  number  of  crossings  to  failure. 
Figure  51  is  a  fatigue  curve  for  facing  sheet  cracking  determined  in  a  manner  similar  to 
that  for  the  panel  edges.  Ninety-five-and  fifty-percent  confidence  limits  were  calculated 
for  both  fatigue  curves  to  show  tne  residual  variance  about  the  regression  line. 


82 


FIGURE  28.  SINE-SWEEP  EXCITATION  -  $TR  4B 
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FIGURE  29.  SINE-SWEEP  EXCITATION  -  STR  17A 


FIGURE  30.  SINE-SWEEP  EXCITATION  -  STR  26A 


FIGURE  31.  BROAD-BAND.  TEST  SPECTRUM  -  STR  4B 


86 


FIGURE  32.  BROAD-BAND  TEST  SPECTRUM  -  STR  17A 


37 


SPECTRUM  LEVEL,  DB  (RE  0.0002  MICROBAR) 


-  4+4 I4-444+-4444-H-H  H  -  jff  44,, 

^FREQUENCY  -  Hz' 


FIGURE  33.  BROAD-BAND  TEST  SPECTRUM  -  STR  26A 
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TABLE  VIII 

SUMMARY  OF  SKIN-STRINGER  FATIGUE  TESTS 


Desig¬ 

nation 

1 

■  rest 
spectrum 

level,  db 

Response 
frequency 
Hz.  ' 

STR-1A 

123 

200 

B 

j  22 

185 

5TR-2A 

1.23 

12.5 

B 

1 22 

125 

STR-3A 

119 

245 

B 

117 

230 

STR-4A 

127 

350 

B 

126 

460 

STR-5A 

131 

305 

8 

129 

275 

STR-6A 

130 

185 

B 

128 

180 

STR-7A 

118 

225 

B 

118 

210 

STR-8A 

140 

490 

B 

138 

460 

STR-9A 

120 

150 

8 

120 

150 

STR-10A 

123 

22C 

B 

121 

280 

STR-11A 

119 

195 

B 

118 

200 

STR-12A 

123 

2.50 

B 

1 21 

270 

STR-13A 

123 

260 

8 

122 

300 

STR-14A 

128 

130 

B 

126 

150 

STR- 15A 

1  23 

575,750 

B 

122 

575,750 

STR-16A 

133 

270 

B 

131 

250 

STR-17A 

134 

220,290 

B 

133 

240,310 

STR-I8A 

137 

300 

B 

135 

340 

STR-19A 

138 

340 

B 

137 

360 

STR-20A 

137 

165 

B 

135 

175 

STR-21A 

125 

220 

B 

125 

220 

STR-22A 

131 

270 

Number 
of  crossings 
to  fcilure; 
Nc  x  10*6 


17.50 
9.10 

13.50 
11.30 
27.00 
16.00 
12.00 
27.00 

1.70 

10.00 

3.65 

3.00 

21.50 

20.20 

9.70 
11.60 

1.61 
1.51 
1 . 20 
1.20 
21.60 
22.00 

21.50 
23.80 
28.00 
33.00 

12.40 

14.90 
210.00 
208.00 

10.00 

9.70 

5.50 
6. CO 

1.90 
3.10 

4.40 

4.50 
4.50 

4.70 
21 .70 
23.00 

nn 


Overall 

Nominal 

nns 


strew,  ksi 


Remarks 


5.22 

4.91 

4.53 


lipped  z 
type  stringer 


4.00 


5.05 

4.10 
5.23 

5.81 
5.50 
4.71 
6.40 
6.74 
4.00 
3.35 

5.11 

3.81 
6.20 
6.52 


1 , 2  mode 
response 


1 , 2  mode 
response 


4.41 
5.10 

5.10 
4.07 
5.22 

4.80 
4. 15 
5.05 
3.25 
2.90 

6.10 
6.88 
6.21 

5.81 
7.65 
8.45 

7.35 
7. 17 
6.85 

4.35 

3.82  1,2  mode 

3.92  response 

n  ns 


results 

questionable 
1 , 2  mode 
response 


TABLE  VIII  (Cont'd) 


Desig¬ 

nation 

Test 

spectrum 
level,  db 

Response 

frequency, 

Hz. 

Number 
of  crossings 
to  failure, 

Nc  x  10-6 

| 

Overall 

Nominal 

rms 

stress,  ksi 

S7R-23A 

123 

260 

17.00 

3.75 

B 

121 

240 

15.50 

3.41 

STR-24A 

120 

150 

14.50 

3.62 

B 

119 

145 

16.20 

4.61 

STR-25A 

133 

350,600 

18.50 

5.52 

B 

132 

350,600 

21.60 

4.75 

SYR-26A 

137 

345 

14.00 

4.27 

8 

135 

390 

13.00 

6.10 

STR-27A 

137 

240,400 

11.50 

5.45 

B 

136 

240,400 

29.00 

6.88 

STR-28A 

137 

310 

10.00 

6.61 

B 

137 

290 

13.00 

6.61 

STR-29A 

130 

220 

9.20 

6.91 

B 

129 

240 

8.10 

7.10 

STR-30A 

140 

165 

100.00 

3.15 

B 

138 

165 

145.00 

2.80 

Remarks 
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400  600  800  1000  ]200  1400 

FREQUENCY  -  Hz 

FIGURE  35.  NARROW-BAND  STRAIN  ANALYSIS  -  SIR  4B 
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FIGURE  36-  PSD  STR, 


FIGURE  37.  NARROW-BAND  STRAIN  ANALYSIS  -  STR  21 B 


FREQUENCY  -  Hz 


FIGURE  38.  PSD  STRAIN  ANALYSIS  -  STR  21B 
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FIGURE  40.  AMPLITUDE  DISTRIBUTION  OF  STRAIN  -  STR-4B 
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RMS  RATIO 


TABLE  IX 


SUMMARY  OF  HONEYCOMB  SANDWICH  MODAL  FREQUENCY 
AND  DAMPING  INVESTIGATION 


Designation 

Modal  frequency*  Hz 
1,2  1,3 

2,1 

Average  damping 
for  (1 , 1)  mode,  % 

HC-1 

229 

451 

925 

632 

1  .5 

150 

HC-2,  HC-3 

250 

521 

- 

- 

1 .8-2.1 

HC-4 

230 

439 

895 

677 

1 .3 

HC-5,  HC-22 

205 

416 

863 

550 

2.4 

145 

HC-6 

190 

335 

780 

536 

1 .5 

HC-7 

255 

475 

- 

- 

1.5 

HC-8 

515 

- 

- 

- 

1.7 

HC-9 

240 

452 

816 

692 

1.5 

HC-10 

275 

359 

900 

- 

1 .6 

215 

HC-1 1 

197 

366 

785 

605 

1.8 

HC-1 2 

282 

635 

- 

- 

2.6 

195 

HC-1 3 

256 

503 

1015 

705 

2.4 

HC-1 4 

248 

513 

950 

715 

1.6 

HC-1 5 

201 

425 

805 

576 

2.0 

HC-1 6,  HC-25 

233 

535 

935 

637 

1.8 

HC-1 7 

206 

417 

885 

582 

1.9 

HC-1 8 

493 

_ 

- 

- 

2.7 

HC-1 9 

2  77 

464 

- 

- 

2.0 

HC-20 

262 

397 

702 

- 

1 .5 

202 

HC-21 

246 

509 

1120 

711 

1 .5 

HC-23 

286 

- 

- 

- 

2.2 

HC-24 

213 

480 

905 

595 

2.2 

HO- 26,  HC-28 

191 

312 

647 

610 

2.2 

HC-27 

213 

413 

722 

635 

1.3 

HC-29 

217 

405 

742 

682 

1.7 

HC-30 

118 

165 

- 

- 

1.9 

95 

*Modes  higher  than  (1,3)  couid  not  be  excited. 
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FIGURE  45.  SINE-SWEEP  EXCITATION  -  HC  20A 
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PEAK  STRAIN  -  MICRO  IN/IN 


FIGURE  46.  SINE-SWEEP  EXCITATION  -  HC  22A 
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SPECTRUM  LEVEL,  DB  (RE  0.0002 


FIGURE  48.  BROAD-BAND  TEST  SPECTRUM  -  HC  20A 
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TABLE  X 


SUMMARY  OF  HONEYCOMB  SANDWICH  FATIGUE  TESTS 

Overal  I 


No.  of  rms  stress 

crossings  at  failure 


Designation 

Test  spectrum 
level,  db 

Response 
frequency,  Hz 

to  failure, 
Nc  x  10  "6 

point 

ksi 

Remarks 

HC-1  A 

138 

230 

0.83 

6.12 

edge  crack 

B 

136 

240 

0.87 

5.98 

edge  crack 

HC-2  A 

138 

250 

1.80 

7.21 

edge  crack 

B 

137 

250 

1.80 

5.45 

edge  crack 

HC-3  A 

134 

270 

1.90 

5.36 

edge  crack 

B 

134 

270 

1.90 

5.11 

edge  crack 

HC-4  A 

134 

240 

18.60 

2.72 

edge  crack 

B 

132 

250 

14.20 

3.21 

edge  crack 

HC-5  A 

133 

210 

0.37 

8.35 

edge  crack 

B 

131 

220 

0.39 

7.25 

edge  crack 

HC-6  A 

134 

200 

0.54 

7.35 

edge  crack 

B 

134 

200 

0.72 

6.81 

edge  crack 

HC-7  A 

140 

270 

5.80 

4.61 

edge  crack 

B 

139 

280 

6.00 

4.12 

edge  crack 

HC-8  A 

138 

470 

14.40 

edge  crack 

B 

136 

500 

13.50 

edge  crack 

HC-9  A 

140 

250 

4.50 

3.82 

edge  crack 

B 

138 

240 

4.30 

4.55 

edge  crack 

HC-10  A 

138 

220 

0.27 

8.23 

edge  crack 

B 

138 

230 

0.28 

8.11 

edge  crack 

HC-1 1  A 

138 

195 

11.50 

5.75 

facing  sheet  crack 

B 

137 

190 

11.80 

4.25 

facing  sheet  crack 

HC-1 2  A 

138 

290 

18.20 

4.50 

facing  sheet  crack 

B 

137 

290 

17.20 

5.35 

facing  sheet  crack 

HC-1 3  A 

140 

250 

11.70 

6.80 

facing  sheet  crack 

B 

139 

265 

14.30 

5.10 

facing  sheet  crack 

HC-1 4  A 

140 

255 

29.50 

5.80 

facing  sheet  crack 

B 

138 

260 

29.50 

5.00 

facing  sheet  crack 

HC-1 5  A 

140 

200 

1.08 

10.79 

core  shear  failure 

B 

140 

200 

1.90 

10.50 

facing  sheet  crack 

HC-1 6  A 

140 

220 

1.60 

8.00 

data  not  valid 

B 

138 

240 

17.00 

6.25 

facing  sheet  crack 

HC-1 7  A 

140 

210 

17.40 

5.00 

facing  sheet  crack 

B 

138 

210 

18.10 

5.15 

facing  sheet  crack 

HC-1 8  A 

138 

260 

11.40 

5.50 

facing  sheet  crack 

B 

137 

285 

13.40 

5.05 

facing  sheet  crack 

HC-1 9  A 

140 

290 

22.50 

4.10 

facing  sheet  crack 

B 

139 

260 

20.50 

4.00 

facing  sheet  crack 

HC-20  A 

138 

210,220 

6.80 

5.50 

facing  sheet  crack 

B 

138 

210,220 

13.60 

5.15 

facing  sheet  crack 

HC-21  A 

131 

220 

16.20 

5.50 

facing  sheet  crack 

B 

129 

220 

16.20 

6.05 

facing  sheet  crack 

HC-22  A 

134 

270 

3.78 

7.61 

edge  crack 

B 

132 

275 

4.80 

7.33 

edge  crack 
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TABLE  X  (Cont'd) 


Designation 

Test  spectrum 
level,  db 

Response 
frequency,  Hz 

No.  of 
crossings 
to  failure 

Nc  x  10 

Overal  1 
rms  stress 
at  failure 
point 
ksi 

Remarks 

HC-23  A 

136 

260 

13.00 

6.20 

facing  sheet  crack 

B 

136 

250 

18.00 

5.95 

facing  sheet  crack 

HC-24  A 

136 

150 

3.45 

8.12 

edge  crack 

B 

334 

145 

3.96 

7.53 

edge  crack 

HC-25  A 

138 

350,600 

6.50 

6.50 

facing  sheet  crack 

B 

137 

350,600 

14.70 

4.80 

facing  sheet  crack 

HC-26  A 

137 

345 

1.38 

8.15 

edge  crack 

B 

135 

390 

1.11 

8.36 

edge  crack 

HC-27  A 

138 

240,400 

6.70  &  8.70  5. 

81  &  4.90  edge  crack  - 

facing  sheet  crack 

B 

137 

240, 400 

5.70  &  7.70  5. 

15  &  5.03 

edge  crack  - 
facing  sheet  crack 

HC-28  A 

138 

310 

1.10 

8.42 

edge  crack 

B 

137 

290 

1.44 

7.65 

edge  crack 

HC-29  A 

137 

220 

10.20 

3.13 

edge  crack 

B 

135 

240 

11.20 

2.76 

edge  crack' 

HC-30  A 

133 

1.30 

6.18 

edge  crack 

B 

130 

1.00 

6.32 

edge  crack 

109 


FIGURE  50  HONEYCOMB  SANDWICH  EDGE  MEMBER  FATIGUE  CUkVE 


FIGURE  51.  HONEYCOMB  SANDWICH  FACING  SHEET  FATIGUE  CURVE 


A  number  of  selected  tape  recorded  strain  gage  signals  were  analyzed  to  determine  the 
nature  of  their  statistical  properties.  Figures  52  through  55  present  the  statistical  analysis 
of  edge  strains  and  facing  sheet  strains.  Probability  density  of  instantaneous  strain  and 
probability'  distribution  of  strain  peaks  are  included. 

Edge  and  facing  sheet  strain  recordings  were  analyzed  with  a  2- Hz  bandwidth  filter  to  de¬ 
termine  the  characteristics  of  the  strain  response  to  the  random  acoustical  excitation.  The 
narrow-band  analyses  were  also  compared  with  the  discrete  frequency  sweeps.  Figures  56 
through  59  are  typical  analyses. 

f.  Edge  Doubler  Tests  -  The  edge  of  a  honeycomb  sandwich  panel  is  generally  comprised  of 
a  Facing  sheet  thickness,  a  doubler,  and  a  closing  member,  sometimes  called  a  pan.  The 
doubler  is  necessary  to  maintain  the  required  strength  at  the  edges  where  the  load  in  the 
smaller  facing  sheet  is  transferred  to  the  larger  facing  sheet.  Obviously,  the  width  of  the 
doubler  is  dependent  upon  the  portion  of  the  panel  span  required  to  transfer  load  across 
the  full  depth  core.  Tne  lengtn  of  span  required  to  transfer  this  load  is  a  function  of  the 
bevel  angleof  the  core.  It  is  very  important  that  a  minimum  width  doubler  be  used  since 
it  contributes  significantly  to  the  total  panel  weight.  It  is  believed  that  rules  of  thumb 
based  on  static  tests  have  been  used  heretofore  to  establish  doubler  width. 

A  set  of  simple  experiments  were  made  in  an  attempt  to  determine  what  width  doubler  should 
be  used  for  a  given  core  bevel  angle.  Four  simple  beams  of  the  configuration  shown  in 
Figure  60  were  fabricated  and  tested.  Strain  gages  were  attached  to  the  facing  sheets  and 
were  spaced  one  half  inch  apart  at  the  ends  and  one  inch  apart  elsewhere  along  the  center- 
line  of  the  beam  to  measure  bending  strain.  Figures  61a  and  61b  show  a  typical  beam  be¬ 
fore  and  after  strain  gage  installation.  The  beams  were  mounted  one  at  a  time  over  the 
electro-mechanical  speaker  system  used  for  modal  frequency  studies,  with  the  speaker  pro¬ 
ducing  discrete  frequency  sound  tuned  to  the  fundamental  frequency  of  each  beam.  Strain 
measurements  were  recorded  at  two  arbitrarily  selected  sound  pressure  levels,  140  and  145 
db. 

Figure  62  is  a  piot  of  the  strain  distribution  of  all  configurations  normalized  to  the 
maximum  strain  magnitude  recorded  at  zero  percent  span  of  Configuration  A  (1  .25-inch 
wide  doubler) . 

The  strain  distributions  shown  in  Figure  62  indicate  that  the  doubler  significantly  influences 
the  spanwise  strain  distributipn .  The  greatest  influence  occurs  at  the  end  of  the  doubler,  a 
very  critical  spot  where  fatigue  cracks  frequently  form.  Also,  it  should  be  noted  that  the 
strain  experiences  a  change  of  sign  (tension  to  compression  or  vice  versa)  near  the  line 
where  the  thin  edge  ends  and  the  beveled  core  begins.  Also  in  this  same  region,  the 
strain  was  reduced  approximately  32  percent  with  the  2.15-inch  doubler  as  compared  to 
the  1 .25-inch  doubler.  The  2.60-inch  doubler  only  reduced  the  strain  36  percent  as  com¬ 
pared  to  the  1  .25-inch  doubler,  however.  In  every  case,  the  strain  experiences  a 
minimum  at  16.3  percent  of  the  span. 

The  curves.  Figure  62,  indicate  that  the  design  with  the  2. 15-inch  doubler  was  almost  as 
effective  as  the  design  with  the  2.60-inch  doubler.  As  a  result  it  was  decided  that  narrow 
doublers  would  be  used  on  test  specimens  HC-15  through  HC-30.  These  test  specimens  were 
subsequently  tested  with  no  observable  reduction  in  fatigue  resistance,  but  with  a  signifi¬ 
cant  reduction  in  weight.  Therefore,  it  is  recommended  that  the  doubler  width  for  design 
be  established  as  shown  in  Figure  63  * 


112 


PROBABILITY  OF  PEAKS  EXCEEDING 


o’  o’  o’  o 


A1ISN3G  Ainavgoad 


114 


FIGURE  53.  AMPLITUDE  DISTRIBUTION  OF  STRAIN  -  HC-12B 
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OF  STRAIN  PEAKS  -  HC-20A 


PROBABILITY  OF  PEAKS  EXCEEDING 


FIGURE  55.  PROBABILITY  DISTRIBUTION  OF  STRAIN  PEAKS  -  HC-22A 
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RMS  STRAIN  -  MICRO  IN/lN/VHz  1 


FIGURE  56.  NARROW-BAND  STRAIN  ANALYSIS  -  HC  I2B 
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RMS  STRAIN  -  MICRO  IN/IN/Viiz 


SPECTRUM  LEVEL:  134  db  I  I  I  !  i  !  |  I 

•j  FILTER  BANDWIDTH:  2  Hz,  Avg  Time  Constant:  2  sec. 
i  SCANNING  SPEED:  0.5  Hz/sec,  Sample  Time:  10  sec 
! OVERALL  STRAIN:  761  u  in/in  1  r  T~~1  X 


MEAN  SQUARE  STRAIN  DENSITY  - 


■Hrrm: 


FIGURE  59.  PSD  STRAIN  ANALYSIS  -  HC  22A 
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DESIGN  DETAILS 


SPANWISE  STRAIN  DISTRIBUTION 


(a)  Before  Strain  Gage  Installation 


(b)  After  Strain  Gage  installation 


FIGURE  61.  TYPICAL  HONEYCOMB  SANDWICH  BEAM 
FOR  DOUBLER  INVESTIGATION 
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FIGURE  62.  SPAN  WISE  STRAIN  FOR  HONEYCOMB  SANDWICH  BEAM 


WHERE  Ld  IS  THE  WIDTH  OF  DOUBLER  NECESSARY . 

L  |S  THE  WIDTH  OF  DOUBLER  FROM  THE  EDGE  OF  THE  PANEL 
dl  TO  BEGINNING  OF  FULL  DEPTH  CORE. 

L  IS  THE  ADDITIONAL  WIDTH  NEEDED  AND  EQUAL  TO 
d2  THE  CORE  THICKNESS. 

FIGURE  63.  SKETCH  OF  RECOMMENDED  DOUBLER  WIDTH  TO  MEET  SONIC 
fatigue  requirements 
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E.  Discussion  of  Test  Results 


1 .  Skin-Stringer 

a .  Mode  Shapes  -  Some  of  the  modes  could  not  be  excited.  From  Table  VI,  it  can  be  seen 
that  the  (2,  3),  (3,  1),  (3,  2)  and  (3,  3)  modes  were  unresponsive  in  a  number  of  instants. 

The  lack  of  response  is  attributed  to  the  inability  of  the  electro-mechanical  speakers  to  sup¬ 
ply  sufficient  excitation  energy  at  the  higher  frequencies. 

A  comparison  of  the  skin-stringer  modal  frequencies  determined  for  the  modal  frequency 
investigation  and  the  frequencies  resulting  from  the  frequency  sweeps  show  that  the  frequen¬ 
cies  for  the  first  mode  do  not  correlate  too  well .  it  is  believed  that  the  differences  can  be 
attributed  to  ambient  temperature  differences.  The  temperature  of  the  test  specimens  was 
approximately  70°F  during  the  modal  frequency  investigation,  which  was  conducted  in  a 
well -insulated  anechoic  room.  Test  specimen  temperature  during  frequency  sweeps  was 
measured  to  be  as  low  as  58°F  and  as  nigh  as  70°F,  depending  on  the  test  chamber  ambient 
temperature  and  the  mass  flow  of  air  through  the  progressive  wave  test  section.  It  is  well 
known  that  modal  frequencies  for  thin  plates  can  vary  greatly  with  small  changes  of  temper¬ 
ature. 

b .  Frequency  Sweeps  -  Strain  activity  above  the  first  mode  was  less  than  originally  expected 
All  but  seven  of  the  skin-stringer  designs  showed  only  one  significant  mode  above  the  funda¬ 
mental,  usually  the  (1,  2)  mode.  Figure  29,  the  sweep  response  for  STR-17A,  is  an  illustra¬ 
tion  of  strain  response  occurring  in  a  number  of  modes.  STR-26A  demonstrated  a  response 
similar  to  that  of  STR-17A. 

As  described  above,  significant  frequencies  found  during  the  modal  frequency  studies  did 
not  accurately  compare  with  the  frequency  sweeps.  This  lack  of  correlation  is  attributed 
to  dissimilar  ambient  temperature. 

In  some  instances,  stringer  and  rib  torsion  and  stringer  and  rib  bending  influenced  the  re¬ 
sponse  of  the  skin.  Figure  30  presents  the  stringer  and  rib  torsion  and  bending  and  rib  in¬ 
fluence  on  skin  response  for  STR-26A. 

c.  Fatigue  Tests  -  Three  types  of  failures  are  commonly  experienced  by  skin-stringer 
construction.  The  most  common  failure  is  the  formation  of  fatigue  cracks.  These  generally 
form  around  the  fasteners  where  stress  raisers  exist.  The  cracks  propagate,  if  the  acoustical 
excitation  is  allowed  to  persist,  until  they  include  a  number  of  fasteners.  Figure  64a  is  a 
photograph  of  this  type  of  fatigue  crack  which  grew  quite  rapidly  with  continued  exposure 
to  the  random  acoustical  excitation. 

The  prying  action  of  the  attached  stringer  or  rib  flange  on  the  fasteners  induces  high  stresses 
in  them,  and  fastener  fatigue  and  popped  heads  result.  This  type  of  failure  was  also  expe¬ 
rienced  during  the  sonic  fatigue  tests.  Figure  64b  is  a  photograph  of  fastener  fatigue  and 
popped  rivet  heads  typical  of  this  type  of  fatigue  failure. 

The  prying  action  of  the  attached  flange  on  the  fasteners  also  cause  stresses  in  the  attached 
flange,  which  may  be  critical  if  the  stiffener  thickness  is  smaller  than  the  skin  thickness, 
or  if  the  flange  is  quite  wide  with  a  large  bend  radius.  This  type  of  fatigue  failure  was  not 
experienced  at  anytime  during  the  test  program,  indicating  that  stringer  and  rib  flange 
thicknesses  determined  from  Figure  31,  Reference  1  ,  result  in  conservative  designs. 

Table  VIII,  the  summary  of  skin-stringer  fatigue  tests,  shows  that  the  test  spectrum  was  not 
the  same  for  specimens  A  and  B  in  every  case.  Average  sound  pressure  level  measurements 
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r-  % 


(b)  Typical  Fatigue  Crack  and  Popped  Rivet  Heads 


FIGURE  64.  SKIN-STRINGER  FATIGUE  FAILURES 
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made  during  the  tests  showed  a  maximum  variation  of  2  decibels.  A  1  decibel  difference 
was  common.  The  difference  of  sound  pressure  at  specimens  A  and  B  was  probably  due  to  the 
normal  modes  of  the  progressive  wave  test  section. 

Four  of  the  skin-stringer  configurations  showed  predominant  strain  response  in  the  (1,  2) 
mode,  instead  of  the  (1,  1)  mode.  The  (1 ,  1)  mode  strain  response  of  STR-3,  STR-7,  STR-11 
and  STR-21  was  considerably  less  than  the  (1,  2)  mode  strain.  A  stroboscope  was  used  in 
each  case  to  determine  which  mode  was  predominant.  The  stroboscope  study  confirmed  the 
belief  that  the  (I,  2)  mode  was  the  strongest  in  each  case.  No  explanation  is  offered  as  to 
why  this  happened,  except  possibly  the  (1,  1)  mode  damping  was  more  than  the  (1,  2)  mode 
damping.  Figure  37,  a  narrow-band  analysis  of  panel  strain  from  the  fastener  line  of 
STR-21B,  is  an  illustration  of  this  type  response. 

The  number  of  positive  crossings  to  failure  was  determined  from  unfiltered  oscillograph  traces. 
Ten-second  samples  of  the  type  shown  on  Figure  65  were  used  to  determine  the  number  of 
positive  crossings  per  second. 

The  fatigue  curve.  Figure  34,  shows  the  characteristic  data  scatter,  especially  above  107 
crossings  to  failure.  Between  the  ±95-percent  confidence  limits,  the  scatter  factor  bounded 
by  the  ± 95-percent  confidence  limits  is  17.85. 

The  fatigue  curve  for  the  skin-stringer  tests  should  be  compared  to  fatigue  curves  for  similar 
material,  loading  and  stress  raisers.  Royal  Aeronautical  Society  fatigue  data.  Reference  19  , 
were  selected  for  the  comparison.  Table  XI  is  a  summary  of  the  comparison. 


TABLE  XI 

COMPARISON  OF  RAS  FATIGUE  DATA  WITH  SKIN-STRINGER  FATIGUE  DATA 


Description 

Inverse  negative  slope 

Regression  line 

Scatter  factor 

of  fatigue  curve 

rms  stress,  ksi 

bounded  by  ±95% 

106 

107  108 

conf.  limits 

RAS 

4.59 

8.95 

5.39  3.24 

23.10 

Skin-Stringer 

Plating 

4.55 

8.80 

5.25  3.18 

17.85 

The  data  tabulated  above  indicate  that  the  skin-stringer  fatigue  data  and  RAS  fatigue  data 
compare  favorably. 


Ninety-five-percent  confidence  limits  were  established  for  the  mean  fatigue  curve  of  the 
skin-stringer  design  nomograph.  Figure  31,  Reference  1  .  This  was  done  by  assuming  that 
the  fatigue  curve  of  Figure  31,  Reference  1,  had  the  same  mean  stress  and  standard  devia¬ 
tion  as  the  fatigue  curve  for  reversed  bending  across  a  rivet  line  contained  in  Reference  20. 
This  assumption  may  not  have  been  extremely  accurate,  but  it  was  the  only  known  way  of 
getting  an  indication  of  the  95-percent  confidence  limits.  The  confidence  limits  were 
needed  to  establish  a  go/no-go  test  for  determining  whether  or  not  additional  sets  of  a 
given  design  should  be  tested  to  establish  structural  reliability. 
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The  rule  established  was  that  if  one  of  the  two  specimens  tested  at  a  given  condition  exhibit 
a  life  which  fell  outside  the  ±95-percent  confidence  limits,  the  status  of  the  test  results 
relative  to  the  design  nomograph  was  no-go.  The  fatigue  life  of  both  replicates  had  to  fall 
within  the  ±95-percent  confidence  limits  for  the  status  to  be  go.  Only  4  of  the  30  designs 
fell  in  the  go  category.  In  fact,  only  15  percent  of  the  test  specimens  tested  had  lives 
which  fell  in  the  ±95-percent  confidence  limits.  Furthermore,  6.65  percent  of  the  60  test 
specimens  tested  fell  below  the  lower  95-percent  confidence  limit  and  78.4  percent  were 
above  the  95-percent  confidence  limit,  these  results  indicate  that  the  design  nomograph 
for  skin-stringer  construction  in  Reference  1  is  conservative  in  the  range  of  test  conditions 
reported  herein. 

STR-1  and  STR-2  were  identical  designs  except  the  stringers  on  STR-2  had  bent  flanges. 

They  were  tested  under  the  same  conditions  to  determine  if  the  added  stringer  flange  stiff¬ 
ness  influenced  the  fatigue  resistance  of  the  design.  Apparently  it  did  not  since  the  fatigue 
life  in  each  case  was  about  the  same  (see  Table  VIII). 

STR-3  design  was  retested  as  STR-21  because  the  first  test  resulted  in  one  specimen  falling 
in  the  go  category  and  the  other  specimen  in  the  no-go  category.  STR-21  tests  decided 
the  issue  and  placed  the  design  into  the  go  category;  both  STR-21A  and  STR-21B  demon¬ 
strated  lives  which  fell  within  the  ±95-percent  confidence  interval  described  above. 

STR-5  was  retested  as  STR-22  because  the  fatigue  life  of  STR-5B  was  approximately  ten 
times  that  of  STR-5A.  The  tests  of  STR-22  showed  conclusively  that  the  time  to  failure 
experienced  by  STR-5A  was  not  representative  of  the  fatigue  resistance  of  the  design. 

The  results  of  the  tests  involving  STR-10  were  questionable  due  to  poor  test  specimen 
design.  This  design  should  have  had  a  greater  fatigue  resistance  than  STR-1,  STR-2,  and 
STR-3.  It  did  not!  The  adjacent  bays  of  the  test  specimen  were  slightly  larger  than  the 
most  important  center  bay.  It  is  believed  the  adjacent  bays  were  excited  in  their  first 
mode  because  the  bandwidth  of  the  sound  pressure  excitation  was  wide  enough  to  develop 
their  full  (1,1)  mode  response.  The  cracks  which  formed  on  the  fastener  lines  probably 
resulted  from  the  vibration  of  the  slightly  larger  adjacent  bays.  STR-23  was  designed  as  a 
retest  of  STR-10.  The  adjacent  bays  of  STR-23  were  the  same  size  as  the  center  bay. 

STR-23  test  results  showed  that  the  STR-10  specimen  design  was  poor.  The  time  to  failure 
was  comparable  to  STR-1,  STR-2,  and  STR-3,  within  acceptable  scatter  of  fatigue  data. 

The  experimentally  determined  probability  distribution  of  fastener  line  instantaneous  strain 
and  strain  peaks.  Figures  39  through  42,  show  that  both  probability  distributions  for  each 
of  the  designs  deviate  somewhat  from  the  Gaussian  and  Rayleigh.  The  deviation  from  the 
theoretical  distributions  is  attributed  to  nonlinear  response  and  multimodal  influences. 

2.  Honeycomb  Sandwich 

a.  Mode  Shapes  -  From  Table  IX  it  can  be  seen  that,  in  general,  modes  above  the  (2, 1) 
were  unresponsive.  In  some  instances  the  (1,3)  mode  could  not  be  detected  by  the  cork 
particles.  The  lack  of  response  is  attributed  to  insufficient  excitation  energy  from  the  elec¬ 
tromechanical  speakers  at  the  higher  frequencies  where  the  generalized  force  is  reduced. 

HC-1 ,  HC-5  and  HC-22,  HC-10,  HC-12,  HC-20,  and  HC-30  showed  two  (1,1)  mode 
Chladni  patterns  at  two  frequencies.  This  is  attributed  to  the  tapered  edge  design,  because 
at  the  lower  frequency  the  mode  lines  formed  very  close  to  the  fastener  row.  At  the  higher 
frequency,  the  mode  lines  formed  near  the  shoulder  of  the  bevel  where  the  tapered  edge 
begins.  Apparently  the  two  (1,1)  Chladni  patterns  are  the  result  of  this  relative  stiffness 
phenomenon. 
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A  comparison  of  the  first  mode  frequencies  from  the  modal  studies  with  the  frequency  sweeps 
shows  that  the  correlation  is  much  better  for  the  honeycomb  sandwich  panels  than  it  was  for 
the  skin-stringer  test  specimens.  Vibration  response  of  honeycomb  sandwich  is  less  sensitive 
to  ambient  temperature  changes  than  thin  skin-stringer  designs. 

b.  Frequency  Sweeps  -  Significant  strain  response  activity  above  the  first  mode  was  negii- 
gTETeT  hC-20A,  however,  showed  a  two-peak  response  with  the  frequencies  about  20  Hz 
apart  (see  Figure  45)  •  This  two-peak  type  response  was  detected  more  during  the  modal  fre¬ 
quency  studies  than  at  any  other  time. 

c.  Fatigue  Tests  -  Three  types  of  sonic  fatigue  failures  were  experienced  by  the  honeycomb 
sandwich  test  specimens.  Tne  most  common  was  fatigue  cracks  which  formed  around  the 
fastener  holes.  Figure  66  is  a  photograph  of  a  typical  fatigue  crack  of  this  type.  This  type 
of  failure  occurred  in  56.8  percent  of  tne  specimens. 

Facing  sheet  cracks  occurred  at  the  center  of  a  panel .  The  cracks  were  found  to  occur  for 
the  most  part  at  the  midspan  of  the  long  dimension  of  the  panel .  Figure  67  is  a  photograph 
of  this  type  failure.  This  type  of  failure  occurred  in  43.4  percent  of  ths  specimens. 

Core  shear  fatigue  was  the  third  type  of  failure  that  was  experienced.  This  type  occurred 
near  the  edge  (20-30%  of  short  span)  of  HC-15A  where  the  core  shearing  stress  was  the 
greatest.  Figures 68a  and  68b  are  photographs  of  this  core  shear  fatigue  failure. 

Bubbling-mode  failures  and  facing  sheet  to  core  bonding  failures  were  not  detected. 

As  discussed  previously,  the  test  spectrum  level  was  not  always  the  same  for  specimens  A 
and  B.  The  variation  seemed  to  be  a  function  of  test  specimen  frequency  and  test  specimen 
size.  However,  the  maximum  variation  in  level  was  3  db,  experienced  by  the  large  37- 
inch  by  61 -inch  panels .  One  to  two  db  variations  were  common . 

Like  the  skin-stringer  tests,  a  go/no-go  test  was  set  up  for  checking  the  honeycomb 
sandwich  design  nomograph  in  Reference  1  against  the  test  results.  Honeycomb  sandwich 
fatigue  data  from  Reference  21  were  used  to  establish  the  ±95-percent  confidence  limits. 

Three  of  the  thirty  designs  fell  in  the  go  category.  Ten  percent  of  all  specimens  tested  fell 
in  the  go  category.  Only  6.7  percent  of  all  the  specimens  tested  fell  below  the  lower  95- 
percent  confidence  limit  and  81 .7  percent  were  above  the  upper  95-percent  confidence 
limit,  indicating  that  the  design  nomograph  in  Reference  1  is  conservative  in  the  range  of 
test  conditions  reported  here. 

HC-16A  was  damaged  when  the  test  fixture  broke.  The  results  from  this  test  specimen  were 
omitted  from  all  analyses. 

HC-27  was  designed  with  an  aluminum  alloy  pan  0.025-inch  thick  instead  of  fiberglass. 
There  were  no  outstanding  increases  in  fatigue  resistance  or  response  characteristics  as 
compared  to  the  fiberglass  closure. 

HC-29  was  designed  with  an  aluminum  alloy  pan  0.030-inch  thick  and  a  tapered  doubler. 
Like  HC-27,  there  were  no  outstanding  changes  observed  in  fatigue  life  and  response 
characteristics  as  compared  to  the  conventional  doubler  and  fibeiglass  closure. 

HC-28  was  made  with  the  honeycomb  core  crushed  to  form  the  beveled  edge  instead  of  the 
conventional  machined  edge.  There  were  no  observable  changes  in  dynamic  response 
characteristics  or  differences  in  fatigue  life  as  compared  to  the  test  specimens  with 
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FIGURE  66.  TYPICAL  HONEYCOMB  SANDWICH  EDGE  FATIGUE  CRACKS 


FIGURE  67.  TYPICAL  HONEYCOMB  SANDWICH  FACING  SHEET  FATIGUE  CRACK 
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(a)  Typical  CORE  Shear  Fatigue  Failure 


(b)  Section  of  Honeycomb  Panel  Showing  Core  Shear  Fatigue 
FIGURE  68.  HONEYCOMB  CORE  SHEAR  FATIGUE  FAILURES 
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machined  beveled  edges.  The  test  specimens  were  sectioned  in  a  number  of  places  to 
examine  the  crushed  core  for  shear  fatigue.  The  examinations  showed  that  the  integrity  of 
the  core  had  not  been  altered  by  crushing  it. 
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IV  -  CORRELATION  OF  TEST  DATA  AND  DYNAMIC  STRESS  RESPONSE  PARAMETERS 


A.  Introduction 


The  theoretical  and  analytical  material  presented  in  Section  II  suggests  ways  of  correlating 
the  test  data  contained  in  Section  III  and  dynamic  stress  response  parameters  for  the  two 
types  of  aircraft  structure  relevant  to  this  research  program.  For  example,  the  first  mode 
response  was  observed  to  be  the  most  active  and  most  fatigue  damaging  for  both  types  of 
structure.  Also,  the  location  and  type  of  fatigue  cracks  experienced  in  the  experimental 
phase  of  the  program  suggest  there  are  bending  moments  at  the  edges  of  the  panels,  or 
along  the  fastener  line  in  the  case  of  the  skin-stringer  configuration.  Modal  frequency 
studies  and  frequency  sweep  investigations  indicate  strongly  that  both  the  skin-stringer  and 
honeycomb  sandwich  panel  edge  fixity  lies  intermediate  of  fully  clamped  and  simply 
supported.  See  Figures  AV-1  and  AV-2,  Appendix  V. 

It  is  the  purpose  of  this  section  to  show  the  correlation  of  test  data  with  the  dynamic  stress 
response  parameters  as  suggested  by  the  theory  presented  in  Section  III . 


B.  Skin-Stringer  Construction 


In  aircraft  structure,  the  individual  bays  (panels)  are  connected  to  adjacent  bays  or 
support  structure  at  their  boundaries  and  have  elastic  restraints  at  their  edges.  As  shown 
in  Section  II. D,  a  combination  of  stringer  torsion  and  bending,  rib  torsion  and  bending, 
and  phasing  among  adjacent  bays  can  result  in  a  very  complicated  dynamic  system.  For 
this  analysis,  however,  the  significant  mode  of  response  is  considered  to  be  that  of  a  plate 
having  an  edge  fixity  someplace  between  fully  clamped  and  simply  supported;  the  fixity 
to  be  determined  empirically. 

Intuitively,  Equation  (36)  contains  the  significant  parameters  for  describing  the  skin  stress 
response  of  a  bay  of  skin-stringer  structure.  This  will  be  tested  as  follows: 

First,  Equation  (36)  can  be  written  in  terms  of  root-mean-square  (rms)  stress  as 


H. .  „ 


(36a) 


where  Kj  is  a  proportionality  constant. 

Second,  the  first  mode  frequency  for  clamped  edges,  Equation  (20b),  can  be  expressed  as 


f.  = 

i 


i*,/2 


(20c) 


where  l<2  is  a  proportionality  constant  dependent  upon  the  edge  fixity. 

Third,  Equation  (20c)  is  substituted  into  Equation  (36a)  to  obtain  two  dimensionless  groups 
as  fo  I  lows: 
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(36b) 


where  is  a  proportionality  constant. 

Equation  (36b)  suggests  that  nominal  overall  stresses  measured  at  the  point  of  maximum 
stress  on  the  fastener  line  where  the  fatigue  cracks  developed  should  be  divided  by  the 

groupings  (c|E /y)  ^ ,  (t/ab)^2,  and  and  plotted  versus  the  dimensionless  group 

(b/t)2/(R8  451/72)  to  empirically  determine  K^.  This  was  accomplished.  Table  XII  is  the 
tabulated  results  and  Figure  69  is  a  rectangular  coordinate  plot  of  the  data. 

The  regression  lines  drawn  through  the  data  points  of  Figure  69  were  determined  by  the 
method  of  least  squares  described  in  Appendix  II.  A  correlation  coefficient,  a  measure  of 
the  "goodness  of  fit"  of  the  equation  to  the  test  data  was  computed  to  be  0.86  for  the 
linear  equation  and  0.88  for  the  exponential  equation.  From  a  table  of  correlation 
coefficients  in  Reference  22,  it  was  determined  that  the  probability  of  getting  correlation 
coefficients  of  the  magnitudes  described  above  in  the  absence  of  any  correlation  is  less 
than  1  in  1 ,000. 


The  exponential  regression  line  equation  for  the  data  plotted  in  Figure  69  is 


a  =  !  .62  x  10' 


ia 

Y 


ll/4a1.25*;f)  ^  % 

^■.75  ^0756  F(ar) 


where  F(ar)  -  (b/a) 1  ‘  75/R°  * 84  ■ 

For  aluminum  alloy,  (Eg/y)  '74pa443-and  Equation  (36c)  becomes 


a1-25*]/2 


(36c) 


(36d) 


The  regression  lines  from  the  physical  viewpoint  should  pass  through  the  origin  of  the  axes. 
It  did  not  in  the  case  of  the  linear  equation  because  of  test  data  scatter  and  possible 
response  non-linearities. 

The  arrangement  of  parameters  suggested  by  the  simply  supported  edge  equations.  Equations 
(24b)  and  (42),  was  aiso  analysed  as  described  above  because  the  edge  fixity  lies  between 
clamped  and  simply  supported.  The  correlation  coefficients  for  these  regression  lines  were 
significantly  lower:  0.54  for  the  linear  fit  and  0.58  for  the  exponential  fit. 

Comparing  the  correlation  coefficient,  it  is  apparent  that  Equation  (36d)  should  be  used 
to  develop  the  design  nomograph  for  aluminum  alloy  skin-stringer  piating. 
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TABLE  XII 


DATA  FOR  CORRELATING  SKIN-STRINGER  PANEL  RESPONSE  PARAMETERS 

WITH  MEASURED  STRESS 


Designation 

Dimensionless  group 

(t)2  „  ,0-3 

Dimensionless  group 

a  ,~-l 

\Y  /  \ab/  (f) 

^4775*1° 

1A 

15.25 

13.62 

B 

15.25 

14.36 

2A 

13.93 

11.82 

B 

13.93 

11.70 

3A 

15.79 

20.58 

B 

15.79 

21.30 

4A 

8.30 

5.74 

B 

8.30 

7.15 

5A 

4.66 

3.85 

B 

4.66 

4.10 

6A 

8.54 

7.50 

B 

8.54 

9.95 

7A 

26.26 

19.62 

B 

26.26 

16.43 

8A 

1.47 

1.07 

B 

1 .47 

1  .00 

9A 

20.64 

30.70 

B 

20.64 

32.28 

11A 

31 .67 

24.20 

B 

31 .67 

31 .64 

12A 

18.65 

12.55 

B 

18.65 

12.22 

13A 

16.30 

11.49 

B 

16.30 

12.03 

14A 

12.83 

7.30 

B 

12.83 

11.18 

16A 

6.19 

3.74 

B 

6.19 

5.38 

17A 

7.99 

4.25 

B 

7.99 

4.46 

18A 

4.66 

2.38 

B 

4.66 

3.31 
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TABLE  XII  (Continued) 


DATA  FOR  CORRELATING  SKIN-STRINGER  PANEL  RESPONSE  PARAMETERS 

WITH  MEASURED  STRESS 


Designation 

Dimensionless  group 

(b\2 

yj  . .  k  i 

Dimensionless  group 

0  irfl 

rstWVW2,'/ 2 

ly)  u)  % 

19A 

4.07 

2.41 

B 

4.07 

2.67 

20A 

5.75 

3.55 

B 

5.75 

3.16 

21A 

15.26 

7.89 

B 

15.26 

8.10 

22A 

4.48 

2.11 

B 

4.48 

2.67 

23A 

12.76 

8.48 

B 

12.76 

9.41 

24A 

21 .97 

17.92 

B 

21.97 

25.46 

25A 

3.49 

2.35 

B 

3.49 

2.26 

26A 

2.77 

1 .38 

B 

2.77 

2.47 

27A 

5.83 

2.08 

B 

5.83 

2.95 

28A 

5.91 

3.01 

B 

5.91 

3.01 

29A 

10.97 

6.70 

B 

10.97 

7.72 

30A 

11.93 

1.24 

B 

11.93 

1  .39 
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C.  Honeycomb  Sandwich  Construction 

Facing  sheet  and  edge  fatigue  cracks  are  by  far  the  most  common  failure  in  ho  jycomb 
sandwich  construction .  Core  shear  fatigue  is  considered  to  be  secondary.  All  three  of  the 
potential  failure  modes  were  investigated  from  the  test  data  correlation  point  of  view. 

I  .  Edge  Stress  Correlation 

The  observed  response  characteristics  of  the  flat  honeycomb  sandwich  panels  were  those  of 
a  panel  having  edge  fixity  between  fully  damped  and  simply  supported.  Therefore,  the 
arrangement  of  parameters  suggested  by  Equations  (82)  and  (86)  was  intuitively  believed  to 
be  the  arrangement  which  would  show  the  better  correlation.  These  equations  were 
arranged  into  two  dimensionless  groups  as  described  in  the  previous  section,  with  the  end 
result  as  follows: 


(86a) 


where  p  in  Equation  (82)  was  approximated  by  2ty/g,  ybeing  the  specific  weight  of  the 
facing  sheet  material. 


Table  XIII  is  a  tabulation  of  the  values  computed  as  suggested  by  Equation  (86a)  and 
Figure  70  is  a  plot  of  the  data.  The  regression  lines  shown  on  Figure  70  were  determined 
by  the  method  of  least  squares.  The  correlation  coefficients  for  the  linear  equation  and 
the  exponential  equation  are  0.86  and  0.88,  respectively. 

The  arrangement  of  dimensionlss  groups  suggested  by  the  combination  of  Equations  (83)  and 
(92)  for  simply  supported  edges  was  also  studied.  Correlation  coefficients  were  computed 
to  be  0.56  for  the  linear  fit  and  0.59  for  the  exponential  fit.  These  are  considerably  less 
significant  than  the  coefficients  for  fixed  edges. 


Consequently,  the  most  accurate  empirically  determined  regression  line  equation  for 
honeycomb  sandwich  edge  stress  design  is 


o  -  1 .46  x  10 


0.51 

a 

,h% 

1/2 

J.51 

r~ 

F(ar) 


(86b) 


where  F(ar)  =  (b/a)^  '^/R^ 

For  aluminum  alloy.  Equation  (86b)  becomes 


a  =  0.65 


0.51 

□ 

TTsT 


h  0, 


3? 


1/2 


F(ar) 


(86c) 


The  above  equation  will  be  used  in  Section  V  to  develop  the  design  nomograph. 


139 


TABLE  XIII 


DATA  FOR  CORRELATING  HONEYCOMB  SANDWICH  PANEL  EDGE  STRESS 
RESPONSE  PARAMETERS  WITH  MEASURED  STRESS 


Designation 

Dimensionless  group 
....  (w.  10-5 

Dimensionless  group 

li  ^  ,  r  ]  n"  ^ 

"3/4  '  1/2  *  10 

IA 

4.15 

2.65 

B 

4.15 

5.50 

2A 

5.63 

2.63 

B 

5.63 

2.22 

3A 

5.21 

4.30 

B 

5.21 

4.10 

4A 

2.27 

2.41 

B 

2.27 

2.92 

5A 

13.12 

8.06 

B 

13.12 

7.20 

6A 

8.00 

6.99 

B 

8.00 

6.48 

7A 

3.19 

1  .33 

B 

3.19 

1  .33 

9A 

3.03 

1.25 

B 

3.03 

1 .87 

10A 

2.93 

3.15 

B 

2.93 

3.10 

22A 

13.12 

7.24 

B 

13.12 

7.18 

24A 

10.00 

6.84 

B 

10.00 

6.56 

26A 

9.81 

3.84 

B 

9.81 

7.64 

27A 

4.50 

2.18 

B 

4.50 

2.31 

28A 

9.81 

3.55 

B 

9.81 

3.60 

29A 

2.79 

1.31 

B 

2.79 

2.25 

30A 

14.01 

7.43 

B 

14.01 

8.52 
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2.  Facing  Sheet  Stress  Correlation 


Equations  for  simply  supported  and  fully  clamped  edges  were  considered  in  the  correlation 
of  the  facing  sheet  stresses  and  flat  panel  stress  response  parameters.  The  arrangement  of 
parameters  suggested  by  Equations  (83)  and  (92)  and  arranged  in  two  dimensionless  groups 
as  described  previously  demonstrated  the  better  correlation  coefficients.  These  were  0.86 
for  the  linear  equation  and  0.87  for  the  exponential  equation.  Correlation  coefficients 
computed  for  the  arrangement  of  parameters  suggested  by  the  combination  of  Equations  (82) 
and  (89)  were  considerably  less  significant  than  those  listed  above. 


A  combination  of  Equations  (83)  and  (92)  produces  the  following  equation: 


(92a) 


where  Kq  is  a  proportionality  constant  and 


f(ar) 


Table  XIV  is  a  summary  of  the  values  computed  as  suggested  by  the  above  equation. 
Figure  71  is  a  plot  of  these  data  which  resulted  in  the  following  exponential  regression 
line  equation: 


arc  .1/4  0.46 

5=1.8x10-3|^1 


$1/2 

(f)  „  >0.73 

0.73,0.23  ■  XB  f(ar> 

t  h  5 


(92b) 


For  aluminum  alloy.  Equation  (92b)  becomes 


0.46  $ 


a  =  0.90 


0,50 


t0.73  h0.23  ^0.36 


f(ar) 


0.73 


(92c) 


Equation  (92c)  above  is  used  infection  V  to  develop  a  design  nomograph  for  honeycomb 
sandwich  facing  sheets. 


D.  Simple  Panel  Curvature  Effects 

Correlation  of  the  theory  developed  in  Section  II  for  evaluating  panel  curvature  effects  on 
stress  response  and  the  test  data  contained  in  Reference  12  will  be  described  in  this  section. 

The  specimens  tested  in  Reference  12  were  defined  by  the  parameters 
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TABLE  XIV 


DATA  FOR  CORRELATING  HONEYCOMB  SANDWICH  PANEL  FACING  SHEET 
RESPONSE  PARAMETERS  WITH  MEASURED  STRESS 


Designation 

Dimensionless  group 

Dimensionless  group 

&  in-1 

(fAy/2*;o/2~ 

11A 

2.62 

2.18 

B 

2.62 

2.00 

12A 

2.24 

1.44 

B 

2.24 

1.86 

13A 

2.32 

1.72 

B 

2.32 

1.45 

14A 

1 .43 

1.62 

B 

1 .43 

1.76 

15B 

6.05 

3.16 

16A 

3.77 

2.39 

B 

3.77 

2.35 

17A 

1.68 

1.26 

B 

1.68 

1.64 

18A 

1.64 

1.19 

B 

1.64 

1.19 

19A 

2.08 

1.27 

B 

2.08 

1.38 

20A 

3.08 

2.19 

B 

3.08 

2.06 

21 A 

5.36 

3.50 

B 

5.36 

3.24 

23A 

3.71 

2.39 

B 

3.71 

2.29 

25A 

3.78 

2.45 

B 

3.78 

1 .94 
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A  -  0.82 
L  =  344 

ft  -  0,  0.094,  0.125,  and  0. 188  radians 


Equation  (111),  rewritten  below,  was  used  to  evaluate  and  Cy-  The  specimen 
parameters  listed  above  and  ratios  of  rms  curved  panel  stress  to  rms  flat  panel  stress  shown 
in  Table  XV  were  used  to  accomplish  the  evaluation. 


ac 


=  |  1  + 

'  l 

The  stress  ratios  are: 


r  C2(ALft): 

1  1  L _ 4 


A4  +  0.61A2  +  1 


-3/4 


1  +  C, 


A2  4  0.034 
i  A4  +  9.62A2  +  1 


ALpf 


(Hla) 


TABLE  XV 

EXPERIMENTAL  CURVED  AND  FLAT  PANEL  RMS  STRESS  RATIOS 


0 

ac/af 

0.00 

1.00 

0.094 

0.80 

0.125 

0.45 

0.188 

0.23 

The  constants  C2  and  were  determined  to  be  0.006  and  0.453. 

Even  though  Equation  (111a)  is  not  exact,  a  method  is  available  for  predicting  the  effects 
or  radius  of  curvature,  aspect  ratio,  and  length  to  thickness  ratio.  Equation  (ilia)  will 
be  used  in  Section  V  to  develop  a  design  nomograph  for  assessing  curvature  effects  of  skin- 
stringer  panels. 


E.  Honeycomb  Sandwich  Curvature  Effects 

Strain  data  were  measured  on  one  curved  honeycomb  sandwich  test  specimen.  These  data 
were  used  to  check  (1)  the  ratio  of  the  facing  sheet  stresses  at  the  center  of  the  panel  and 
(2)  the  ratio  of  the  y  component  of  the  stress  to  the  x  component  of  the  stress  at  the  same 
point  on  the  curved  panel .  In  an  attempt  to  validate  the  effect  of  curvature  on  facing 
sheet  stress  at  the  center  of  the  panel  and  also  the  effect  on  edge  stress,  strain  data  were 
measured  on  a  similar  flat  sandwich  panel.  The  strains  were  transformed  by  the  theory 
developed  herein  to  accomplish  the  comparison. 

In  the  experimental  deteimination  of  strain  ratios,  both  flat  and  curved  panels  were  tested 
at  a  discrete  frequency  sound  pressure  level  of  130  decibels.  The  frequency  in  each  case 
was  set  at  the  (1,1)  mode  response.  Absolute  strain  was  not  determined;  however,  the 
strain  for  both  panels  was  normalized  to  the  same  base. 

The  flat  and  curved  panel  parameters  are  listed  in  Table  XVI .  The  non-dimensional 
parameters  are  given  in  Table  XVII,  and  the  relative  strain  data  are  tabulated  in 
Table  XVIII. 
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TABLE  XVI 


TEST  SPECIMEN  PHYSICAL  PARAMETERS 


Curved  Panel 

Flat  Panel 

a 

23.0" 

20.0" 

b 

29.0" 

30.0" 

nl 

.125" 

.  145" 

*2 

.008" 

.025" 

R 

84" 

00 

G 

y2 

9,000  psi 

15,000  psi 

G 

xz 

18,000  psi 

30,000  psi 

E 

16.2  x  106  psi 

10  x  10^  psi 

t 

A 

.045" 

.101"  -4  ply 

TABLE  XVII 


TEST  SPECIMEN  NON-DIMENSIONAL  PARAMETERS 


Parameter 

Curved  Panel 

Flat  Pune! 

0 

.345  Radians 

0 

L 

2875. 

800 

b/te 

511. 

198 

A 

1.26 

150 

g 

15.6 

5.8 

s 

.000495 

.00133 

c 

2.0 

2.0 

TABLE  XVIII 

RELATIVE  MEASURED  FACING  SHEET  STRAIN  (NO  SCALE) 
Outside  Curved  Panel  Flat  Panel 


f  y  center 
f  x  center 
*  y  edge 
f  x  edge 


34.2 

25.3 
38.0 
26.5 


14.5 

31.0 

15-5 

U1.0 


Inside  (z=-h  j) 

( y  center 
fx  center 


10.0 

13.5 


17.0 

27.5 


The  stress  values  for  the  flat  panel  stress  data  were  derived  using  Equation  (137)  for  °  y  with 
0  set  zero.  In  addition.  Equations  (149),  (155),  and  (156)  were  used  to  define  the  general¬ 
ized  displacement,  W  .  To  obtain  the  stress  transformation  between  flat  panel  designs  the 

equation  is 


A2L2S, 


a 


yl 


_L 


g  s 


center 


2 


y2 


center 


(158) 


for  converting  from  a  design  2  to  a  desired  design  1 .  A  similar  transformation  for  edge  stress 
was  derived  from  the  relationship 


a 

y 


Et 


2(1-  *2) 


,2 

a  w 

9y 4 


.  2 

+  ) 

dx 


(25a) 


3 

The  term - s-  is  zero  along  the  edge,  y  =  0. 

3  x 

used  to  determine  the  generalized  deflections. 


Again  Equations  (149),  (155),  and  (156)  were 
The  stress  transformation  ratio  for  the  edge  is. 
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=  > 


comparison  of  the  calculated  and  measured  data  is  presented  in  Table  XIX. 
tudes  do  not  correlate  exactly  the  trends  ore  apparent. 


TABLE  XIX 

STRESS  RATIO  COMPARISON 

Calculated  Measured 


Zxs. 

°xc 

center 

1.71 

1.16 

ffyc, 

z=hl 

1.86 

2.93 

ayc, 

Z=hl 

center 

°YC 

0.873 

0.337 

V 

center 

V 

edge 

2.04 

1.23 

(159) 


Although  mag- 


148 


V  -  DEVELOPMENT  OF  DESIGN  EQUATIONS  AND  NOMOGRAPHS 


A .  I  ntroduction 

The  information  required  to  produce  a  design  nomograph  for  a  specific  type  of  structure  is 
an  analytical  expression  for  the  stress  which  is  caused  by  the  random  acoustical  excitation 
and  the  allowable  fatigue  stress  and  life  relationship.  The  analytical  expression,  as 
derived  in  Section  IV,  has  empirically  determined  edge  fixity  and  built-in  stress  concen¬ 
tration  factors. 

To  be  useful  and  practical,  the  design  nomograph  also  must  be  as  free  as  possible  of 
complexities,  easy  to  use,  and  give  an  accurate  representation  of  the  test  results  from 
which  it  was  derived. 

B.  Skin-Stringer 

1  .  Skin  Design 

Figure  72  is  a  sketch  showing  the  nomenclature  for  a  simple  flat  panel  that  is  representative 
of  a  single  bay  of  skin  stringer  construction.  The  analytical  expression  used  for  developing 
the  skin-stringer  nomograph.  Figure  73,  is  Equation  (36d)  and  the  fatigue  curve.  Figure  34. 
An  example  problem  is  also  presented  to  illustrate  the  use  of  the  nomograph. 

2.  Stringer  Flange  Stresses 

The  prying  action  of  the  attached  flange  on  the  fasteners  (rivets)  causes  stresses  in  the 
attached  flange  which  may  be  critical .  Since  there  were  no  stringer  or  rib  flange  failure 
experienced  during  the  experimental  phase  of  this  program,  the  design  nomograph  is 
identical  to  the  one  developed  in  Reference  1  .  The  stress  in  the  flange  of  the  rib  was 
expressed  as; 


a 


r 


CC 


Pa 


r 


(160) 


One  factor  entering  the  rib-moment  equation  is  the  stringer  flange  width.  In  Reference  I 
it  was  stated  that  this  dimension  varies  only  with  stringer  thickness  for  good  design  practice 
and  is  evaluated  through  the  test  results.  Figure  74  is  a  sketch  of  skin-stringer  construc¬ 
tion  details  and  Figure  75  is  the  stringer  and  rib  flange  design  nomograph. 


3.  Design  Nomograph  for  Supporting  Structure 

In  order  to  develop  the  design  nomographs  for  stiffened  flat  panels,  it  was  convenient  to 
restrict  considerations  to  a  panel  with  constant  rib  spacing  and  constant  stringer  spacing 
covered  by  a  uniform  thin  sheet.  The  ribs  and  stringers  are  considered  to  be  uniform  mem¬ 
bers,  each  characterized  by  its  bending  rigidity  and  its  torsion  rigidity.  The  above  situa¬ 
tion  is  illustrated  in  Figure  76. 

As  mentioned  previously,  there  are  three  modes  of  interest  for  sonic  fatigue  considerations. 
These  modes  are  broadly  characterized  as  rib  bending  -  stringer  trosion,  rib  -  torsion  - 
stringer  bending,  and  coupled  rib  -  stringer  bending  -  torsion.  Hence,  the  nomograph 
can  be  used  to  predict  three  frequencies  for  each  design.  The  purpose  of  this  analysis  is 
to  consider  the  motion  of  the  structure  as  a  whole  and  is  not  concerned  with  component 
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FIGURE  72.  NOMENCLATURE  FOR  SIMPLE  FLAT  PANEL 
(REPRESENTATIVE  OF  A  SINGLE  BAY) 


PANEL  WIDTH,  a,  INCHES 


[O  FAILURE 


73.  SKIN-STRINGER  PLATING  DESIGN  NOMOGRAPH 
7075-T6  ALUMINUM  ALLOY 


EXAMPLE:  A  flat  aluminum-alloy,  skin-stringer  structure  is  required  to  withstand  an  esti¬ 
mated  service  noise  spectrum  level  of  120  db.  The  design  life  is  5  x  10^  cycles,  the 
damping  ratio  is  assumed  to  be  0.012,  the  assumed  stringer  spacing  a  =  4.75",  and  the 
aspect  ratio  is  assumed  to  be  1  .5.  Follow  through  the  nomograph.  Figure  73,  as  indicated 
by  the  arrows  and  the  skin  thickness  t  =  0.032".  Next,  follow  through  the  stringer  or  rib 
flange  nomograph,  Figure  75,  as  indicated  by  the  arrows  to  tr«  0.043".  The  fundamental 
frequency  can  be  estimated  using  Figure  AVI-1,  Appendix  VI,  and  in  this  case  is  calcu¬ 
lated  to  be  290  Hz.  At  this  frequency,  the  service  environment  noise  spectrum  level  is 
checked  with  the  noise  spectrum  level  used  above.  If  necessary,  an  iteration  is  made  to 
obtain  agreement. 

If  other  materials  such  as  titanium,  etc.,  are  to  be  considered.  Figure  48,  Reference  1, 
should  be  used  to  make  the  necessary  fatigue  curve  conversion. 


FIGURE  74.  SKETCH  OF  SKIN-STRINGER  CONSTRUCTION  DETAILS 


Blank 


153  preceding  Pa9e 


7075-T6  AND  2024-T3 


FIGURE  75.  STRINGER  OR  RIB  FLANGE  DESIGN  NOMOGRAPH 


modes,  such  as  individual  panel  motion.  If  the  cover  sheet  is  thin,  then  the  sheet  affects 
the  frequency  mainly  by  its  mass  effect.  Thus,  the  stiffness  of  the  cover  sheet  will  be 
neglected  as  small  compared  to  the  stiffness  of  the  supporting  structure. 

From  Equations  (55a)  through  (55f),  i.t  is  seen  that,  if  terms  on  the  order  of  /3  L  and  /3  are 
neglected  as  small  compared  to  /3L  ,  the  stiffness  matrix  reduces  to  consideration  of  only 
[^xx  3  and[Kzz*J  Similarly,  since  the  rib  and  stringer  spacing  is  considered  to  be  uniform, 
the  mass  matrix  uncouples  so  that  the  general  eigenvalue  problem  can  be  reduced  by  con¬ 
sidering  only  the  diagonal  submatrices: 


'  Kzz  ' 
K 

yy 


Equation  (161)  governs  the  rib-bending  stringer- torsion  mode.  Equation  (162)  governs  the 
rib  -torsion  stringer- bending  mode,  and  Equation  (163)  governs  the  coupled  rib-stringer 
bending-torsion  mode. 

The  lowest  eigenvalue  of  each  of  the  above  problems  can  be  approximated  as: 

"2  [  Kn  -  I"  +  Ki3 1  V2  ]  <,64> 


where  Kj,,  Kj2,  K^,  and  |  depend  upon  the  mode  under  consideration, 
bending  stringer  -torsion  mode,  the  constants  are 


For  the  rib— 


Kn  =  2.0  (GJ/L)$+  8.0  (EI/L)r 

K12  =  (GJ/L)S  (165) 

K13  =  2.0  (EI/L)R 

M11  =  ^lr  "  lr  Ls  +  Ls  ^  (2*°  wplr ls  +  W579 

For  the  rib- torsion  stringer-bending  mode  the  constants  are 
Kn  =  2.0  (GJ/L)r  +  8.0  (EJ/L)$ 

K.?=  2.0  (El/L)c 

U  1  (166) 
K,3  -  (GJ/L)r 

Mn=  <lr2  -  44  +  LS>  <2’0  wp44  +  ws4>/579- 


156 


Finally,  for  the  coupled  pb-sfringer  bending-torsion  mode,  the  constants  are 
Kn  =  24.  j  (El/i.3)s  +  (EI/L3)r  | 

K)2  =  12  (EI/L3)s  (,67) 

K13=  12  <EI/L3)r 

M,,-  (2.0  W?LRLS  +  WRLR  +  WsLs)/,93 

Equation  (164)  is  evaluated  for  the  appropriate  K^,  ^12'  ^13'  anc^  ^11  ^  use  ^'9ures 
77  through  78  . 


As  an  example,  consider  a  panel  with  the  following  characteristics 


(ei)r 

=  9; 

x  104  lb.  in.2 

<EI)S 

=  4.5  X  104  lb.  in.2 

(GJ)r 

-•  1 

.5  x  102  lb.  in.  2 

(GJ)S 

=  2.0  x  102  lb.  in. 2 

73 

il 

0.1 

lb. in. 

II 

U1 

£ 

0.05  lb. in. 

r~ 

73 

11 

6  in. 

Ls  = 

8  in. 

^skin 

=  0. 

032  in. 

Wp  = 

0.0032  Ib/in. 2 

The  frequency  for  the  rib-bending  stringer- torsion  mode  is  computed  as  follows: 

1  .  Enter  Figure  77  with  (GJ)^  =  2.0  x  102,  move  to  the  left  until  L<.  =  8.0,  and  pro¬ 
ject  up  to  read  =  25.0.  At  (El)p  =  9x  104  move  down  to  Lp  =  6.0  and  to  the 
right  to  read  =  3.0  x  10  .  Projecting  the  line  downward  and  the  line  to 

the  left,  Kj  ^  =  1 . 1  x  10^  is  read  in  the  lower  left  hand  corner  of  the  chart. 


2. 


Enter  Figure  78  at  the  lower  left  hand  s,de  with  =  Lp  =  6,  project  to  the  right  until 
1-2  =  L<j  =  8  is  reached,  move  vertically  to  Wp  =  0.0032,  and  then  project  to  the  right. 
With  L|  =  6  enter  the  lower  right  hand  side  of  the  chart,  move  left  until  =  Wp  =0.1 
is  reached,  and  project  upward.  At  the  intersection  point  read  the  value  of  Q.  Enter 
Q  in  the  lower  right  hand  side  of  the  chart,  move  vertically  to  the  line  P  =  52,  and 


project  to  the  right  to  read  ^  =  0.09. 

Enter  Figure  79  in  the  left  hand  side  with  K,,  =  1 . 1  x  10^  and  project  horizontally. 

2  2  8 

Movinq  downward  from  +  KT-6  =  9  x  10  ,  read  the  value  at  the  intersection  point 
r~H  I  ^  14.  / — x — - «- 


of  Kj  j  -  +  Kjg  =  1  x  IQ5.  On  the  right  hand  side  enter  -  y  ^]2  f  ^ 

move  left  to  the  Mn  =0.09  line,  and  project  vertically  to  read  f  =  180  Hz. 


The  frequency  for  the  rib-torsion  stringer-bending  mode  is  computed  in  a  similar  fashion  with 
only  slight  changes.  Enter  Figure  77  with  (GJ)R  and  (El)<j  and  proceed  as  before  to  read 


K| 2'  ^13/  and  (follow  the  dashed  line).  Enter  Figure  78 
=  Wj  and  read  ^  .  Enter  Figure  79  with  ^  and  + 

as  before. 


with  L1  =  Ls,  i-2  =  Lr,  W1 
and  compute  the  frequency 


The  frequency  forcoupled  rib-stringer  bending-torsion  mode  iscomputed  using  Figure  80  through 
82.  Enter  Figure  80  with  (EI)R  and  LR  as  the  scales  to  the  right.  Connect  (EI)R  and  LR  to 
read  12(EI/L^)R  =  Similarly,  with  (El)<-  and  L^,  read  12(EI/L^  =  K^-  In  the  left 

hand  portion  of  Figure  80  ,  enter  and  and  at  the  intersection  point  read  ^  .  On 
Figure  81 ,  enter  the  three  scales  on  the  left  by  connecting  WR  with  LR  and  with  to 
read  the  numbers  WRLR  and  W<jLg.  In  the  center  portion  connect  LR  and  to  read  LrL<j. 


Connect  bRL^  to  Wp  to  read  2  WpLRLg.  Enter  2  WpLRL<.  +  WRLR  +  Wj-L^  on  the  scale  to 
the  right  and  move  laterally  to  read  Mj  ^ .  Finally,  enter  Figure  82  with  and  + 
as  in  Figure  79.  Enter  K  ]]~  yj  ^]2  +  ^3  =  7.8  x  10^  on  the  right  side,  move 
laterally  to  =  6.7  x  10  and  project  vertically  to  read  f  -  160  Hz. 


C.  Honeycomb  Sandwich 


1 .  Edge  Design 

The  mathematical  expression  used  for  estimating  the  honeycomb  sandwich  edge  stress  is 
Equation  (86c).  It  was  used  in  combination  with  the  appropriate  fatigue  curve.  Figure  50, 
to  derive  the  honeycomb  sandwich  edge  design  nomograph.  Figure  83  is  the  design  nomo¬ 
graph  and  Figure  84  is  a  drawing  of  the  construction  details  with  an  example  problem  solved. 


2.  Facing  Sheet  Design 

Equation  (92c)  was  used  to  derive  the  honeycomb  sandwich  facing  sheet  design  nomograph. 
Figure  51,  the  facing  sheet  fatigue  curve,  was  used  to  relate  facing  sheet  stresses  to  number 
of  positive  crossings  to  failure.  Figure  85  is  the  resulting  honeycomb  sandwich  facing  sheet 
design  nomograph.  Figure  84  is  a  sketch  of  the  construction  details  with  an  example  problem 
solved . 


3.  Core  Shear  Design 

The  maximum  shear  stress  in  the  core  of  a  flat  honeycomb  sandwich  panel  occurs  in  the 
neighborhood  of  the  edges  regardless  whether  the  panel  is  simply  supported  or  clamped. 
The  mathematics  expressing  the  core  shear  for  simply  supported  edges  is  less  complicated 
and  slightly  conservative  compared  to  clamped  edges  and  can  be  used  for  purposes  of 
checking  the  core  shear  fatigue  resistance. 
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FIGURE  77.  NOMOGRAPH  FOR  COMPUTING  <1]f  K]2,  AND  K]3  FOR  RIB  BENDING-STRINGER  TORSION  AND 
RIB-TORSION  STRINGER  BENDING  MODES  OF  SKIN-STRINGER  PANELS 


RiB-TORSION  STRINGER-BENDING  MODE 


RIB-TORSION  STRINGER-BENDING  MODES 


FIGURE  82.  NOMOGRAPH  FOR  COMPUTING  FREQUENCY  OF  COUPLED  RIB-STRINGER  BENDING  TORSION  MODE 


EXAMPLE:  A  flat  aluminum-alloy/  honeycomb  sandwich  structure  is  required  to  withstand 
an  estimated  service  environment  noise  spectrum  level  of  130  db.  The  design  life  is 
5  x  10®  cycles,  the  average  damping  is  considered  to  be  0.019,  the  assumed  panel  width 
a  -  20",  an  aspect  ratio  of  1 .2,  and  a  honeycomb  core  depth  h  =  0.45".  Follow  through 
Figure  03  as  indicated  by  the  arrows  and  obtain  an  edge  thickness,  t~  =  0.064".  Move  to 
Figure  85,  using  the  same  design  conditions  specified  above,  and  follow  the  arrows  to 
obtain  the  required  facing  sheet  thickness,  t  =  0.015".  The  fundamental  frequency  can  be 
calculated  using  Figure  AVI-2,  Appendix  VI,  and  is  found  to  be  290  Hz.  At  this  frequency 
the  service  environment  noise  spectrum  level  is  checked  with  the  noise  spectrum  level  used 
above.  If  necessary,  an  iteration  is  made  to  obtain  agreement. 

If  other  materials  such  as  titanium,  etc.,  are  to  be  considered.  Figure  48,  Reference  1, 
should  be  used  to  make  the  necessary  fatigue  curve  conversion. 


FIGURE  84.  HONEYCOMB  SANDWICH  CONSTRUCTION  DETAILS 
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CYCLES  TO  FAILURE 
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The  shearing  stress  in  the  core  of  a  flat  honeycomb  sandwich  panel  subjected  to  static  load 
can  be  approximated  by  the  following  equation,  derived  in  Reference  17: 
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W  D  n3 
o 

h  a 


1  +  1 


(168) 


where  Wq  is  evaluated  by  Equation  (40). 

The  maximum  static  shear  stress  resulting  from  a  uniform  pressure  of  unit  magnitude  is 


T 


16  1 


(169) 


The  rms  core  shear  stress  for  a  lightly  damped,  simply  supported  flat  honeycomb  sandwich 
panel  exposed  to  random  acoustical  excitation  can  be  estimated  by  combining  Equations 
(5),  (83),  and  (169).  The  result  is 

1/2 


0(f) 

hS 
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(170) 


where  the  damping  was  assumed  to  be  an  average  value  of  1 .85%.  (See  Table  IX.) 

The  required  core  density  (lbs/ft^)  can  be  estimated  using  Equation  (170)  and  the  core  shear 
fatigue  curve  extracted  from  Reference  23.  This  fatigue  curve  was  derived  from  flexural 
tests  of  bonded  sandwich  panels  subjected  to  bending  conditions  producing  high  core  shear 
and  facing  sheet  stresses.  The  constant  amplitude  core  shear  stresses  were  converted  into 
rms  core  shear  stresses  using  the  method  recommended  by  Crede  and  Lunny  (Reference  24). 
This  simple  relationship  is 


(171) 


where  x  is  the  rms  shear  stress. 

x  is  the  constant  amplitude  shear  stress. 

a  is  the  inverse  negative  slope  of  the  fatigue  curve  equal  to  6.2. 
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There  are  only  two  core  shear  fatigue  points  available  to  check  the  accuracy  of  the  fatigue 
curve:  HC-15A  and  one  from  a  recently  conducted  C-5A  sonic  fatigue  test.  (The  latter  is 
reported  in  Reference  25.)  The  core  shear  for  each  panel  was  estimated  using  Equation 
(170)  and  the  parameters  listed  in  the  table  below: 


TABLE  XX 

CORE  SHEAR  STRESS  PARAMETERS 


Designation 

a 

b 

h 

SPL 

db 

Core 

density 

Damping 
ratio,  % 

HC-15A 

20 

24 

0.27 

140 

3.4 

2.0 

C-5A  panel 

23 

35 

0.50 

139 

2.4 

2.2 

The  rms  shear  stress  for  HC-15A  was  computed  to  be  60.1  psi  and  it  developed  a  core  shear 
failure  in  1 .08  million  zero  crossings.  The  rms  shear  stress  for  the  C-5A  panel  was  com¬ 
puted  to  be  35.1  psi  and  it  developed  a  shear  failure  in  2.1  million  zero  crossings.  These 
points  are  plotted  on  Figure  87  to  show  a  comparison  of  predicted  life  and  test  lift. 

Equation  170  and  its  companion  fatigue  curve.  Figure  86,  can  be  used  to  select  the 
required  core  density  (Ib/ft3)  for  a  given  design. 


D.  Skin-Stringer  Curvature  Effects 

The  effect  of  curvature  on  stress  ratio,  expressed  by  Equation  (1 1  la).  Section  V-D,  was 
used  to  develop  the  curvature  effects  nomograph  for  skin-stringer  construction. 

The  curvature  effects  nomograph  is  designed  to  be  used  in  connection  with  the  design 
nomograph  for  flat  skin-stringer  plating.  Figure  73  is  used  to  estimate  the  panel  dimen¬ 
sions.  Using  these  dimensions,  a  stress  reduction  ratio  due  to  curvature  is  obtained  from 
Figure  87.  The  procedure  is  repeated  until  a  satisfactory  convergence  is  obtained.  Two 
to  three  iterations  are  generally  sufficient  for  satisfactory  convergence. 


E.  Design  Equations  for  Curved  Honeycomb  Panels 

The  equations  below  were  developed  as  a  design  aid  for  curved  honeycomb  panels.  Although 
the  equations  are  accurate  within  the  assumptions,  it  was  deemed  impractical  to  design 
nomographs  for  solution  and/or  evaluation  of  the  equations  presented. 


To  determine  the  effect  of  curvature  on  face  sheet  (z  -  h^)  tensile  stress  at  the  center  of  the 
panel,  in  the  direction  of  curvature  (y),  use  the  following  equation: 
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DATA  POINT 


FIGURE  86.  5052-H39  ALUMINUM  ALLOY  CORE  SHEAR  FATIGUE  CURVE 


The  effect  of  curvature  on  face  sheet  (z  =  hj)  tensile  stress  in  the  longitudinal  direction  (x), 
at  the  center  of  the  panel,  is  determined  by  the  relationship  which  follows: 


1  +  .0396  (1  +  .103  A^) 


'i  (6te\ 1  ^ 
\  9  /  J  Wf 


The  ratio  of  a*  to  at  the  center  of  the  panel  in  the  outer  face  sheet  (z  =  hj)  is  found 
from  the  following  relationship,  which  holds  for  flat  as  well  as  curved  panels. 
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The  ratio  of  tensile  stress  at  the  panel  center  in  the  outer  face  sheet  (z  =  h,)  to  tensile  stress 
in  the  inner  face  sheet  (z  -  -h^)  is 
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The  effect  of  curvature  on  edge  stress  at  the  center  of  straight  edge  (y  =  0  or  a,  x  -  b/2) 
is  given  by  the  following  equation: 


1  +  .445  (Az  +  .108)  ALg0  Wc 

?  A  W 

1+9.62AZ  +  A*  vvf 


(176) 


The  effect  of  curvature  on  transverse  shear  stress  is  simply  stated  as  the  ratio  of  generalized 
coordinates 
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If  the  excitation  is  broad-band  random  noise,  the  ratio  Wc/W^  is 
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1 

s2 
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If  the  excitation 

is  sinusoidal. 

the  ratio  becomes 
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Wf  = 
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It  is  assumed  in  Equations  (172),  (173),  (176),  and  (1 77)  that  the  damping  ratio  for  the  first 
mode  does  not  change  as  curvature  increases  (frequency  increases).  If  this  assumption  is  not 
acceptable  for  certain  designs,  then  the  ratio  of  W^/W^.  is  expressed  as 

(a)  for  random  excitation 


(b)  for  sinusoidal  excitation 


(180) 


(181) 


The  parameters  to  are  determined  by  the  following  relationships 


S 

S 


1 

2 


S 

S 


3 

4 


1  +  A4  +  40.7  l- 

)  (1  +  cA2) 

\  scL^  / 

1  +  40.7/ 

-Jy  ) 

(]+?-)+  1655  /  — ' 

\ 

scL^  / 

\  A2  /  \  sc 

1  +  .33A2 

+  38.0 

(  0  +  -143  c) 

\  scL^  / 

1  +  3.0A2 

+  16.3 

(*)  (,  +  7-0c) 

(182) 


176 


Due  to  simplifying  assumptions  made  in  the  derivation  of  these  equations,  certain  restrictions 
are  placed  on  the  variables.  These  limits  are 

0  <  .35 

g  >  5 

s  <  .01 

L  >  100 

.3  <  A  <  3 

L  >  100 

e 

The  equations  derived  in  this  section  are  too  complicated  for  graphical  solution.  If 
extensive  calculations  are  to  be  made,  the  computer  program  described  in  detail  in 
Appendix  IV-B  can  be  used.  While  the  program  is  coded  in  Fortran  IV,  slight  modifica¬ 
tions  may  be  necessary  for  compatibility  with  various  computers.  This  specific  program 
was  developed  for  use  on  the  IBM  360  system . 
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VI.  CONCLUSIONS 


A.  The  results  of  the  analytical  and  experimental  investigations  reported  here 
indicate  that  the  overall  rms  stresses  induced  in  flat  skin-stringer  plating  and  honey¬ 
comb  sandwich  construction  by  broad-band  noise  can  be  estimated  with  sufficient 
accuracy  for  design  purposes  by  a  development  involving  plate  theory  and  Miles' 
response  theory , 

B.  The  skin-stringer  plating  and  honeycomb  sandwich  design  nomographs  presented  in  this 
report  give  estimates  of  fatigue  life  well  within  the  ±95-percent  confidence  limits  of 
their  respective  fatigue  curves.  Obviously  the  accuracy  of  the  nomographs  can  be 
established  only  in  the  range  of  the  test  conditions.  Additional  testing  is  required  to 
define  the  accuracy  of  the  extrapolated  portion  of  the  fatigue  curves. 

C.  Skin-stringer  and  honeycomb  sandwich  construction  designed  using  the  nomographs 
and  procedures  presented  here  will  be  lightweight  and  sonic-fatigue-resistant. 

D.  Back-up  structure,  ribs,  and  stringers,  designed  using  the  nomographs  and  procedures 
described  in  this  report,  will  be  conservative.  None  of  the  30  skin-stringer  designs 
that  were  tested  experienced  rib  or  stringer  flange  cracks. 

E.  The  dynamic  analysis  of  back-up  structure  indicates  that  the  motion  is  highly  coupled 
and  that  the  phase  relations  between  adjacent  bays  is  extremely  important  for  fatigue 
considerations.  Hence,  in  order  to  realize  the  degree  of  constraint  along  a  boundary 
of  an  individual  panel,  the  coupled  effect  cf  rib  and  stringer  motion  must  be  taken 
into  account.  This  approach  is  necessary  since  the  back-up  structure  represents  some 
60  percent  of  the  weight  of  skin-stringer  construction.  The  back-up  structure  must 
then  be  considered  as  an  integral  part  of  the  fatigue  analysis  in  order  to  design  a 
system  which  is  optimum  from  a  weight  standpoint.  In  addition,  consideration  must  be 
given  to  multi-bay  effects  for  structure  which  deviates  from  the  model  considered  here. 

F.  The  experimental  data  contained  in  Reference  12,  when  compared  with  the  various 
analytical  studies  reported  here,  indicate  that  aircraft-type,  skin-stringer  construc¬ 
tion  does  not  realize  the  degree  of  in-plane  restraint  represented  in  the  classical 
clamped  panel  analysis.  An  analytical  solution  that  more  nearly  determines  the 
fundamental  frequency  of  an  aircraft-type  curved  panel  can  be  obtained  by  relaxing 
the  in-plane  restraint  along  the  length  through  appropriate  changes  in  the  energy 
representation  and  the  assumed  eigenmodes. 

G.  The  effects  of  curvature  in  honeycomb  sandwich  construction  evaluated  by  state-of- 
the-art  analysis  techniques  are  (1)  the  vibratory  stresses  in  the  outer  facing  sheet,  in 
the  direction  of  curvature,  are  of  greater  magnitude  than  the  stresses  in  the  inner 
facing  sheet,  when  the  facing  sheets  ere  of  equal  thickness;  (2)  the  vibratory  facing 
sheet  stress  at  the  center  of  a  panel  increases  with  curvature,  because  of  membrane 
influences;  and  (3)  honeycomb  core  shear  stress  decreases  with  panel  curvature. 

H.  The  design  nomographs  and  computer  program  for  skin-stringer  and  honeycomb  sand¬ 
wich  construction  which  account  for  structural  curvature  effects  will  increase  sonic 
fatigue  design  capabilities. 
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APPENDIX  I 


TEST  DETAILS  AND  INSTRUMENTATION 
A-  Modal  Frequency  Investigation 


The  test  arrangement  for  the  modal  frequency  investigation  is  shown  in  Figure  A I  - 1  .  This 
was  accomplished  by  placing  the  frame-mounted  specimens  over  two  15-inch  loudspeakers, 
one  of  which  was  equipped  with  a  phase-reversal  switch  so  that  the  test  specimen  could  be 
excited  with  an  in-phase  or  out-of-phase  forcing  function.  The  mode  shapes  (Chladni 
patterns)  were  obtained  by  sprinkling  cork  particles  on  the  surface  of  the  vibrating  panel, 
and  then  varying  the  input  frequency  until  the  cork  was  aligned  along  the  mode  lines. 
When  a  well-defined  pattern  v/as  observed,  the  input  frequency  was  recorded.  After  a 
study  of  the  mode  shapes,  it  was  decided  that  the  strain  gages  should  be  located  in  the 
positions  shown  in  Figures  22  and  23,  Section  III. 

B.  Damping  Studies 


The  test  setup  for  performing  damping  studies  is  shown  in  Figure  AI-2.  The  test  specimens 
were  mounted  over  two  loudspeakers  in  the  same  manner  as  was  done  for  the  modal  frequency 
investigation.  An  oscilloscope  and  a  Polaroid  Oscilloscope  Camera  was  used  to  photograph 
the  decaying  strain  oscillations.  The  photographs  were  obtained  by  using  a  3-way  switch 
which  turned  off  the  loudspeaker  and  at  the  same  time  triggered  the  oscilloscope.  When  the 
trace  appeared  on  the  oscilloscope  a  photograph  was  made. 


The  logarithmic  decrement  was  found  using  the  following  equation 


In  X  /X 
o 


n 


where  D  =  the  logarithmic  decrement 

Xq  =  the  amplitude  of  the  damped  wave  at  point  0 
X^  =  the  amplitude  of  the  damped  wave  after  n  cycles 
n  =  the  number  of  cycles . 

The  damping  ratio  5  was  determined  using  the  log  decrement  and  the  equation 
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In  order  to  get  an  average  value  for  each  panel,  several  decaying  strain  oscillaTions  were 
analyzed.  The  mean  damping  ratio  was  calculated  using 


n 


1 


where  <5^  =  the  average  damping  ratio 

N  =  the  number  of  decoying  strain  oscillations  taken  into  account  for  a 
specific  panel . 
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TEST  SPECIMEN 
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FIGURE  AM.  TEST  SET-UP  FOR  MODAL  STUDIES 


C  .  Discrete  Frequency  Sweeps 


The  test  arrangement  for  conducting  the  sinusoidal  frequency  sweeps  is  shown  in  Figure  AI-3. 
The  test  specimens  were  instrumented  with  strain  gages  and  mounted  in  the  test  facility. 

Each  strain  gage  was  connected  in  a  four-arm  bridge  circuit  and  subsequently  connected  to 
a  6-channel  signal  conditioning  unit.  The  test  panels  were  then  excited  by  using  two  trans¬ 
ducers  for  a  frequency  range  of  80  Hz  to  2000  Hz.  The  strain  gage  signal  was  amplified  and 
then  analyzed  using  a  Spectral  Dynamics  Analyzer  with  a  narrow  band  (2  Hz)  filter.  The 
output  was  plotted  on  an  x-y  recorder. 

During  the  frequency  sweeps,  accelerometers  mounted  on  the  test  fixture  were  monitored 
to  see  if  test  fixture  resonance  was  interfering  with  the  test  specimen  response.  There  was 
no  indication  of  significant  test  frame  resonances  in  the  frequency  range  of  80  Hz  to  2000  Hz . 

D.  Panel  Response 

The  strain  gaged  specimens  were  mounted  in  the  test  facility  as  shown  in  Figure  25,  Section 
III,  and  exposed  to  broadband  noise  excitation.  The  response  of  the  panels  were  obtained 
by  using  the  same  procedure  that  was  used  for  the  frequency  sweeps.  A  simplified  diagram 
of  the  set-up  is  shown  in  Figure  AI-4.  From  the  response  plots  the  most  active  gage  was 
determined  for  each  set  of  test  specimens.  This  gage  was  then  recorded  on  magnetic  tape 
using  the  set-up  shown  in  Figure  AI-5.  Before  each  recording  was  made,  a  calibration 
signal  from  the  strain  gage  was  recorded  on  tape. 

A  tape  loop,  which  included  the  calibration  signal,  was  made  from  the  recorded  data. 

The  loop  was  used  to  obtain  PSD  and  peak  strain  probability  analyses.  The  data  analysis 
systems  are  shown  in  Figure  AI-6,  and  a  block  diagram  of  the  set-up  is  shown  in  Figure 
AI-7. 


F. .  Spectrum  Set-up 

The  block  diagram  for  plotting  of  the  spectrum  is  shown  in  Figure  AI-8.  The  spectrum  for 
each  set  of  panels  was  shaped  to  concentrate  the  acoustical  energy  in  a  desired  range  of 
frequencies.  The  bandwidth  of  the  test  spectrum  was  chosen  to  be  at  least  three  times  the 
bandwidth  of  the  panel  response. 

F.  Calibration  of  Test  Section 

A  detailed  survey  of  the  sound  pressure  level  distribution  in  the  test  section  of  the  High 
Intensity  Structural  Test  Facility  was  made.  Discrete  frequency  sweeps  were  used  to  make 
the  evaluation.  Discrete  frequency  sound  pressure  levels  were  recorded  at  the  locations 
shown  on  Figure  AI-9.  This  figure  is  also  a  plot  of  the  sound  pressure  level  distribution 
for  145  decibels  reference  level  and  frequencies  of  125  Hz,  150  Hz,  and  275  Hz.  This 
range  of  frequencies  encompasses  the  (1,1)  mode  response  of  most  of  the  specimens  tested 
for  this  program.  It  can  be  seen  that  the  distribution  of  sound  pressure  level  over  the  test 
section  at  125  Hz  is  uniform  with  a  maximum  positive  deviation  of  3  db  and  a  maximum 
negative  deviation  of  1  db  relative  to  the  reference  145  db  at  the  center  of  the  test 
section.  The  distribution  is  somewhat  better  at  150  Hz  with  a  maximum  positive  deviation 
of  3  db  and  no  negative  deviation.  The  distribution  at  275  Hz  showed  a  maximum  positive 
deviation  of  2  db  and  a  maximum  negative  deviation  of  3  db. 

G  .  Microphone  Calibration 

The  microphones  were  calibrated  before  each  test  spectrum  was  shaped.  This  was  accom¬ 
plished  by  using  a  Photocon  Model  PC-125  calibrator,  which  is  traceable  to  a  secondary 
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FIGURE  AI-3.  BLOCK  DIAGRAM  OF  FREQUENCY  SWEEP  TEST  SET-UP 


FIGURE  AI-5.  BLOCK  DIAGRAM  OF  TEST  SET-UP  FOR  DATA  RECORDING 
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a.  Block  Diagram  of  PSD  Analysis 


b.  Block  Diagram  of  Peak  Strain  Probability  Analysis 

FIGURE  AI-7.  BLOCK  DIAGRAMS  FOR  STATISTICAL  ANALYSIS 
OF  TEST  DATA 
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X-Y  RECORDER 
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FIGURE  AI-8 .  BLOCK  DIAGRAM  OF  TEST  SET-UP  FOR  PLOTTING  SPECTRUM 


FIGURE  AI-9.  SOUND  PRESSURE  LEVEL  DISTRIBUTION  IN  TEST  SECTION 


standard.  The  signal  from  the  microphone,  proportional  to  the  calibration  sound  pressure 
level,  was  sent  though  the  analyzing  equipment  to  the  x-y  recorder.  The  displacement 
of  the  recorder  pen  was  equal  to  the  level  applied  to  the  microphone. 

When  the  same  spectrum  was  used  for  several  days  of  testing,  the  microphones  were  cali¬ 
brated  at  the  beginning  and  the  end  of  each  day  to  see  if  there  was  any  change  in  the  sen¬ 
sitivity. 


H.  Strain  Gage  Curcuit 

A  four  arm  bridge  circuit  shown  in  Figure  Ai  —  10  was  used  for  all  strain  measurements.  The 
mathematical  expression  for  this  bridge  circuit  is  as  follows: 


E  = 


VB  RG 
RG  +  R2 


where  Vg  =  excitation  voltage 

Rq  =  strain  gage  resistance  (  120  ohm) 

R^.  =  dummy  gage  120  ohm  (for  temperature  compensation  ) 
Rj  =  R2  =  fixed  resistance  (  120  ohm) 

F  =  gage  factor  (F  =  2) 

Eo  =  gage  output  (  volts) 


Differentiating  with  respect  to  Rq 

dEo  _(RG  +  R2n-VBRG 
dRG  (Rg  +  R2)2 


and 


dEo  =  VB 


dR 


(Rg  +  r2)' 


G  * 


The  gage  factor  is 
c  _  d  R/R 
dL/L 


then  d  Rq  =  F  Rq  dU/L 


and 


R« 

dE  =  VD  - — —  F  d  L/L 

O  B  r. 


(Rq  +R2) 
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where  Rj  -  Rq  . 


Therefore 


R?  Rg 


(Rg»-R2)2  (2R2)2 


then  dE  =  Vc  1  F  d  L/'L 
o  B  4 

and  d  L/L  -  f  =  strain. 

VB 

Consequently,  dE  -  F  t 

°  4 


Example: 

Foj 

100  \i  in/in  strain 

V.,  ~  1 1  volts 

F  =  2 

dE 

=  "<2>x  looxi.r6 

o 

4 

dE 

=  —  X  100  X  10"6 

o 

4 

dE 

=  550  x  !0~6  =  .550  millivolts 

for  100  n  in/in  stiain 


I .  Strain  Gage  Cal  i bration 


The  strain  gages  were  calibrated  using  the  shunt  calibration  method.  This  was  accomplished 
by  pressing  a  switch  on  the  Signai  Conditioning  equipment.  This  switched  a  resistor  in 
parallel  with  the  strain  gage.  A  simplified  diagram  of  the  circuit  is  shown  in  Figure  Al-ll  . 
The  parallel  combination  of  the  strain  gage,  R.-.  and  the  calibration  resistor  R  ,  caused  a 

AR  change  in  the  resistance  of  the  active  arm  of  the  bridge  circuit.  Four  calibration  re¬ 
sistors  were  chosen  so  that  the  AR  change  in  resistance  of  the  parallel  combination  would 
be  equivalent  to  500,  1000,  1500  &  2000  micro  in/in  of  strain.  When  the  strain  gages 
were  being  calibrated  the  output  of  the  signal  conditioning  equipment  was  o  d-c  voltage. 

To  facilitate  the  data  analysis,  this  d-c  signal  was  changed  to  an  a-c  signal.  This  was  ac¬ 
complished  by  monitoring  the  d-c  voltage  on  an  oscilloscope  and  then  using  an  oscillator 
to  provide  an  a-c  signal  with  a  peak  to  peak  value  equal  to  the  d-c  voltage.  To  compen¬ 
sate  for  possible  gain  or  attenuation  due  to  analyzing  equipment  rhe  a-c  calibration  signal 
was  sent  through  the  same  equipment  as  the  strain  gage  signal.  The  calibration  signal  pro¬ 
vided  a  deflection  of  the  x-y  recorder  pen  which  corresponded  to  units  of  strain.  Each  plot 
of  panel  response  was  preceded  by  a  calibration  which  gave  a  very  accurate  measure  of  the 
panel  response  in  terms  of  strain. 
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An  example  of  the  method  used  for  determining  the  values  of  the  calibration  resistor  R 
Figure  AI-11,  is  given  below. 

Example: 

Find  Rcq|  for  an  equivalent  strain  of  1000  ft  in/in 


AR/V 

F  = - 5  & 

AL/L 

AR./R  =  F  AL/L  =  Ft 


For  =  2000  x  10  6  where  R^  —  120  ohm 

rg  G 


AR  =  c  Rg  ~  ( 2000  x  10"°)  x  120  -  .2-100  ohm 


AR  =  R„  -  !S_!i2L 

RG+Rc.ai 


RG  Rcol 
RG  +  Rcal 


-  Rq  -  .  240  ohm 


120  R 


cal 


=  120  -  .24 


120  HR 


cal 


120  R  ,  «=  (120)2  +  120  R  ,  -  .24  R  . 
cal  '  cal  cal 


r  .  020/ 

“• '  _2T 


<=«60, 000  or  60K  ohms 


J .  Test  Equipment  List 


Microphones: 

Microphone  power  supply: 
Microphone  calibrator: 

Signal  conditioning: 

Data  amplifier: 

X-Y  recorder: 

Sweep  oscillator: 

Power  amplifier: 

Transducer: 

Random  signul  generator: 

True  R.M.S.  voltmeter: 

Spectrum  Shaper: 

Tape  recorder: 

Time  code  generator: 

Log  convertor: 

Power  amplifier: 

Speaker: 

Counter: 

DC  power  supply: 

Multi-channel  averaging  control: 
Accelerometers: 

Accelerometer  amplifier: 

Strain  Gages: 

Level  recorder: 

Oscillator: 

Narrowband  analyzer: 

Octave  band  analyzer: 
Oscilloscope: 

Oscillator: 

Analyzer: 

Multiplier: 


Bruel  &  Kjaer  type  4136 

Bruel  &  Kjaer  type  2801 

Photocon  model  PC- 125 

Endevco  model  4400. 10 

Hewlett  Packard  model  8875A 

Moseley  model  2FRA 

Spectral  Dynamics  model  SD- 104-2 

Ling  model  RP-3/5 

Ling  model  EPT  200 

Bruel  &  Kjaer  type  1 024C 

Bruel  &  Kjaer  type  2416 

Bruel  &  Kjaer  model  1612S/2 

Ampex  FR-1300 

Ampex  TCS-100 

Technical  Products  model  TP-662 
McIntosh  model  Me  240 
Altec  model  605B 
CMC  model  800A 
Kron-Hite  model  RS-3610SR 
Spectral  Dynamics  model  5D-22 
Endevco  model  2223 
Endevco  271 1 A 

BLH  type  DLB-A1 2-4AS1 3  &  Den- 
tronics  type  202C13 

Bruel  &  Kjaer  type  2305CS 

Hewlett  Packard  model  202C 

Spectral  Dynamics  model  SD-101A 

Bruel  &  Kjaer  type  2112C 

Tektronix  type  502A 

Technical  Products  TP-626T 

Technical  Products  TP-627T 

Technical  Products  TP-645 
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APPENDIX  II 


STATISTICAL  ANALYSIS  OF  DATA 

Frequently  it  is  necessary  to  determine  whether  an  apparent  relation  between  two  variables 
is  significant  and  having  shown  it  to  be  significant,  to  determine  the  best  form  of  represen¬ 
tation. 

When  pairs  of  observations  on  the  two  variables  plotted  on  rectangular  coordinate  or  func¬ 
tional  curve  paper  give  something  that  looks  approximately  like  a  straight  line,  there  are 
some  reasons  to  suspect  that  the  relationship  could  have  the  forms: 

y  -  a  +  bx 

or  if  transformed 

Y  =  AXB  . 

In  the  fitting  of  the  best  straight  line  to  data,  the  definition  of  the  best  line  is  that  line 
which  makes  a  minimum  the  sum  of  squares  of  the  deviations  from  the  line  of  the  measure¬ 
ments  of  y;  the  independent  variable  x  is  assumed  to  be  free  of  errors.  When  this  process 
is  applied  to  a  pair  of  variables  transformed  in  the  manner  indicated  above  what  is  being 
minimized  is  not  the  deviations  of  Y  but  the  deviations  of  log  Y. 

To  test  for  the  significance  of  an  apparently  linear  relation  a  correlation  coefficient  r  de¬ 
fined  by 

r  =  -&--!> _ _ 

2(x-x)2  2(  y-9)2 

where 

s(x-s)2=  s 

n 

Z(y-y)2  =  Sy2-iM.2 

n 

2  (x  -  x)  (y  -  y)  =  2  x  y  -  —  x  -  S  Y~ 


n  =  number  of  pairs  of  observations 

The  correlation  coefficient  has  the  characteristics  such  that  if  the  relationship  between 
the  data  can  be  represented  exactly  by  a  straight  line  then  r  =  ±  1,  positive  if  the  straight 
line  has  positive  slope  and  negative  if  the  line  has  negative  slope.  If  on  the  other  hand 
there  is  no  relation  at  all  between  the  variables  than  r  is  zero. 
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If  r  proves  to  be  significantly  greater  than  zero  the  matter  can  be  taken  farther.  The 
values  of  a  and  b  can  be  evaluated  from 


b  =  2  frc  -  x)  (y  -  y) 

2(x  -  x)2 

The  regression  line  wili  have  an  equation  of  the  type 
y  =  a  +  bx 


etc. 

It  is  assumed  that  the  regression  line  is  established  with  points  scattered  about  on  either 
side  of  it.  The  standard  deviation  of  this  scatter  E,  measured  in  units  of  y  parallel  to 
the  y  axis,  is  given  by 


■  ■  ^  • 

Suppose  it  is  desired  to  draw  95%  confidence  limits  on  either  side  of  the  regression  line, 
within  which  95%  of  all  points  should  lie.  First,  the  value  of  t  for  n  2  degrees  of  free¬ 
dom  and  the  5%  level  of  significance  is  found  from  a  Student's  t  table.  The  product  of 
tE  is  determined  next.  Last,  two  lines  are  drawn  parallel  to  the  regression  line,  but  one 
displaced  tE  units  of  y  downwards  and  one  displaced  tE  units  of  y  upwards.  Thus,  if  the 
regression  line  is  used  to  predict  the  value  of  y  from  a  known  value  of  x,  these  confidence 
limits  give  the  approximate  limits  between  which  there  is  a  95%  chance  of  obtaining  a 
correct  prediction. 
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APPENDIX  III 


BRIEF  DERIVATION  OF  THE  PLATE  ELEMENT  STIFFNESS  MATRIX 
Consider  the  thin  plate  element  shown  below 


z 


Assuming  a  deformation  of  the  form 
w,(x)  =  X(x)  along  lines  of  z  =  constant 

i 

w2(z)  =  Z(z)  along  lines  of  x  -  constant 

3  2 

where  X(x)  -  A^x  +  A^x  +A^x+A 


0) 


Z  (z  )  =  ®3Z^  +  +  BjZ  +  Bq  j 

the  coefficients  A.  and  B.  w ;  1 1  be  determined  from  the  boundary  conditions  at  the  corner  of 
the  plate  element. 

Along  the  edge  z  =  0 

X(o)=Ao  =  S3 

X  (b  )  =  Agb^  +  A2b^  +  Ajb  +  Aq  =  5^ 

X’(o)=  Aj  =  ®z3  '  (2) 

X'(b)=  3A3b2  +  2A2b  +  Aj  =  -e  4 

4 

where  §.  is  the  displacement  at  rhe  i^  comer 

9^ .  is  the  rotation  about  the  z  axis  at  the  i^  comer 

0^.  is  the  rotation  about  the  x  axis  at  the  i  comer 
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Solving  Equations  (2)  for  the  coefficients,  A.  , 

A  =  S 
o  3 


A1  "  °z3 


2e  +Ie  - s  +.1  a, 

^  z3  ^  z4  2  ^  2  4 

^  b  b 


A,,  =  J-0  ,  -  JL  0  +  i_  a  -  JL  fi 

w  ft  z3  a  z4  #■»  3  ^  4 

b2  b2  b,J  bd 


(3) 


Then,  along  the  edge  z  -  0 


2 x3  3x2 


3,2  3  2 

X  (x)  «  (^~^  +  x)6  ~-)9  4  +  ( 

o  2  b  Zw  2  b  24  3  2 

b  b  b  b 


+  i)a. 


^  2x3  3 x2 ^  q 

,3  .2  4 

b  b 


(4) 


Similarly,  along  the  edge  z  =  a 


Xa(x)  =( 


3,2 

!L-^L-  +  x ) 9  . 

_  I  Zl 

b2  b 


-( 


3  2  3  2 

iL.  -  )9z?  +  (~  -  +  1)5 


2  b 


-  5iL  )  5 

b3  b2  ' 


(5) 


along  the  edge  x  =  0 


Z  (z) 
o 


3,2  32  ,3,2 

(—  +z)9 —  )0x1  +(—  +  1)  S, 

2  a  2  a  xl  3  2 

a  a  a  a 


,3  ,2 

(iz_  -  )  § 

3  2 

a  a 


(6) 
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and  along  the  edge  x  =  b 


3  0  2  32  ,  3  ,2 

zbM  .  <i-  -S-)Qx2+(2i-  -Si.  +  i)*4 

2  a  2  a  3  2 

a  a  a  a 


-  (id -2i1)s7 

3  7  1 

a  a 

For  thin  plates  the  strain  energy  expression  is 


(7) 


U 


°  I  f 


(w,  +  w,  r-2(l-i')  (w,  w, 

XX  zz  '  XX  'zz 


-  w,  )|dxdz  (8) 
xz 


Letting  w(x,z)  =  Wj(x)  +  w2(z) 

the  assumed  deformation  pattern  will  be  taken  as 


w(x,z)-Xo(x)  +  Zq(z) 
w(x,  z)  =  XQ(x)  +  Zb(z) 
w(x,z)  =  Xq(x)  +  Zq(z) 
w(x,z)  =  Xq(x)+  Zb(z) 


0  <  x  <  b/2,  0  <  z  <  a/2 

0  <  x  <  b/2,  a/2  5;  z  ^  a 

b/2  <  x  <  b,  0  <  z  <  a/2 

b/2  <  x  <  b,  a/2  5-  z  <  a 


(9) 


(10) 


Equations  (10)  are  used  to  evaluate  the  bending  contributions  to  the  strain  energy  through  the 
terms  w,  and  w ,  .  The  contribution  to  the  strain  energy  for  the  torsion  of  the  plate  is 

obtained  oy  considering  the  curvature  to  be  defined  as 


w'xz  =  *]  +  S2+  S3  ”  S4)/ah  HD 

Then,  substituting  Equations  (10)  and  Equation  (11)  into  the  strain  energy  expression  (Equa¬ 
tion  (8))  and  performing  the  integrations,  the  strain  energy  is  obtained  in  terms  of  the  corner 
displacements  <5  .,  0  .,  and  0  ..  The  stiffness  matrix  is  obtained  by  application  of  Casti- 
gliano's  Theorem. 


m 


From  the  definitions  given  in  Equation  (  58  )  the  plate  stiffness  matrix  is  ( see  Equation  (  43 ) ) 
for  the  i  plate  element. 


K 


xx 


2a. 

0  a . 

i 

0 

0 

2a.  0 

i 

a . 

i 

(12) 

a , 

i 

0  2a. 

i 

0 

0 

a.  0 

i 

2a. 

i 

-P. 

30(5/fc> ). 

a. 

i 

-6(8  A ); 

30(S/b ). 

-  P. 

i 

-6(8/b ).  ff. 

i  i 

03) 

-  o. 

6(S/b). 

P. 

i 

-30(S/b  )j 

6(8/b). 

-a. 

i 

-30(  8/b ).  P. 

2*i 

*i 

0 

0 

a. 

i 

2^i 

0 

0 

(14) 

0 

0 

23, 

0 

_ 

0 

0  i 

23;  _ 
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K 

yy 


r 


X. 

1 

-  T  . 

1 

-30  (5/a). 

6(5/a). 

T  j 

-  X. 

-6  (5/a ). 

30  (5 /a). 

-30  (5  /a). 

6  (  S ,/a  ). 

X; 

-  x . 

i 

-6  (S/a) 

30(5 /a ). 

t  . 

1 

-X. 

i 

(1 

5) 

K. 

1 

-X; 

-  n. 

V .  +  72  ^5/ab ). 

-X. 

1 

K. 

1 

Y.  +72  (5/ab). 

-^i 

-"i 

y.  +  72(  5/ab ); 

K  m 

i 

-X, 

Y.  +  72(5/ab). 

-M. 

K, 

i 

(16) 

where  only  the  upper  diagonal  sub  matrices  have  been  shown. 
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APPENDIX  IV 
COMPUTER  PROGRAMS 


A.  Substructure  Dynomic  Analysis 

The  computer  program  givers  below  uses  the  equations  developed  in  Section  1 1  — D  and  the 
method  of  solution  given  in  Section  ll-D-3. 

The  only  option  is  whether  or  not  eigenvectors  are  desired.  For  print  out  of  eigenvectors 
set  MV  =  0.  The  input  data  are  defined  as 

E(l)  Young's  Modulus  for  the  i^’  stiffener 

J 1  ( I)  St.  Venant's  torsion  constant  for  the  i^  stiffener 

l!(l)  second  area  moment  of  Inertia  for  the  i^  stiffener 

L 1  ( I)  length  of  the  i"1  stiffener 

EP(I)  Young's  Modulus  for  the  i^h  panel  element 

HP(I)  thickness  of  the  i™  panel  element 

AP(I)  length  of  the  i*^  panel  element  in  the  z  direction 

BP(I)  length  of  the  i^  panel  element  in  the  x  direction 

WP  weight  of  panel  element  per  unit  area 

WR  weight  of  rib  per  unit  length 

WS  weight  of  stringer  per  unit  length 

All  dimensions  are  in  inches  and  weight  in  pounds. 

B.  Effect  of  Curvature  on  Stress  in  Honeycomb  Panels 

The  program  uses  the  equations,  summarized  in  Section  V-2  in  the  calculation  of 
stress  ratio . 

There  are  two  input  options.  OPT1  selects  the  type  cf  excitation;  0  for  random 
noise,  1  for  sinusoidal  excitation.  OPT2  allows  the  inclusion  of  variations  of  damping  due 
to  change  in  natural  frequency  with  curvature;  0  leaves  out  damping,  1  includes  variable 
damping . 

There  are  three  lines  of  input  data.  The  first  two  lines  contain  the  physical  data.  The  data 
are 

A  arc  Jength 

B  panel  length 

E  modulus  of  elasticity  of  skin 

GXZ  Transverse  shear  modulus  for  core 
GYZ  Transverse  shear  modulus  for  core 
H  core  thickness 

R  panel  radius  of  curvature 

T  skin  thickness 

TE  edge  thickness 

RO  skin  mass  density 

DELC  damping  ratio  for  first  mode  of  curved  panel 

DELF  damping  ratio  for  first  mode  of  flat  panel 


The  data  are  prepared  as  follows; 
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c  THIS  PRC CHAM  COMPUTES  NATURAL  FREQUENCIES  AND  MODE  SHAPES 

C  OF  A  NINE  BA Y  0KTH060NALL Y  STIFFENED  PANEL  ARRAY  WITH  Cl.AMHED 

c  EUGtS.  SUBROUTINES  REuUIheD  M<E  LOC#EI"gEN  # AND  NROOT. 

REAL  E*  1<*>  »U1(  12)  »  U  (12)  »L1  (12>  *B(12>  #T(  J2)  »Kl  (12#12> 

REAL  El’  ( 9 )  ♦  HP  ( 9  >  » AP  ( 9)  # dP  < 9  >  » L  ( 9  )  »  AL  (  9  >  *  BE  (9 )  » 6M  (  q ) » DE ( 9 ) 

REAL  Pi<9)  »P2(9)  *P3«9)  »P4(9)  #F’b(9)  »P6(9lrJ57(9) 

REAL  Kp(12#12) #K2(12*12)  *  MM (12*12) 

_ DOUBLE  PRECISI ON  EVL(12) #EVC ( 144) # AA( 144) #Bd( 144) _ 

lFnZG£p~oTs 

C  IF  eigenvalues  and  eigenvectors  are  desired#  set  MV=0  . 

C  . IF  ONLY  EIGENVALUES  ARE  DESIRED#  5FTW=i  ". .  ' 

C  SEE  COMMENTS  IN  SUBROUTINE  •  EIGEN*. 

I  t  1  REAL)  ( S ,  1 )  MV 

RE AD ( b »  2 )  WP»wR*WS 
Nr  12 

Sr  9  __ 

READ  t  5 »  3 )  (E( 1) » 11 (I ) » Jl(I) #Ll (I ) *  I  =  1#N) 

RE  AU  ( b »  3 )  (EP(  I )  #HP(  I )  .API  I  )  ,£iP<  I)  »  I=l»S> 

1  FORMAT (313) 

_ 2  FORMAT  (3El2.b) _ _ 

3  FORMAT <4E12.b) 

WRITE! o*4)  __  _  

4  FORMAT (1H1#BX. • VIBRATION  ANALYSIS  OF  A  STIFFENED  PANEL*)' 

C  COMPUTh.  STIFFNESS  OF  SUPPORTS  G  STRUCTURE 

DO  S  Orl »  N 

_ I1(0J £__i  (u)*U(Q)/Ll(Q)**3 _ 

T(u)r  u . 384 ♦t (Q)*J1(Q)/L1(Q) 
b  CONTINUE  __ 

K 1 ( 1 . 1 ) r  T ( 2 ) +  T ( 3 ) +  4 . 0  * ( d ( 4 ) *  L 1 ( 4 ) *  *  2 + B ( 1 2 ) ♦ L 1(12 ) *  *  2 ) 

Kl(lf2)=-T(3)  _  _  __ 

M ( 1 # 3 ) r  2.0*B(12)*Ll(12)**2 

K 1  ( l .  9  )=  6.  ujdi  ( 12)*Ll(12)-d(4)*Ll (4)  ) _ 

K 1  ( 1 » 1 1  j  r-fa .  0*ii  (12  )  *L1  ( 12  ) 

K1  ( 2  #  2  (  -  T(3)+T(7)44.U*(b(5)*Ll(b)**2+B(6ULl_(J2.).*t2.) _ 

K1 (2#4 ) r  2.0*B(t>)  *Ll  (fo)**2 

Klt2#l  i)r  6.0*(d(6)*Ll(b)-0(S)*Ll(5> ) 

K1  (2#  1  )=-6.0*h<6)  *H  (6> 

Kl(3#3)=  T(l)+T(9)+4.u*(G(li)*Ll(ll)**2+B(12)*Ll(i2)**2) 

K 1  IT-  4  )  r'-T  (  V) . . 

Kl(3#9)=  6 • 0*B ( 12 ) *L1 ( 12 ) 

K1(3»U)=  S.  0*  ( d  ( 1 1 )  *L1  ( 11 )  -p  ( 12  )  *El  ( 12  )  ) 

M (4 #4 ) r  T(a)  +  r(9)+4.U*(B(6)+Li(b)**2+B(10)*Ll (10)**2) 

K 1  ( 4  #  1  c  )  r  t> .  0  *ii  ( 6  )  *L1  ( b  ) 

K 1  ( 4 » 1  )  r  6.0*  (u(  10  )  *H  (  10  ) -t '  (b)  *L1  (6)  ) 

K 1 <  b#  b ) r  f ( 4) +  T (12 ) +4 . 0 * ( B ( 2 ) * L 1 ( 2 ) *  *  2 +G ( 3 ) *L 1 ( 3 ) *  *  2 ) 

K 1  (  5 » t> )  r  2.0*B(3)*L1(3)**2 
K 1  (  b »  7  ) =—T ( 12) 

K 1 ( b » 9 ) r  b.0*(B(3)*Ll(3)-B(2)*Ll(2)) 

Kl(b»lu)r-b.U*b(3)*Ll (3) 
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Kl  (Of  6)=  T  i  b>  +1  (6  5  +4. 0*<H(  3)  *U  (  3)  **2+H(  y>  *L1  (  7)  *  +  ?.) 

. . . K  l  16,  t V)  =  -T  (  6) . . . .  .  “ .  . 

Kl<o»9)=  b.0*B(3)+Ll  (3) 

K 1  ( o  f  i  u  )  =  6.0*(13(7)*L1(7)-H(31*L1<3)) 

K 1  <  7,1)-  T<ll)  +  TU2)*4.U*(ti<  I  )*L1(1)**2+B<9)*L1(9)**2) 

Kl.(7.b)  =  ii » 0  *  ii  ( 9 )  *  L 1  ( 9 )  *  * 

K 1 (  / » 1 1 ) -  b. 0*  <  b<  9) *L1 ( 9  > -H ( 1 ) *L1 ( 1 )  > 

- XTT77 17 Ts=5nr*ifnjT iot9T‘  . * - - - 

Kl <b;8)=  T (o) +T ( 10) +4.0* (B(b) *L1 (6)**2+B(9) *L 1 ( 9 ) * *2. ) 
Kl(a*U)=  6. U*b ( 9 ) *L 1 <  9 ) 

Kl  (b*  1^)=  b,O*(U(0)*LHa)-H(O)*LK9)) 

'  k II 9 7 9 ) ~  1 2.0 ♦Id ('27  +G73 r+oTu )  +  d (lTH 
Kl  (9»lu)  --12.0*13(3) 

Kl(9r  U)=-I2.0*tf(12> 

Kl(lo»iO)=  12.O<MB(3>+B(b)+B(o)+0<7>  ) 

KlUU»i2)=-U'.U*B(b) 

_ Kl < 11.  1 1 )  =  1 2_.  0  *  ( 6_(  1  )  +  b  <  9 )  til  ( 1 1  >  +  0(12)  ) _ _ 

K  iTU*'  i  2T=-l2  ,'0*6'(  9  > 

Kl(12»i2)=  12 • 0* ( B ( 6 ) +B ( b) *b ( 9) +B ( 10 )  ) 

C  COMPUTl  stiffness  of  COVER "plates  "  . . 

Du  o  Q^1»S 

u  (  Q )  =  U  .  0  9 1 2  *EP  ( 9 )  *  HP  ( Q )  **3 
AL(Q)=U(Q)*bP(Q)  /AP(Q)**2 
BE  ( W )  =  )  ( a )  *AP  (  QTToT3!  Q )  **2 

GM<Q>=1.4*D(Q)/(AP(0)*BP(G) )  _  _ _ _ 

Dt ( 0 ) r  0.01875*10(0) 

Pl(u)  =  3.0*(AL(O?-H0.0*UE(Q)/dP(Q)  )  _ 

P2(0)=  3 . 0* ( BE (Q ) +10*  0*DE ( Q)7aPTQ) ) 

_ _P3^0 )_•!_ 3 .  u*  (  AL  (JJ )  +2 , 0 *DL(Q)/RP(Q)  ) _ _ 

P4  ( O )  =  3".  0*  <  BE  ( Ql  +2. 0*DE  ( Q )  /AP  ( Q )  ) 

P5(0)=  GM(O)+72.0*UE(O)/(AP (0) *BP(Q) ) 

P6(0)=  b • 0*AL ( Q) ZAP ( G ) 

P7(u)=  b.0*BE(O)/bP(Q) 
b  CONTINUE 

_ K  P  ( 1 » 1 )  -  2  «  0  * ( AP ( 5)»AL(5)+AP(4)*AL(4)+AP(2>*AL(2)+AP(1)*AL(1) ) 

KP (1.3)-  AP ( 5) *AL ( 5) +AP ( 4 ) *AL ( 4 ) 

KP<1.9)=  PKD+Pl  (2)“P1(4)-Pn5) _ _ _ _ 

KP(l*lo)=  3U.b*(DE(5)/bP(5)-De(2)/BP(2) ) 

K  (-'(l.ll)—  P3(4)+P3(5) 

KP  ( 1 . 1,'i )  ~-6. 0*Dt  (  5)  /BP ( b) 

_ K  P ( g.2)  =  2 «0*( AP (6 )»AL(b)»AP(5)*AL(5)+AP(3>*AL(3)+AP(2)*AL(2) ) 

K  p  f  2  (4)-  AP(5)  *A"L(S)+AP(b)  *AL(6> 

KP (2.9)-  30.0*(D£(b)/OP(5)-DE(2)/BP(2) ) 

KP  (2.  if, )  =  PI  (2  >  +P1  (  3 )  -PI  (  5)  -Pi  (6 ) 

Kp ( 2 . 1 1 ) =-6. Q*DE ( 5) /OP  (5) 

KP(2»1,  )=  P3(5)+P3(6) 

KP(3.3)=  2  -  0* ( AP ( fi) *AL ( b)  4  AP  (  7 )  *AL  (  7 )  +AP  (  5 )  +AL  ( 5 )  ■♦  AP  (  4 )  *AL  ( 4 )  ) 
K  fT3 . 9 )  -  -F3T4  r^PT(T)  ~ 
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KP(3.lO>=  6,0*DE(S>)/BP(5) 

KP(3»li  )  =  Pi  (4 )  +P1( 5 ) -Pi  (7) -Pi  ( ST 
KP(3#l2>=  30  •  0* ( DE (8) /BP ( 8 ) -DE ( 5 ) /BP ( 5 ) ) 

KP (4 » 4 )=2. O  ’*  ( A  P (9 ) *AL (9 ) +AP (ft ) *AL (8 >  + AP ( 6 ) * AL (6) +AP  ( 5 ) *AL ( b  )T 
KP(4.9)=  6.0*DE(5)/BP(5) 

KP<4»lui=-P3<5)-P3<6) 

KP(4»li)=  3O.O+(DE(0)/BP(8>-Dfc(5)/6P<5) ) 

KPT4 .  ly )  =  Pi(5i +PKbl -Pi ( 8) -Pi (9T 

KP(b»5)“  2.0*(BP(l>%BE(i>+BP(2>*B£<2)+BP(4)*BE<4>+RP<5)+BE(5>  > 
KP  < b  .  6 )  =  BP  (5 )  V8E  ( 5 )  +BP'(  2 )  *bE  ( 2 ) 

_ KP(b.9)=  P2  ( 2 ) +Pg ( b ) -P2 ( 1 ) -Pg ( 4 ) _ 

KP ( 5 ► 1 u  > =-P4 ( 2 ) “P4 ( 5  J 

KP(5#li)=  30.0*(DE(4)/AP(4)-DE(5)/AP(5) ) 

KP  fb  Via) =  b . 0*0£( 5  >7APT5T 

KP ( 6 » b )  =  2 . 0  * ( BP  <  2 ) ♦BE ( 2 ) +BP ( 3 ) *BE ( 3 ) +BP ( 5 ) *BE ( 5 ) +BP ( 6  > *BE « 6 ) ) 
KP (o » 9)=  P4  (2 )  +P4  <  5 ) 

Kp«b»l0>=  P2(3)+P2(6)-P2(2)-P2(5) 

KP ( to* TI i =-6 . 0 ♦DE ( 5T7 AP ( b ) 

KP( bn  12)=  30.0 ♦ i DE ( 5 ) / AP ( 5 ) -DE (6 ) /AP (6)  ) _ 

KP(7r7)=  2.0*(BP(4)»BE(4)+BP(b)*BE(5)+BP(7)*BE(7)4BP(8)*BE(8) ) 

_ KP  (7.8)=  SP  ( 8 )  *BE  ( 8 )  +BP  ( 5 )  *BE  (  5 ) _ _ 

KP(7.9)=  30 .0* (DE (4 ) /AP (4 ) -DE ( 5) /AP ( 5 ) ) 

KP(7.1'j)  =  6. 0*D£ ( 5) /AP ( b ) _ 

KP(7»U>  =  P2(5)4-P2(8)-P2(4)-P2(7> 

_ KP(7.1P)=-P4(5)-P4(8) _ _ 

KP(8»8)=  2.0+(BP(5)*BE(b)+BP("b>*BE(6)+BP(8)+BE(8)+BP(9)*BE(9)  ) 

_ KP (8.9 ) =-6«  0*DE ( 5) /AP (5) _ 

KP(8.1u>=  30.0+(DE(5>/AP(5)-DE(6)/AP(6) ) 

KP  ( 8. 1 1 )  =  P4(5)4~P4(8) _ 

KP(8.12)=  P2(6)+P2(9)-P2(5)-P2(8f 

_ KP(9.9)=  P5(  1 ) +P6 ( 1 )  4P7  ( i  )  +P5 (2)  +P6(2)  +P7 (2 ) _ 

i  +P5  ( 4 )  +P6  ( 4 )  4-P7  ( 4 )  +P5  ( 5 )  +P6  ( 5 )  +P7  ( 5 ) 

_ KP(9.iui=-P5(2)-P7(2>-P5(5>-P7(5) _ 

KP (9* li ) =-P5  4)-P6(4)-Pb(5)-P6(5) 

KP(9.iP)=  P5CS) _ 

KP(10» 10)=  P5(2) +P6 (2 ) +P7 (2) +P5 (3) +P6 ( 3) +P7 ( 3 ) 

_ J. _ 4-P5 (5)  +P6 ( b )  +P7  (5) +P5 (6 )  +P6 (6) +P7 (6) _ 

KP(10»il)=  P5(5) 

KP(1Q.12)=~P5(5)-Pb(5)-P5(6)-P6L6)  _ _ 

KPUl.  il  )=  Pb(4) +P6(4) +P7(4 )  +p5 (5)  +P6 ( 5)  +P7 (5) 

_ 1 _ ♦P5(7)+P6(7)+P7(7)+P5(8)+P6(8)+P7(8) _ 

KP (  i  1 .  i2 )  =-Pb  ( b)  *-P7 ( b)  -P5  ( S )  — P 7(8) 

_ _ KP (i2» 12 )=  P5(5)+Pb(5)+P7(5)+P5(6) +P6(6)+P7(6) _ 

1  ~  +Pb( 8) +Pb( 8) +P7 ( 8 ) +P5 (9) +P6 (9)+P7 (9) 

C  COMPUTE  STIFFNESS  OF  STIFFENED  PANEL  ARRAY _ _ 

DO  7  1=1 .N 
DO  7  J=1*N 
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K1(X*J)=  K1  ( I »  J)  +KP  ( I » J ) 

K 1  ( I » J )  _  _  _ _ 

K2(I»J)= 

7  CONTINUE _ 

c  compute  maSS  Matrix 

rtP=WP/^86.0  _  _ _____ _ _ ___ _ 

WR=WH/ 386.0 
WS=WS/38fo.O 

OPi=WP*(Li( 12) +L1 (4)  >*0.1(2) +U1 ( 3TT 
DRl^WR* (Ll ( 12) ♦Ll (4) ) 

DSl-WSt ( Ll (2)+Ll (3T) 
Pll=Ll(2)**2-Ll(2)*Ll(3)*Ll(3)»*2 
P21=  Li ( 12 ) ♦*2-Ll ( 12T *Ll ( 4 ) +Ll (4) **2 

P31=  Lt(3)-Ll(2)  _ _ _ _ 

P4i=  Li (4) “Ll ( 12) 

MM(1»1)=  P2l*(DPl+DRl)/3.0 

- MM  (  i  VS  j  =-P3 1  *'P4 IWRLTO - 

MM(l»9)=-P41*(DPl+DRl)/2.0 
'  MM ( 5 i 5 ) sPl 1 * ( DFi +DS lT 7 3.0 
MM(5#9)=P31*IOP1+DS>1)/2.0 
MM « 9 » 9 ) =  DPI ♦DR 1+DSl 
DP2=  Wp*(Ll(b)+Ll(5)  )*(Ll(3)+lU7)  ) 

- DR2="U;r*"TLX(  b  )"+LK  55  J - 

OS2=  WS*<L1(3)+L1(7) ) 

~Pl2=  L  i  (  31  **2-Cl  ( 3F*XlT7T «Xl  ( 

P22=  Ll(6)**2-L1(6)*L1(5)4L1(5)**2 
P'32=  L i ( 7 ) -Ll ( 3 ) 

P42=  L  i  ( 5 ) -Ll ( 6 ) 

MMl2 » 2  )  s'  P22*  <  bP2+t)R2 )  /3  ♦  0 
MM(2.6)=-P32*P42*DP2/4.0 
"AM (2 » 1 0 ) =-P42  * ( OP2+DR2 )/2 . 0 

_ MM(6>fe)=  P12* (0P2+DS2) /3»  0 _ 

MM(6«lo)=  P32*(DP2+DS2)/2.0 
MM(10»i0)=  DP2+DR2+DS2 
DP3=  Wp*(H(ll)+H(i2)  )*(L1  ( 1 )  +L1  (9)  ) 
DR 3—  WR*(L1(U)+LK12) ) 

DS3=  WS*(L1 (1)+Ll (9) ) 

Pl3=  Lj ( 1 ) ♦♦2-Lj ( 1 ) *Lj ( 9) ♦Ll ( 9 ) **2 
P23=  Li (11)**2-L1(11)*L1(12)+L1{12)**2 

P33=  L i ( 9)-Ll ( 1) _ 

P43=  Li ( 12) -Ll (11) 

_ _ MM(3»3)=  P23*(PP3+DR3)/3.0 _ 

MM(3»7)=-P33*P43*QP3/4,0 

MM(3»U  )=-P43»<DP3+DR3)/2.0 _ 

MM  ( 7  »  7 ) =P13* ( QP3+DS3 ) /3  *  0 

_ MM(7»ll)=  P33* (OP3+DS3) /2 • 0 

MM(il»il)r  DP3+DR3+DS3 
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DP4=  top* ( LI ( 10 ) +L1 ( 6 ) ) * (Ll (9) +  L1 (8)  ) 
DR4=  WR*(L1(10)+L1(6)  ) 

DS4  =  to‘.,*( L 1  ( 9)  »L1  (  a  )  ) 


P14=  Li(9)**2-Ll  (9)  ♦LKBM-Ll  (fl)**2 

P24=  Li (10) **2-Ll (10) *L1 (6)+Ll <6)**2 

P34=  Li (8)-Ll (9) 

P44=  Li(b)-Ll(10) 

MM  ( 4  »  4 j=  P24* ( UF4+DR4 }/3.0 

MM ( 4 »  8  > =-P34*P44*DP4/4 . 0 

MM(4»l,>)=-P44*(DP4+DR4)/2.0 

MM ( 8 »  B ) -  P14* ( 0P4+DS4 ) / J, 0 

MM(8#l;:)=  P34*  (DP4«-0S4)/2,Q 

MM‘12»i2)=  0P4+0R4+DS4 

DO  8  Irl « N 

DO  8  J=1»N 

MM(J»I)=  MM ( I » J) 

6 

CONTINUE 

DO  9  I  =  1»N 

DO  9  J=1»N 

CALL  LoC(i»U»lJ»N#0) 

AA(IJ>=  M ( I » J) 

9 

C 

OB(IJ)-  MM ( I # J ) 

CONTINUE 

COMPUTE  EIGENVALUES  AND  EIGENVECTORS 

CALL  NpOOT ( N#  AA . BQ »  EVL » EVC ) 

DO  13  I=1»N 

IF  (EVL  ( I )  )  10»H»11 

in 

ll 

£VL(I)=  DABS (EVL ( I ) ) 

EVL (I)r  0 * 1591*DSQRT (EVL ( I ) ) 

12 

WR I TE ( b » 14 )  EVL ( I ) 

IF(MV)  12 » 12’ 13 
toRlTE(,j#15) 

IX=  1  +  rj*  ( 1-1 ) 

IY  =  N+e* ( I— 1 ) 

13 

14 

WRlTE(o’ib)  (EVC { J) t  J=IX.IY) 

Continue .  " 

FORMAT (5X» «FREQUENCY=*»1X’E12,5»1X’ ’HZ.* ) 

15 

10 

FORMAT ( 5X ,  ’ MODE  SHAPE’) 

format { sx ’ 4e 12 . b ) 

D0  17  J=X7l2 
DO  17  u=l»12 
Kl (1 t J)=  0.0 

_ K2 ( I t J) s  0 .  u _ 

kpTTTJj  =  o.o 

CALL  LyC (I»JfIJ»l2#o) 
AAUJT;  0.0 

_  Bb( IJ) -  0.0 

17  CONTINUE 

00  TO  HI _ 
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■  noon  do  n|n  no  dn  r.;n  n  nolo  n  n  nn  n  n  n  n  Pio  r  o  n  n  dr.  f>rnr,  rir, 


C 


“Sub»T6uTTNE"EI gEn  ' 


PUKi  JOSE 


COMPUTE 

!ATKlX 


EIGENVALUES  AND  EIGENVECTORS 


OF  A  REAL  SYMMETRIC 


usa,;t 

CALL  EIGLN(A»R#N»MV) 


description  of  parameters 

a  -  ORIGINAL  MATRIX  ( SYMMETRIC  T* . DESTROYED . I N~C OMPUTATION. 

RESULTANT  EIGENVALUES  ARE  DEVELOPED  IN  DIAGONAL  OF 

'matrix ‘a’  IN'  descending  Oiicror. 

-  RESULTANT  MATRIX  OF  EIGENVECTORS  (STORED  COLUMNWISE# 

IN  SAME  SEQUENCE  AS  Ei  GEN  VALUES)  . ' .  . “ 

,  -  ORDER  OF  MATRICES  A  AND  R 

iV~  INPUT  CODE.'  . . . . . . .  . . .  . . . 

0  COMPUTE  EIGENVALUES  AND  EIGENVECTORS 

' .  I  CbMPuTcTTr I GE N V ACuITsT- ONLY  (R  NEED  NOT  BE 

DIMENSIONED  BUT  MUST  STILL  APPEAR  IN  CALLING 
SEQUENCE )  - . .  .  . 


REMARKS  . .  . .  '  .  - 

original  matrix  a  must  be  real  symmetric  (storage  mode=i> 
IaTrIx  7T  CaKINOTBEIN  ThE  SAME  LOCATION  AS  MATRIX  R - 

SUBROUTINES  and  funC T  ION  SUBPROGRAMS  ~  RE QU TreET  . 

NONE 


method 

■“  (TAGONALIzatION  method"  ORIGINATED  by  JACOBI  AND  ADAPTS - 

BY  VON  NEUMANN  FOR  LARGE  COMPUTERS  AS  FOUND  IN  *  MATHEMATICAL. 

methods  for  digital  Computers* »  edited  by  a,  ralSton  and 

H*S.  WlLF#  JOHN  w I LEY  AND  SONS#  NEW  YORK#  1962#  CHAPTER  7 


P 


SUBROUTINE  EIGEN(A»R#N»MV) 
D I MENS  I ON  A ( 1  f » R ( 1 ) 


© 


IF  A  DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED#  THE" 

_  _ G_Ih_COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PRECISION 

C  STATEMENT  WHICH  FOLLOWS.'  '• 
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O  OjO  j  1  !  i  ojnoj  J  l  j  jo  ojo  r>  o  ojo  o|r>  < 


DOUBLE  PRECISION  A ♦ R * ANORM * ANRMX ► THR * X » Y * SINX *5lNx2 * 

_ 1  _ _ _ CQSX2*SINCS*C0SX _ 

c 

C  THE  C  MUST  ALSO  BE  REMOVED  FROM  DOUBLE  PRECISION  STATEMENTS 

Z  APFe'aRTRG-!^  OTHEFTT^TJOTTnes  used  IN  C ON JUNCTTSN~wTtTr THIS  "' . . . 

C _ ROUTINE. _ _ 

r  ’  ~  ~ 

THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 
CONTAIN  DOUBLE  PRECISION  FORTRAN  FUNCTIONS.  SORT  IN  STATEMENTS 

40*  66*  75*  AND  78  MUST . bE  CHANGED  TO  DSQRT .  ABS  IN  STATEMENT 

62  MUST  BE  CHANGED  TO  DABS. 


GENERATE  IDENTITY  MATRIX 


IF(MV-li  10*25*10 

10  I QF-N _ 

DO  20  J=1*N 

Ift=IQ+N 

DO  20  I=1»N 

IJ=I(S+1 

R ( I J) -0  •  0 

IFU-U)  20*15*20 

15  rTTJT=T7o  “ 

20  C0NTINuE 


COMPUTE  INITIAL  AND  FINAL  NORMS  (ANORM  AND  ANORMX) 


25  ANORM=CJ.O 

- - 

DO  35  J=I*N 
TF(T-J)  30*35*30 
30  I A~I+ ( J )  /?. 

ANOfTMiANORM+A TI A )  *A<  TAT 
35  CONTINUE' 

- 1 F  ( ANCTRWTT55V  T65V¥{J - 

40  AN0RM=l.414*DSQRT< ANORM) 

AN'fiMX= aNOrA*  1 .0Y-67FUOAT(  N ) 

I^I^-iAl1Z£  indicators  and  compute  threshold*  Thr 
IND=I) - 

THR=AN0RM 

45  T H RL T  HrTFl  OATCFTJ 
50  L=1 
55  M=L+1 

Z  COMpuTE  SIn  and  CoS 

r 

'W 
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60  MQ=<M*m-M)/2 

lQ=(L*L-U/2  _ _ 

LM=L+Mq  ' 

62  I FJOABSJ  A  (LM))~THR)  1 30*65  *65 _ 

""  65  lND=i 

_LL=LH-0 _ _ 

MM— M+Mq 

X=Q.5»{A(U-)~A(MM)  ) _ 

6fl  Y=~A<LM)/DSaRT<AU.M)*A(LM)+X*X> 

IF ( X )  70*75 * 75 _ 

7  b  "r =-r~ 

75  5INX~ Y/DSQRT (2. 0* ( 1 . Q+( DSQRT ( i »0-Y»Y) ) )  ) _ 

SINX2=SINX*SINX 

78  COSX=DsQRTU«0~SINX2) _ 

COSX2~COSX*COSX 

_ SINCS  rSlNX*COSX _ 

_  ROTATE  l  and  m  columns _ _ 

ILQ=N»(U-l) _ 

IMQ=N*(M-1) 

_ DC__125  _I=1#N _  _ _ 

Tq=(i*i-i>/2 

_ I FJ  I-Li  80  *1X5 *60  _ 

80  IF ( 2~M)  85* 115*90 
85  IMsI+Mcj 
IQ  TO  95 
90  IMsM+Iq 

95  IF  ( I-L )  1 00  *105*105"  ' 

100  IL=I+LQ 
"""  60  TO  HO' 

105  IL=L+IQ 

TTcTT=A  ( I  l  )  *cosx-A  Cm  X*s1nx 

A(  IM)=a(IU)*SINX+AUM)*COSX 

A(IL)=X  .  .  . ~ . “  . 

115  IF(MV-I)  120*125*120 
120  ILR=ILvj+I 
IMR=IMo+I 

X=R(ILk)*C05X"R(  IMR)*5ir^X 
R(IMR)=R(ILR)*S1NX+R(IMR)*C05X 

K(ILR)=X  . . 

125  CONTINUE 

X=2.0*A(LM)*SINCS 
Y=A(LL)*C0SX2+A(MM)*SINX2-X 
X=A(LL) *SINX2+A<MMT*tOSX2+X 

A ( LM ) = ( A ( LL ) -A ( MM ) ) *S1NC5+A ( |_M ) * ( C0SX2-SINX2 ) 
A(LL)=y 

A ( MM) -x 
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o  ojo  I  i  o  o  oj  Ojnqnoi 


TESTS  for  completion 


test  FOR  M  =  LAST  COLUMN  _ 

130_If«_^-Nj  i35»  140 » 135  _  _  _ 

135  M=MU 

GO  TO  qQ _ t _ 

TEST  FOR  L  =  SECOND  FROM  LAST  COLUMN 

14Q  IFa-(N-l))  145.  150  * 145 _ 

145  L=L+1 

_G0  TO  Jj5  _  _ 

150  IF(lNO-l)  lbO  . 155* 160 

155  IND=Q  _  _ _ _ __ 

60  "TO  50 


COMPARE  THRESHOLD  WITH  FINAL  NORM 
160~  I fTTHR-A NRMxT  lb 5 #165 * 45  ~  “ 

Z  - SORT  E 1 GENVALUE S“AW~O^En7E  CTORS- 


16b  IQ=-N 

_ DO  185  1  =  1 »N _ 

lQSlQ+ig 

_ LL= 1+ ( i*I-I)/2 _ 

JQ=N*( 1-2) 

DO  165  J~I » N _ 

«JQ— *JQ  +  |  j 

MM=JHj»J-J)/2 _  _ 

I F ( A ( LlV-ATmm ) )  170. 1 8b . 185 

170  X=A(LL) _  _ 

A  <  LL ) =A ( MM ) 

A(MM)=X _ 

IF(Mv-i)  175*185.175 
175  D0_180_  K=1#N 
ILR=IQ+K 

IMR=JO+K  _  _ _ 

X=R  (  ILr  ) 

H(ILR)=R(IMR1 
180  R(IMH)=X 
105  CONTINgE 
RETURN 
END 
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n  nn  n|  n|nnnn!nn!  n  n|n  ojn  r-  nn;n  n;r,njnor,  doom  do  n(n  n  nr,  n  n|n  o  n  o  n 


SUBROUTINE  NROOT 


PURPOSE 

COMPUTE  EIGENVALUES  AND  EIGENVECTORS  OF  A  REAL  NONSYMMETRIC 
A 1 R I X  OF  THE  FORM  B-IMvERSE  TIMES  A.  THIS  SUBROUTINE  IS 
'NORMALLY  CALLED  "BY  SUbROUTINE'CANOR  IN  PERFORMING  A 
CANONICAL  CORRELATION  ANALYSIS. 


USAcE 

CALL  NKOOT  ( M  *  A  *  ti#  XL  #  X  ) 

uL'suRTFnun'ar  ftirtoeters 

M  -  ORDER  OF  SQUARE  MATRICES  A*  B.  AND  X. 
a  -  INPUT  MATRIX  (M  X  .  " 

H  -  INPUT  MATRIX  (M  X  M )  • 

<L  -  OUTPUT  VECTOR  OF  LENGTH  M  CONTAINING  EIGENVALUES  OF 
B-INVERSE  TIMES  A. 

— a  - . output  matrix  tm  x~m)  'contain  rNG  ih'genvectoes  “column- 

M ISE  •  _  _ 


REMARKS 

nONST 


SUBROUTINES  AND' FUNCTION  SUBPROGRAMS  REQUIRED 

. 1 1 GEN  _ _ 

METHOD _ 

WEFER"  TO  W.  W.  COOLEY  AND  P,  R.  LQHNES*  ’MULTIVARIATE  PRO¬ 
CEDURES  FOR__THE  BEHAVIOrAJ.  SCIENCES’*  JOHN  wlLEY  AND  SONS* 
1962 *  CHAPTER  3. 


SUBROUTINE  NROOT  ( M*  A  *  6*  XL  *  X ) 

D I MENS  I ON  All )  » B ( 1 ) , XL ( 1 ) »X( 1 ) 


IF  A  OOUbLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED*  THE 
C  Tn  COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PRECISION 
STATEMENT  WHICH  FOLLOWS,  _  _ _ _ _ 

DOUBLE  PREC  IS  I  ON  A  *  B  *  XL  *  X,  SUMy _ _ _ _ _ _____ _ 

THE  C  MUST  ALSO  BE  REMOVED  FROM  DOUBLE  PRECISION  STATEMENTS 
"APPEARING  IN  OTHER  ROUTINES  USED  IN  CONJUNCTION  WITH  THIS 
ROUTINE.  _ _  _ _ _ _ 


FIGURE  A-IV-1.  (CONTINUED) 


213 


c 

C  THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 

C  CoNtATNOOUBLE  PRECISION  FORTRAN  FUNCTIONS r  "SORT  IN  STATEMENTS 

C  110  AND  175  MUST  BE  CHANGED  TO  DSQRT »  ABS  IN  STATEMENT  110 

C - RDTf  "Br  CHARED  "TO  "DAB'S . - : - 

C 


C 

C  COMPUTE  EIGENVALUES  AND  EIGENVECTORS  OF"B~ 

: 

K-l 

DO  100  J=2»M 
U=M*(J-ir  ” 

00  100  1  =  1. J 

L=1>1 ' 

K=K+1 

100  B(K)sB(L) 

C 

THE- MATRIX  B  IE  A  'PEAU  SYMMETRIC- MATRIX* 

C 

MV=0 

CALL  El GEN  (B.X.M.MV) 

T - 

C  FORM  RECIPROCALS  OF  SQUARE  ROOT  OF  EIGENVALUES.  THE  RESULTS 
T  “AR^PReMULTIPlIED  by  the  associated  e igEnvEc  t orST  ~  ““ 
c 

L=0  ' 

DO  110  J=1»M 

110  XL<J)=1.0/DSQRT(0ABS(B(L)  )) _ 

K=0 

DO  115  J=1.M _ 

DO  115  I=1»M 
K=K+1 

115  B(K)=X(K>*XL(J> 

C _ 

C  FORM  (B**(-l/2) )PRIME  *  A  *  <B**<-l/2>) 

C 

DO  120  I=l» M 

_ N2=0 _ _ 

DO  120  J=1#M  ~ 

N1=M*<I-1) 

L=M*(J-1)+I 
_ X ( L) =0 . 0 

DO  120  K=1 » M  -  ~ 

_ Nl=Nl+i _ 

N2=N2+x 

120  X ( L)  =X ( L ) +B ( N1 ) *A ( N2 ? _ 
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o  ojo 


LrO 

DO  13 0  J= 1»M _ 

00  130  1=1 »J 

Nl-l-M  __  _ _ _  _____  _ _ 

L=L+1 _ 

A(L)=Q,0  '  . . . 

_ DO  130  K=1>M _ 

N1=N1+m 

__  N2=N2+1 _ _ _ _ _ 

130  A(L)=A(U+X(N1)*B(N2) 

COMPuTF' EIGENVALUES ^W“EI6EWECT6RS  of H 

CALL  El  GEN  (A»X*MiWJ 
L=0 

DO  140  1=1  »M 
LrL+I 

1 4  o"  xltt)  =  a  mr 

^oMPuTe  The  normalized- eigenvectors 
00“  1 5CTT=i7M  — 

N2=0 

DO  lbO  J= 1 » M 
N1=I-M 
LrM* ( J-l ) +1 
A (L)=0,0 

DO  150  K=i»M  '  ' 

Nl=Ml  +  'i 

N2=N2+ 1  ~ 

150  A(L),=A(U+B(Ni)*X(N2) 

c=0 

KrO 

Do  lao  J=i»M  ' 

SUMV=0.0 

DO  170  1  =  1 » M 

L=L+1 

170  S0MV=5uMV+A<L)*A<U - 

175  SUMV=DsQRTtSUMV) 

DO  180  1=1 »M 
K~K  +  1 

T0O"  X(KJ=A  (K)/SUMV - - 

RETURN 
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_  kUaHOUTlNE  LOC(I. 

IX=I 

JX=J 

IF (MS- l > 10 » 20 #30 

10  IUX=N*{.JX-1)  +  IX _ 

GO  TO  36 

_ 20  IF  ( IX-JX )  22 »  24*24 _ 

22  IJX=IX+  (viX*JX-JX)/2 
GO  TO  36 

24  IJ X=JX+lI X  * IX-IX i/2 
GO  TO  36 
30  IJX=G 

IF ( I X-jX  >36*32*36 
"  "  32  TjX=I X 

_ 36  IJ=lOX _ _ _ 

RETURN 

END 


FIGURE  A-IV-1  .  COMPUTER  PROGRAM  FOR  FREQUENCY  ANALYSIS  OF 
SUBSTRUCTURE  (CONCLUDED) 
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/JOB 

1 


2 

3 
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G 
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8 

9 

10 

11 

12 

13 

14 

15 
11. 

17 

18 

ia 

50 

20 


FIGURE  A 


GO 

REAL  L,  LE 

READ! 5, 100)A,B, E,GXZ..GYZ 

HE  AO  (  5, 1UO)H,R,T,TE,RO 

READ!  5,  100)0011 ,UPT2,1)ELC,DELF 

IHD*0  . . . - - - 

THETA -A/ it 

L-B/T 
All  “A/  B 

G=H/(2.*T) 

ii-.8b'J*GYZ/E 

LE-B/TE  -  -  - - - 

C-GXZ/GYZ 

IF (THETA-. 35) 3, 3,  2 
WR I  TEUi ,  10 1 )  THETA 
I  NO 55 1 

I  F ( L-100. )4,5, 5 

WRITE(6, 102)1. . "" .  . . 

iriD-i 

IF(AR-3.  > 7,7,6 
UR  I  TRIG, 10 3  JAR 

I  ND-1  . 

IF(AR-.3)8,9,9 

WUIYE(6,103)AR  - 

I  ND-1 

IF(G-5, >10,11,11  - 

UKITE(6,1Q4)G 
I  ND-1 

IF(S-. 01)13, 13,12 

UR  I TE( 0,105)8  -  - 

INO-1 

IF(C-3. >15,15,14 
URITE(6,106)C 

INO-1  . 

I F ( C- . 3)10,17, 17 

WRITE! 6 ,106>C  - - 

I  ND-1 

!F(LE-1U0.)18,19,19  - —  - - 

WRITING, 107)  LE 

I  ND-1  - - - - 

I  F  (  I  NO )  2,G  ,  2G,  50 

WRITE(6,108>  -  . — - - 

GO  TO  1 

Cl -AR*  L*TH ETA/ G 
C2~-G/  (S*C*L**2) 

81  “1 .  +AR**4+4C  . 7*C2* ( 1 ,  ♦  C»Aft**2) 

S2-1.+40. 7*C2*(1,+C/AR**2)+1C55.*C2**2*C/AR**2 

S3=l .  ,  33*AR**2+38 .  « 14 €~) - 

f»4  =1 .  +3  .  *  AR**  2  +  iG  „  3»C2*  { 1 .  +7 .  *  C) 

FR2-1.+ .UQ2»S2*C1**2/S1 
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20 

rt 


22 

23 

24 

25 
30- 


FR»SQRT(FR2) 
IF(OPT2>20,20,23 
I F ( QPT1  >21,21,22 


GO  TO  30 
WCF-1./FR2 


GO  TO  30 

I FCOPT1 >24 ,24,25 
WCF-(DELF/DELC)**.5/FR2**.75 


GO  TO  30 

WCF  *L)ELF  /  C  OE  LOFR2 ) 


SISY€F*{  1. 


SIGXCF»C1.+.03G*(1.+.103*AR**2)*S2*C1/S4)*WCF 

S  hQXYO-2.5fi~*<<lrr^vQ^^AR^*2-KeJr*-i9-.-*&3/-S-2->y - 

1CC1.+.0258*AR**2)*C1+1G.  3*S4/S2) 


■0*7052  7  »<  1'  V  »  .  0 49*AR» »2  )  *S2»  Cl/S  3 
SIGH1“(C3+1.)/(C3-1. ) 


StGtFF-C ...... 

1C  1 . +9 . 62*AR**2+AR**4 ) ) *WCF 


TAUXZC»WCF 

TAUYZC-WCF 


- ITRTTFTFTTTtn - 

I FC0PT1)  33,33,  34 

33  - WRfTEt6,121? - 

GO  TO  35 

34  - WR1“TEtC,T22  )  - 

35  I FC0PT2  >32,32,30 

315 - WfUTEC  6,12  3) - 

32  WRI TEC  6, 124  > 

-  - tTRtTt{6-,t23>A,fi,t<,T,TE-,^- - 

WRITEC6,126)R0,EfGXZ,GYZ 
- 1 F  <tH>T2138,  38,37— - 

37  WRITEC6,127)DELC,DELF 

38  - tTRt  TEC  6, 128 ) - 

WRITEC  6, 129  > 

- WRt TECeviSfrHTtETA,  fc ,  AR ,G,  S - 

WRITEC  6,140) 

- wmteyiw-eyte - 

WRITEC  6,131) 

- WTC1  TEC  5 , 132  >  S  I  GY~CF - 

WRITEC6,133)SIGXCF 

- WRI  TECti,  134  ) S  I G XY C - - 

_ W R I TE C  6 , 1 3  5 ) S  j  GH 1 

WrItEC 6^ 137 >TAUXZC 

- WRITEC  6, 138) TAUYZC - 

WRI  TEC  6, 139  >  FR 

- GO  TO  5 - — - 


100  FORMATC8E10.3) 

101  - FORMAT ( *  C  »  *,E12. 5) 
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102 
10  5 


104 

105 
10G 


10  7 
108 


109 

120 


121 
122 
12  3 

124 

125 


12G 


127 


128 

129 


130 

131 

132 

133 

134 


135 

13b 

13  7 

138 

139 

140 


FORMAT ( ' L  » 
FORftAK  'Aft 


FORMAT( ! 0 
FORMAT ( 1 S 

FORMAT ( ' C 


'  ,E12.5) 


i 

i 

i 


, El  2 , 5  ) 
, El  2 . 5  ) 
, E12 . 5  ) 


FORMATCLE  -  '  ,  E 12 .5)  - - - 

FOIU1AT( 'CALCULATION  TERMINATED  FOR  THIS  CASE.' 
It  ‘  VARtABLEm  mttT-EO  AWYf  QUT-OF “RANQE,J~1 - 


FORMAT (2I1,2F7.4) 

FORMAT (  ,1‘//////,3X,'EFFE.CT  OF  CURVATURE  ON  HONEYCOMB', 

1'  SANDWICH  PANEL  STRESS')' 

FORMAT(/18X, 'WITH  RANDOM  EXOFTATION') 

F0RMAT(/1GX, 'WITH  SINUSOIDAL  EXCITATION') 

FORMATS  /l  OX,1  AND— DAM  PI  NG-Wft-CH  V-ARI  ES  W  I TH-  GUR VATU  REM - 

FURI1AT(//21X,  '  ****  INPUT  DATA****  '  ) 

FORMAT ( / 1  ARC  LENGTH  IS  '  ,F14 . 3,  2X,  ' LENGTH  FS  ',F19.3, 
l/'CORE  THICKNESS  IS  ' , F10 . 3, 2X, 'SKIN  THICKNESS  IS  ',F11.3, 


2/' EDGE  THICKNESS  IS  ',  F10 . 3,  2X,  CACTUS  TS”+,FI9 . 1 ) 

FORMAT (/'SKIN  DENSITY  IS  , Fl2 , 7, 2X, 'SKIN  MODULUS  IS  '#F13.0, 

It '  CORE  "  MOOU  LUS  ;  ~OXZ ~t  ~S  *~7  ~FS  .~0 X ,-'-0 Y  Z  I  S  J  7  F  22.0) - - - 

FORMAT(/' CURVED  PANEL  DAMPING  RATIO  IS',3X,F5.3 
I, / ' F LAT  PANEL  DAMPING  RATIO  tS*,5X/F5.3> 

FORMAT ( // 14X, ' ***NON-D I  MENS  I ONAL  PARAMETERS***' ) 

FORMAT ( /, 5X, 'THETA' ,9X, ' L' , 11X, ' AR'  , 10X, 1 G' ,  11X, 

1 '  S  '  ) 


FORMAT(7(2X/E10.3)) 

FORMAT (//2QX,'****5TR ESS  RATIOS****! ) 

FORMAT ( // ,  'STRESSES  AT  “PANEL  CERTER'T//,  " 

1' RAT  10  OF  SIGMA  Y  CURVED  TO  SIGMA  Y  FLAT  I S ' , 5X, E12 . 5  ) 

FORMAT( 'RATIO  OF  SIGMA  X  etffiVgO  TO-* . . . . . . 

I, 'SIGMA  X  FLAT  IS  ' ,4X,E12.5) 

FORMAT <  'RATTO  Of  SIGMA  Y  CURVED  TO  - 

1' SIGMA  X  CURVED  IS  ',2X,E12.5) 

FURMAT( 'RATIO  OF  SIGMA  Y  OUTf R~ TACf-TO — '"7 - — -  — 

1/26X, 'S IGMA  Y  INNER  FACE  IS  ',E12.5) 

FORMAT  (  //  'RATIO  OF  S  I  GMA  Y  OOft^EO-TO^/  /  9  X  ,  'STOMA  Y  L  >- - - -  - 

l'FLAT  AT  CENTER  OF  STRAIGHT  EDGE  IS  ',E12.5) 

FORMAT*  ' RATI-O  OF-TAO-XZ  CURVED  TO  TAU  XZ  -ft AT  IS  '>GX,  E12v5T 
FORMAT( 'RATIO  OF  TAU  YZ  CURVED  TO-  TAU  YZ  FLAT  IS  ’/6X/E12.5) 

FORMAT  (  //  '  THE  FREQUENCY  RATIO  IS  J,£l-2;-5) . . . 

FORMAT ( /GX, ' C',11X, ' LE ' ) 

END 
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/DATA 


29".~ 

.25 

00 

29. 

.25 

01 

0. 

.001 

"25; - 

84. 

0. 

23. 

84. 

.002 

- t-tWO<K)~r 

.008 

17000000. 

.008 

18000-; — 
.045 

18000. 

.045 

— 9000. - 

.000453 

9000. 

.000453 

29 . 

17000000  * 

1  SO 00  • 

- 9000  . . . 

.25 

84. 

.  008 

.045 

.000453 

10 

0. 

0. 

29. 

23. 

17000000. 

18000. 

9000. 

.25 

84. 

.008 

.04  5 - 

—  -.0004  53 

11 

r-errrr 

.001 

.002 

/ETCH 
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EFFECT  OF  CURVATURE  ON  HONEYCOMB  SANDWICH  PANEL  STRESS 
WITH  RANDOM  EXCITATION 


•  -  . . rfftPtJT  -  DA YA"» - - 

ARC  LENGTH  IS  29,000  LENGTH  IS  23 i 000 

CURE  THICKNESS  IS  0.250  SKIN  THICKNESS  IS  0.008 

EDGE  THICKNESS  IS  0.045  RADIUS  IS  84.0 

SKIN  DENSITY  IS  0.0«0453(f  -SKI  N  . IS - 17060003.- 

CORE  MODULUS,  GXZ  IS  18000.  GYZ  IS  9000. 


***NON"U  IMENS  I  GNAL  PARAMETERS11'** 


theta  . — l- - - - -m - tr - sr - 

0.345E  00  0.288E  04  0.126E  01  0.156E  02  Q.471E-03 

C  LE 

0.200E  01  0  .SHE  03  . ~ . . . 


****STR ESS  RATIOS**** 


STRESSES  AT  PANEL  CENTER  . . 

RATIO  “OF  STGTT/r  Y'CtTRVTTD  TO  "S  t  GMA  Y  FIAT  TS - 0tT?OSO1!-Y>T 

RATIO  OF  SIGMA  X  CURVED  TO  SIGMA  X  FLAT  IS  0.4535GE  00 

RATIO  OF  SIGMA  Y  CURVED  TO  SIGMA  X  CURVED  IS  0.1723KE  01 

RATIO  OF  SIGMA  Y  OUTER  FACE  TO 

SIGMA  Y  INNER  FACETS  0.18530E  01 


RATIO  OF  SIGMA  Y  CURVED  TO 

SIGMA  Y  FLAT  AT  CENTER  OF  S  TR  A  TGttTEPfr£"-TS — 0~.2fr80  8  E  -01 
RATIO  OF  TAU  XZ  CURVED  TO  TAU  XZ  FLAT  IS  0.30106E  00 

RATIO  OF  TAU  YZ  CURVED  TO  TAtT  YZ  FtAT— IS - OtSOTOSIt  00 


THE  FREQUENCY  RATIO  IS  0.22202E  01 
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EFFKCT  OF  CURVATURE  OR  HONEYCOMB  SANDWICH  PANEL  STRESS 
WITH  SINUSOIDAL  EXCITATION 


. — - Trwwv~It1''PU  T  ^AT/T*™-™ - 

ARC  tENuTH  IS  29.0<H>  LENGTH  IS  23.000 

CORE  THICKNESS  IS  0,2r>0  SKIN  THICKNESS  IS  0.00  8 

EDGE  THICKNESS  IS  0.045  RADIUS  IS  84.0 

SKTTTITENSTTT IS - D .  T)  OtTASTO — SK  TTP  rtOWCDS-TS - TTtrOOtmtT, 

CORE  MODULUS,  GXZ  IS  18000.  GY Z  IS  9000. 


*** NON* DIMENSIONAL  PARAMETERS*** 

- ITHETA - t - a-r - «— - rS - 

0.345E  00  0.288E  04  0.126E  01  0.156E  02  0.471E-03 

C  LE 

0 .200E  01  0.511E  03 


****STRESS  RATIOS**** 


STRESSES  AT  PANEL  CENTER 


CURVE U  TO"ij-K}ttA  RtAT — F 


RATIO  OF  SIGMA  X  CURVED  TO  SIGMA  X  FLAT  IS 
RATIO  OF  SIGMA  Y  CURVED  TO  SIGMA  X  CURVED  IS 
RATIO  OF  SIGMA  Y  OUTER  FACE  TO 

SIGMA  Y  INNER  FACE  IS 


0*.33079E  00 
0.171J4E  01 

0.18530E  Oi 


RATIO  OF  SIGMA  Y  CURVED  TO 

SIGMA  Y  FLAT  AT  CENTER  Of  STRAIGHT -EDGE-fS — OvE0038E  01 
RATIO  OF  TAU  XZ  CURVED  TO  TAU  XZ  FLAT  IS  0.20177E  00 
RATIO  OF  TAU  YZ  CURVED  TO  TAtT  YZ  F  LAT  TS - 0T2OT77E  00 


THE  FREQUENCY  RATIO  IS  0.22262E  01 


FIGURE  A-IV-2,  COMPUTER  PROGRAM  FOR  CURVATURE  EFFECTS  ON  STRESS  IN 
HONEYCOMB  PANELS  (CONCLUDED) 
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APPENDIX  V 


EMPIRICALLY  DETERMINED  NATURAL  FREQUENCIES 


A .  Ski n-Stringer  Plating 

To  obtain  an  esfimato  of  the  fundamental  frequency  of  the  skin-stringer  plating,  data  from 
the  sine-sweep  tests  were  used  to  establish  an  empirically  derived  frequency  equation. 

The  procedure  used  in  developing  the  frequency  equation  was  as  follows: 

1  .  The  fundamental  frequency  of  each  configuration  was  adjusted  by  dividing  by  the 
aspect  ratio  function,  1  +  (b/a)A 

2.  The  parameter  t/a^,  determined  from  Equation  (24b),  wos  then  computed  for  each 
configuration . 

3.  The  adjusted  fundamental  frequency  was  plotted  versus  the  parameter  t/a^. 

4.  A  "least-squares"  fit  of  the  data  war  determined  to  obtain  a  regression  line 
equation , 

Figure  A-V-l  is  a  plot  of  the  equation  resulting  from  the  correlation  of  test  data  and 
theory.  Clamped  and  simply  supported  edge  conditions  are  shown  in  dotted  lines  to  afford 
a  comparison  of  the  edge  fixities. 


B .  Honeycomb  Sandwich 

The  fundamental  frequency  data  from  the  tests  involving  the  honeycomb  sandwich  test 
panels  were  also  correlated  as  described  in  the  section  above.  The  parameter  (h+t)/a^  was 
plotted  versus  the  adjusted  frequencies. 

Figure  A-V-2  is  a  plot  of  the  "least-squares"  fit  of  data. 


ASPECT  RATIO  b/a 


APPENDIX  VI 

SKIN-STRINGER  TEST  SPECIMEN  DETAILS 


Details  of  the  skin-stringer  test  specimens  are  shown  in  Figure  A-VI-1  and  Table  AVI-T. 
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TABLE  A-VI-I 

BILL  OF  MATERIALS  FOR  SKIN  STRINGER  TEST  SPECIMENS 


Desig¬ 

nation 

Quan¬ 

tity 

No . 
Bays 

Dimensions* 

Remarks 

L 

W 

a 

b 

t 

t 

r 

b/a 

STR-1 

2 

9 

30 

20 

6.00 

12.00 

.032 

.040 

2.0 

lipped  stringer 

STR-2 

2 

9 

30 

20 

6.00 

12.00 

.032 

.040 

2.0 

STR-3 

2 

9 

30 

20 

6.00 

12.00 

.032 

.040 

2.0 

STR-4 

2 

12 

30 

20 

4.00 

8.00 

.032 

.050 

2.0 

STR-5 

2 

9 

30 

20 

5.00 

10.00 

.050 

1  .071 

2.0 

STR-6 

2 

9 

30 

20 

9.00 

18.00 

.071 

.071 

2.0 

STR-7 

2 

9 

30 

20 

5.00 

10.00 

.020 

.027 

2.0 

STR-8 

2 

9 

30 

20 

6.00 

12.00 

.100 

.125 

2.0 

STR-9 

2 

9 

30 

20 

9.00 

18.00 

.040 

.040 

2.0 

STR-10 

2 

9 

30 

20 

6.00 

9.00 

.032 

.040 

1.5 

STR-ll 

2 

9 

30 

20 

8.00 

16.00 

.032 

.032 

2.0 

STR-12 

2 

9 

30 

20 

5.00 

10.00 

.025 

.032 

2.0 

STR-13 

2 

12 

30 

20 

4.03 

8.00 

.020 

.032 

2.0 

STR-14 

2 

9 

30 

20 

9.03 

18.00 

.050 

.050 

2.0 

STR-1 5 

2 

12 

30 

20 

3.03 

6.00 

.025 

.040 

2.0 

STR-1 6 

2 

9 

30 

20 

5.03 

10.00 

.040 

.050 

2.0 

STR-1 7 

2 

9 

30 

20 

7.00 

14.00" 

.050 

.063 

2.0 

STR-1 8 

2 

12 

30 

20 

4. CO 

8.00 

.040 

.063 

2.0 

STR-19 

2  1 

9 

1  30  j 

20 

6. CO 

12.00 

.063 

.080  I 

2.0 

' 

STR-20 

2 

9 

30 

20  1 

1 10. CO 

20.00 

.090 

.090 

2.0 

1 

STR-2 1 

2 

9 

30 

20 

6.00 

12.00 

.032 

.040 

2.0 

same  as  STR-3 

STR-22 

2 

9 

30 

20 

5-00 

10.00 

.051 

.070 

2.0 

same  as  STR-5 

STR-23 

2 

12 

30 

20 

6.00 

9.00 

.032 

.040 

1.5 

redesign  of  STR-10 

STR-24 

2 

9 

30 

20 

9.  CO 

18.00 

.040 

.040 

2.0 

STR-25 

2 

12 

30 

20 

4.00 

6.00 

.040 

.063 

1.5 

STR-26 

2 

9 

30 

20 

6.00 

9.00 

.063 

.080 

1.5 

STR-27 

2 

12 

30 

20 

4.00 

12.00 

.040 

.063 

3.0 

STR-2 8 

2 

9 

30 

20 

6.00 

18.00 

.063 

.080 

3.0 

STR-29 

2 

9 

30 

20 

5.00 

10.00 

.032 

.050 

2.0 

STR-30 

2 

9 

30 

20 

8.00 

16.00 

.050 

.064 

2.0 

Letters  representing  test  specimen  dimensions  refer  to 
Figure  AVI-1,  Skin-Stringer  Test  Specimen  Drawing 
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Errata:  AFFDL-TR-67-156 

"Refinement  of  Sonic  Fatigue  Structural  Design  Criteria" 

P  one  Correction 

xiii  g  hj/t  =  h/(2t) 

xiv  L  b/t 


p  mass  per  unit  area  of  a  plate 

mass  per  unit  volume  -  Section  F 

The  (2,1)  element  in  the  matrix  of  Equation  (56b)  should  read 

(30[(6A)5  “  (S/b)2j) 

The  line  following  Equation  (128)  should  read,  "where  it  has 
been  assumed  that  tg  -  tg  =  t." 

Changes  in  the  definitions  in  Equation  (1 29)  should  be  L  =  b/t, 
9  “  h,/t. 

Following  Table  III  the  data  should  read 
L  =  b/t  =  1031  .25 
g  =  hj/t  =  1 1 .625 

In  Equation  (146)  the  "L"  in  the  numerator  should  be  "L  ". 

In  Equations  (156)  and  (157)  the  subscript  "2"  should  be  deleted 
from  t. 

In  Table  XVI 

change  "a"  to  "b" 
change  "b"  to  "a" 
delete  the  subscript  "2"  on  t 

In  Table  XiX  the  definition  on  the  second  line  should  be 

ayc,z=hj 

yc,z  hj  cenj.er 

Equation  (120)  should  read 


(a  +  bZ) 
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Correct  ion 

The  first  sentence  in  the  next  paragraph  should  read: 

"The  required  core  density  (Ibs/in  )  can  . . 

In  Equation  (178)  the  exponent  on  the  right  hand  side  should  be 
-3/4. 

In  Equation  (179)  the  entire  right  hand  side  should  be  raised  to 
the  power  -1  . 

In  Equation  (180)  the  exponent  on  the  bracketed  term  should  be 
-3/4. 

In  Equation  (181)  the  bracketed  term  should  be  raised  to  the 
power  -1 . 
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